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ABSTRACT 


RUMMLER, DONALD ROBERT. Creep-Rupture Data; Analysis - 
Engineering Application of Regression TecKniques (Under 
the direction of HAYNE PALMOUR' III,). 

The creep and rupture behavior of materials can 
control the design of structures which operate at elevated 
temperatures. In lieu^of an adequate fundamental under- 
standing, current design practice -makes use of a variety 
of empirical techniques to predict creep behavior. 

The results of investigations to apply regression 
techniques to the development of methodology for 
creep-rupture data analysis are presented. Regression 
analysis techniques are applied to the explicit 
description of the creep behavior of materials for space 
shuttle thermal protection systems. A regression analysis 
technique is then compared to five parametric methods for 
analyzing three simulated and twenty real data sets. 
Finally, a computer program for the efficient evaluation 
of creep- rupture data with five parametric methods- is 
presented. 
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GfENEHAL INTRODUCTION - 

The creep— rupture behavior of materials can and does 
control the design of many structural components. Designers 
and analysts in the nuclear power generation, aerospace 
turbine, and chemical processing industries, for example, are 
required to design structural components which must operate 
reliably for periods- up to forty years in complex, high 
temperature environments. Unfortunately, the current state of 
our understanding of the creep process does not allow the use 
of "first principles" for sizing components and predicting 
their service behavior. Consequently, the creep— rupture design 
techniques used today can at best be called "enlightened - 
empiricism." There is no generally accepted method of analysis 
for the prediction of creep-rupture behavior. In fact, a 
method which work-s well for one material very often will not 
work well for a different material. 

The purpose of the investigations reported herein was to 
explore the application of regression analysis techniques to 
the analysis of creep-rupture data of interest in aerospace 
applications. They constitute a part of a continuing effort, 
begun in 1970, to provide the materials related methodology 
necessary to, design efficient aerospace vehicles. 

The first paper deals with the application of regression 
analysis to the creep of space shuttle materials. Regression 





techniques are used as a tool (1) to assess the effects of 
sheet thickness and oxygen partial pressure on the steady-state 
creep behavior, (2) to analytically describe the low creep- 
strain behavior, and (3) to assess the effects of data scatter 
for materials where data are limited. 

The third paper describes the development and use of a 
computer program for parametric analysis of creep rupture data. 
The program includes provisions for the analysis of five 
different parameter methods* Sample problems to aid the user 

i * 

in setting- Up a problem are presented. 
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APPLICATION OP REGRESSION ANALYSIS TO CREEP OF 


SPACE SHUTTLE MATERIALS 
Donald R. Rummler 

NASA Langley Research Center 
Hampton, Virginia 

ABSTRACT 

Regression analysis techniques were used to assess the effects of sheet 
thickness and oxygen partial pressure and to develop constitutive creep equa- 
tions. Application of prediction intervals is emphasized. 

1 SYNOPSIS 

Metallic heat shields for Space Shuttle thermal protection systems must 
operate for many flight cycles at high temperatures in low-pressure air and 
use thin-gage (<0.65 mm) sheet. Available creep data for thin sheet under 
those conditions are inadequate. To assess the effects of oxygen partial 
pressure and sheet thickness on' creep behavior and to develop constitutive 
•creep equations for small- sets of data, regression techniques are applied and 
discussed. 

2 SYMBOLS 

G = creep stWin ' 

t ’= time, hours 
th - sheet thickness, mm 

T = temperature, -K 
2 

a = stress, MN/m 

X, y, z, D, ^ = dummy variables 



5 INTRODUCTION 


Recent Space Shuttle technology research and development studies ((l) 
and (2)) have indicated that the creep behavior of hi^-ten^erature alloys 
may control the design and reusability of metallic heat shields for radiative 
thermal protection systems (TPS). The heat shields function as lightly loaded 
aerodynamic surfaces, and they must efficiently utilize thin-gage sheet to 
avoid weight penalties. Loads are applied at hi^ temperature, when the local 
partial pressure of oxygen is low. In general, creep strains must be liroited 
to less than 0.005 to avoid excessive panel deflections. 

The creep data which exist for candidate superalloys are for steady- state 
creep tests run on relatively thick specimens at atmospheric pressure. These 
data are presented as time to a given strain level for various combinations 
of stress and temperature (see, for example. Refs. (3) and (k)). Attempts to 
use this type of data to predict the cyclic creep deformation of simple tensile 
specimens or for the preliminary design of heat shields underestimated the 
experimental creep strains by as much as a factor of 10 ((l) and (2)). These 
predictions typically utilized one of the parameter methods (5) combined with 
a life fraction approach to sum the cyclically accumulated strains. This 
failure to predict the experimental creep strains could be the. result of one 
or both of the following: 

(l) The data upon which calculations were based were for the creep of 
relatively thick specimens at atmospheric pressure, and may not be applicable 
to thin specimens at low pressure. 


* 


References are given in Appendix 1. 
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(2) No analytic expression was available which could account for both 
the nonlinear primary and linear secondary creep stages. 

The purpose of this paper is to present the results of an investigation 
to determine the applicability of regression analysis techniques to predict 
creep behavior when data are limited. Three applications of regression tech- 
niques which address the aforementioned shuttle TPS creep problems are dis- 
cussed. Regression techniques are used as a tool (1) to assess the effects 
of sheet thickness and oxygen partial pressure on steady- state creep behavior, 
(2) to analytically describe the low creep strain behavior, and (3) to assess 
the effects of data scatter for materials when data are limited. 

k ANALYSIS PROCEDURES 

)+. 1 Development 

To evaluate trends in creep data and to predict creep behavior, explicit 
expressions for -the mean and the expected upper and lower bounds for creep' 
strain data as a function of stress, temperature, and time were desired. 

Little information is available about the form of these expressions for the 
candidate materials at low levels of creep' strain. Consequently, two computer 
programs were written and applied to develop the desired expressions. Both 
programs utilize standard linear regression techniques (6). One program was 
of the form: 

w = b^ + bj^ p . (1) 

where w = log (stress) 

p = log (time) 


> 
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This program was used to generate coefficients, mean value estimates, and 
95 percent prediction intervals for data at specific values of strain and 
temperature. 

The second program was used to develop models for creep strain' as a 
function stress, temperature, and time. For this multiple regression program-, 
the equation form assumed was: 


f(y) = 


= g 


4* 


^ 1^1 


+ c 


A 


2 

V2 


^ V2 


+ c, 




2 ^ 


V5^ 




( 2 ) 


where y, x^, x^, and x^ are, respectively, functions of creep strain, 
stress, temperature, and time. 

Provision for transformation of y, x^, x^, and x^ was included in the 
program. The transformations, which included many of those found useful for 
analysis of creep data ( 7 ) were as follows: 


The prediction interval ( 6 ) is used to msike a statonent about the antic- 
ipated value of the dependent variable (y) for a future single observation at 
a specific value of the independent variable (x) or variables (x- , x^, Xj^ . . . 
for example, y will be between 2 and 6 for 95 percent of all future‘s single 
observations taken at x = 5 . The more familiar confidence interval, on the 
other hand, is used to make statements about the true mean value of yj for 
example, there is a 95 -percent probability that the true mean value of y at 
X = 5 is between 5 and 5- The prediction interval limits are wider since 
these include both the sampling errors and the uncertainties in estimating the 
mean value of y. 


Transformation 
Code (TC^.) 

0 

1 

2 ' 

5 
4 

■ 5 

6 

7 

8 ' 


Transformation 
(0 < i < 5) 


X. 

1 


= z. 


X. = 


log (z^) 


X. = 1/z 


X. = 

Is 


1/^i 

log (l/z ) 


X. 

1 


X. 

1 


Ijl’ (z^) 


(^i) 


1/2 


X. 

1 


z. + 1. 0 
1 


X. 

1 


X. 

■1 


= log (z. + 


(^i) 


^ 1/3 


1 . 0 ) 


where- the z^ are specific values of stress, temperature, or time. Similar 
functional transformations (y .= f(.D))- were used for strain. Each transforma- 
tion combination was assigned a .four^digit transformation n\miber where- the 
digits are the transformation code values for y, x^^, x^, and Xy respectively. 

Thus transformation 1025 -used the following -transformations: 


y = log b = I6g"‘ ('€ j ' 


Xg = l/Zg = 1/T, 

X, = = (t)V= 


Creep data sets usually include .a wide range of times, typically three 
orders of magnitude, -whereas, the ranges for creep strain, stress, and tempera- 
ture are seldom in excess of one .order of magnitude. Early analysis of 
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miltiple regression computer runs revealed that the comhination of the wide 

range in the variables associated with creep data sets and eq^uation forms 

2 

tfhich include terms that can be highly colinear, such as x and x , led to 
ill-conditioned normal equations which were subject to significant round-off 
errors during a matrix inversion operation. In order to minimize these 
errors, the data were scaled from 1 to 10 after transformation of the primary 
variables (y, x^, x^, x^) as follows: 

^i ” ~ ^mlij^t^max 7 ^min) ^ ^ 

^ij “9.0 ^x^j - x^ m±n)j^i max ^ ^i min^ ^ ^ 

where v . and y are the minimvun and maximum values of the transformed 
•'min •'max 

strain. The x. . and x. have similar definitions as they apply to 

1 min 1 max 

the transformed values of stress, temperature, and time. 

After transforming and scaling the primary variables. Equation (2) was 
expanded and new independent variables, defined as follows, were introduced: 

If 

/ 2 2 2\ , -u / 2 2 \ 

y - aj^a^a^ ®'i^2^3 V1^2^5/ 

' J=1 

This procedure results in an equation with 2y terms having linear coefficients 

(0j). ■ ' 

Some values of were set equal to zero so that, in Equation (3)^ the 

order (degree of interaction) for the number of terms in the regression 
analysis could be reduced as follows: 
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Order 

(Allowed term types] 

25 

4-th 

( X, X x^ and 
\ 1 m n 

2 2\ 

17 

3rd 

'(x,x X and 
\ 1 ra n 

x^x ) 
1 m/ 

10 

2nd 

(=<,* ) 
\ 1 m/ 


4 

1st 




(Hote that the reduced form can no longer be factored back to Equation (2).) 

4 . 2 Application 

To perform a multiple 'regression analysis using Equation (5)^ the order 
of the equation (k value) was selected first. Next, the transformations to 
be used on the primary variables were selected. Each observation of the data 
set was transformed, then scaled. The transformed and scaled values for 
strain, stress, temperature, and time were then used to generate values for 
the additional variables in Equation (5) . This data set was then used in 
the regression analysis. The mean values of creep strain were calculated 
from the coefficients derived during a multiple regression analysis. Explicit 
functions for the upper and lower bounds (95 percent prediction intervals) 
were calculated by treating either the upper or lower prediction limit calcu- 
lated for each observed value of strain during the initial regression as 
another set of observed strain values^ two additional regression analyses 
provided the desired coefficients. The residual mean square (EMS) for the 
prediction interval "data" sets were always extremely small (w 10 times that 
of the original data set analysis). This suggests that the errors involved in 
these approximations for the original prediction intervals were not large. 
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After a regression analysis was performed, all variables and residuals 

* t 

were descaled and back- transformed. Several quasi- statistical parameters were 
then calculated to aid model development and "best- equation” selection. These 
parameters are described as they are introduced. 

5 RESULTS MD DISCUSSIOR 

The following examples illustrate how regression techniques were applied 
to three areas of creep behavior which are of interest in Space Shuttle TPS 
creep studies. These areas are typical of those which can, occur during the 
preliminary design phases of any program when extensive creep data are not 
available. 

5.1 Use of Simple Regression (Equation (!))■ 

Haynes alloy H-I 88 is a cobalt base alloy which has excellent oxidation 
resistance and moderate elevated temperature strength. It is a candidate 
material for TPS application up to 1250 K. The creep data base consists 
primarily of the work reported in (t). This work includes creep tests bn 
H-I 88 sheet from 10 production heats and for thicknesses ranging from 0. 5 I to 
, 2.05 mm. All creep tests were run in air at standard pressure. 

Figure 1 presents the data at llh4 K at a strain level of 0.002. A 
regression analysis was perfomaed on the data set with sheet thickness 
< 0.84 mm. These data will be defined herein as the "standard data," against 
which data from future observations will be compared. The regression line 
and the 95-pe3:cent prediction interval for the standard data are also shown 
on the figure. The results shown in Figure I allow the following statements 
to be made: 
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(1) Ninety-five (95) percent of all future observations jaade under the 

lame test conditions are expected to fall within the prediction interval for 
iheet thicknesses between 0. 51 0* 1^" creep data from tests at 

lifferent test conditions generally fall outside of the prediction interval, 
>hen the new test conditions have probably changed the creep behavior of the 
laterial. 

(2) Most of the data for the > 0, 81j- Bim fall well within the prediction 

.nterval for the "standard data." Thus, the e = 0.002 creep strength of 
laynes alloy H-188 at 1144 -K is not significantly different for sheet thick- 
lesses from 0.51 to 2,05 This is in contrast to the results presented 

;.n (4)’ where creep rupture strengths of sheet <1.27 mm thick were lower than 
:hose for sheets >1.27 mm thick. 

The prediction interval and mean line from Figure 1 for the "standard 
lata" are shown in Figure 2. Also shown in Figure 2 are the results of 
ireep tests run in another laboratory on thin-gage H-188 at both standard and 
reduced pressures of air. The focus provided by the prediction interval indi- 
jates that the e = 0.002 creep strength of H-l88 for sheet thicknesses 
between O.51 and 0,64 mm both at standard atmospheric and reduced pressures 
■ras not significantly different from that previously established for 0.5I to 
D.84 mm sheet at standard atmospheric pressure. However, for thinner sheet 
(0.254 mm) at reduced pressure creep, strength was significantly higher as 
Indicated by the many test data points (open circles) above the prediction 
interval. Similar results were observed for other strain levels at ll44 K. 

The conclusions drawn from Figure 2 could have been reached with far 
fewer tests (as few as 2 or 5 for any of the test conditions shown). The use 
of prediction intervals data appears to be an efficient technique to ecplore 
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the eHects of "nonstandard" creep conditions and to compare creep data from 
different sources. This is particularly useful during the preliminary design 
phases of a program shen the consequences of "nonstandard" conditions, such 
as thin gage or low air pressure, must be assessed rapidly and maxlmu, us. of 
existing data base for thicker material at atmospheric air pressure Is necessary. 

5.2 Use of Multiple Regression (Equation (5)) 

To explore the effects of primary creep and various hardening rules, 
such as strain hardening, on the accumulation of cyclic creep strain, it is 
useful to have a constitutive relationship for steady- state creep strain. This 
is particularly true when the data base is limited and does not include a 
large number of test stresses and temperatures. 

The data set (8) for Rene sheet (solution treated at 1450 K and aged at 
1172 K) was selected to demonstrate the application of -multiple-regression 
techniques to develop a constitutive creep equation. ' Creep tests were con- 
ducted at 1005, 1089, and II72 K. Tests were not replicated. For this study, 

142 strain- time data points (observations) with strain levels from 0.0005 to 
0.005 were selected as input for the multiple regression analyses. 

In addition to a normal regression analysis, the program numerically 
solved the resulting equation to estimate the time (t^) required to reach 
each input strain level. To assure- compatibility with a strain- hardening 
cyclic- creep analysis, all equation forms which did not permit efficient 
solutions (less than 500 iterations) for all t^ were rejected. The program 
also rejected all equation forms which calculated either a negative strain or 
time. Early computer runs revealed that the multiple correlation coefficient 
square (e^) and the residual mean squared error (MSE), commonly used (7) to 
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rapidly evaluate a large number of equation alternatives were poor discrimi- 
nators for this data set and these variable transformations. The following 
parameters were determined from the descaled and back-transformed calculated 
values of strain and time: 

EMSE (strain mean squared error) 

e/TO (maximum calculated strain at t = 0.001 h) 

T/EO (maximum calculated time at e - 0,000001) 

AE (average strain error) 

ATP (average time error, percent) 

These parameters have recognizable consequences in the preliminary design 
sense and were considered useful discriminators for the selection of a "best" 
equation. Numerous variable transformations were evaluated in a single com- 
puter run. Typically, 200 different transformations were examined in a single 
600- second computer run. 

Analysis of several "best" equations during early computer runs indicated 
that the equations were often unstable near time = zero. This unstable behav- 
ior is illustrated In Figure 3 for typical values of stress and temperature. 
This failure to predict e = 0 at t = 0 was eliminated by assuming an 
unrecorded data point (e = 0.000001,, t = 0.001 h) for each creep test reported 
in (8). These assumed data points were added to the initial data set to yield 
the 167 data points and were Included in all further regressions. The dashed 
line in Figure 3 shows that a typical predicted creep curve using the addi- 
tional assumed points is reasonable, although the fit to' the original data 
(open circular symbols) is not as good. ■ _ , 

Even with the addition of the assumed data points, none of the variable 
transformations yielded a satisfactory prediction equation for the k = 27 
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version of Equation (3). ilie model was ^unstable •Hben projected on log- stress, 
log-time plots. At the lowest test temperature (IOO5 K) and short test times 
(as 10 h) these equation forms began to predict longer times for a particular 
level of creep strain as the stress was increased. For this particular data 
set, run i4-124 with k = 25 produced the "best" model equation. This run 
produced the lowest values of EMSE, AE, and ATP and computed e/TO < 0. 000001 
and T/EO <0.01 hr. The use of fewer terms in the model (k < 23) signifi- 
cantly increased the EMSE, AE, and ATP values calculated with the 
original 142 observations. This is illustrated in the following table; 


k 



23 

10 

7 

EMSE (x lo"^ 

6.71 

7.25 

9.14 

AE (X 10 ^) 

550 

622 

730 

ATP (0/0) 

55 

57 

74 


Thus for this data set, the inclusion of the higher order interaction terms in 
the model significantly improved the model’s ability to fit the data. 

The degree of fit typically provided by "best" model equation is illus- 
trated in Figure h for g = 0.002. The symbols are the data taken from (8), 
the solid lines are the mean stress and the 95-percent prediction interval 

calculated from a regression of log time on log stress using only those data 

* 

points shown for each temperature. The dashed lines are the mean stress 
values and the 95-pei'cent prediction intervals calculated by run hl2h, k = 23 
which included all of the 167 data points available in the data set. Agree- 
ment between the two calculated mean stress values is considered good. More 
importantly, however, this figure illustrates that the calculated 95-percent 
prediction intervals from run 4124, k = 23 are consistent with those obtained 
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from the linear regressions on the data for each temperatxire. This indicates 
that the model is prohahly as good as the data scatter warrant and that the 
consequences of this scatter can be adequately assessed in a steady- state 
creep analysis by utilizing the coefficients determined by run hl2.h to calcu- 
late mean creep strains and the coefficients determined for the lower bounds 
of the prediction interval shown in Figure 4 to calculate maximum creep 
strains. For instance^ a "best" model equation could be used to calculate 
cre^ strains at intermediate values of temperature to compare with other 
creep data obtained by other investigators. 

Figure 5 illustrates some typical mean creep curves calculated with the 
coefficients determined for the "best" equation. The shapes of these curves 
are consistent with those obtained by fairing through the original data 
points. More importantly, the curvilinear nature of the creep curves donon- 
strate that the model equation applies even when creep strain does not 
accumulate linearly as a function of time. Therefore, the model is function- 
ally capable of accounting for the effects of -primary stage creep in a strain- 
hardening analysis of cyclic creep. 

To further assess the applicability of the' regression analysis, the 
standard deviations for the average percentage time error for strain levels- 
0.001, 0.0015,, and 0.002 were calculated. These standard deviations were 
compared to similar results obtained from three optimized "c" value Larson- 
Miller analyses (5) of the data at these strain levels with the following 
results. 
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Comparison of Standard Deviation of -Percent Tme Error 


e 

Lar son- Miller 

Run 4124, k = 23 

0.001 

44.2 

20.4 

0. 0015 

46.7 

20.8 

0.002 

36.3 

33.5 


This comparison suggests that the "hest” regression equation, which includes 
all strain levels, predicts the observed creep behavior at least as well as 
the family of Lar.son-Miller curves which would be required to cover a similar 
range of strain levels. 

Multiple regression techniques can also be applied to fit "faired" data 
to estimate mean values for creep strain. This is illustrated in Figure 6. 
First, linear regressions of log time on log stress (Eq. (l)) were run on 
the original data set (8) for each level of strain and tenperature. The 
results of 'several of these regressions are shown as solid lines in the 
figure. Next, the mean times to a given level of strain were calculated from 
the regression equations of the solid lines. Finally, these calculated mean 
times and the appropriate values of creep strain, stress, and temperature 
were used as input data for a multiple regression apalysis (Eq. (3))* The 
dashed lines in Figure 6 were calculated from the results of a run 4-121, 
k = 27, using these calculated mean times as input data. The k = 27 version 
of Equation (3) was not unstable with the "faired" data set, whereas, as 
noted before, this version was unstable with the "raw" data. 

Often creep data are presented in the literature as families of faired 
curves for specific levels of strain and temperature. No individual creep 
curves are available for the material of interest. 
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As can be seen from this example, multiple regression techniques can be 
used to obtain a single equation which will coalesce families of curves. 
However, a prediction interval is no longer applicable because the calcula- 
tions are no longer based on scattered data. 

6 CONCLUDIHG REMARKS 

Frequently, creep data are limited during the preliminary design phases 
of a program such as the design of Space Shuttle thermal protection systems. 
The examples presented herein illustrate the applicability of regression 
techniques for (l) evaluating the effects of "nonstandard" creep conditions 
such as sheet thickness or low oxygen partial pressure on creep behavior and 
(2) developing analytical expressions to predict creep behavior from limited 
data. The use of prediction intervals to evaluate the design consequences 
of the data scatter has been discussed. 
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0 th > 0*84 mm (4) 

□ th < 0-84 mm (4) ' , 

lines: mean and bounds of 95 per cent prediction interval from 
linear regression of log t on log a for th < 0-84 mm 


Pig, 1. Creep strength of Haynes airoy H-188 
at 1144K, test pressure = 101 KPa, 
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Fig, 2. Effect of sheet thickness and test pressure 
on creep strength of H-I 88 at K, e = 0,022 







run 4024, ^:= 27 :a - 207 MN/m^, T= 1089 K 


o: 


experimental (8) 

— : calculated, raw data, 142 observations 
— : calculated, raw data + ‘zeros’, 167 observations 


Fig. 


3. Effect of 'zero’ data points on a typical 
calculated creep curve for Renfe 
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MN/m^ 



Symbols: experimental (8) 

: linear regression at each temperature, mean and 95 per cent 

prediction interval 

; run 4124, k = 23, mean and 95 per cent prediction interval 


Pig, 4, Comparison of experimental 
and calculated creep strength 
for Rene 4l, e = 0,002 


21 




MN/m2 



I ! 1 1 

10 , 100 1000 10000 


t hours 

; experimental means from linear regressions on data (8) for 

each T and e shown 

: calculated means from run H 21, = 27 


Pig, 6, Comparison of experimental and calculated 
mean creep strengths for Renfe 4,1 
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STRESS-RUPTURE DATA CORRELATION - 
GENERALIZED RECESSION ANALYSIS 
AN ALTERNATIVE TO PARAMETRIC METHODS 

By ‘ ‘ . 

« 

Donald R. Rummler 

ABSTRACT: The applicability of multiple regression analysis 

techniques to stress-rupture data correlation has been investi- 
gated, A generalized interacting variable (GIVAR) method of data 
correlation is proposed and evaluated. The GIVAR metnod is 
compared to six parameter methods of data correlation on' three 
sets of simulated data and tv/enty sets of real data. In all 
cases, the GIVAR method provided the best data correlation. 
Application of prediction Intervals and correlating variables in 
addition to temperature and stress is also discussed. 



INTRODUCTION 

Since 1952 when the first paper [1] introducing the concept of 
a time-temperature parameter (TTP) was published-, the need to ' 
correlate and extrapolate stress-rupture data has continued 
unabated. The importance of stress-rupture data analysis has led 
to a large number of papers which either propose new parametric 
approaches [2-5], offer detailed comparisons of analysis 
techniques [5~7]j and/or provide state-of-the-art surveys [8-12]., 
Although the development of some parametric methods can be 
related to ’creep behavior and fundamental prooesses, most 
parametric methods have been empirically derived. Most also make 
the assumption that there is a simple functional relationship 
between temperature and tlme-to-rupture which will yield a 
constant value of the parameter at a given level of applied- 
stress , Consequently, the selection of a particular- parameter to 
use for data analysis imposes rigid requirements on the nature- of 
the allowable interactions between time-to^rupture, , applied 

• A 

stress, and temperature. Methods for the selections of a 
particular parameter for the analysis of data, sets are given in 
the previously cited survey papers. The application ^f these, 
methods to real data sets is often difficult-. Often the analyst 

>r * 

is required to use data sets which -are inadequate , in terms of 

, . . f ' ‘ 

stress or temperature range to allow a clear, pe'lection ■ of the 
parametric method best suited for data -correlation. Data scatter 
further compounds the difficulty of selecting an analysis * 
technique and often forces the analyst to "smooth" or approximate 
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his data in order to conform reasonably to the functional 
requirements of a particular parametric representation. 

An attempt to overcome some of the difficulties has led to 
the concept of minimum commitment [7, 10, 13], This method 
(MCM) proposes the use of a general time— temperature functional 
relationship. The MCM method has recently been evaluated during 
an investigation concerned primarily with its extrapolative 
characteristics [ 7 ]. Although the MCM shov;ed promise during the 
®^9-luatlon, its clear superiority over other forms of parametric 
analysis was not demonstrated. In addition, in its present form, 
the MCM does not provide the analyst with an explicit form of 
parametric representation directly noi' is it completely general 
in the allov/ed functional interactions between the primary 
variables of time-to-rupture , stress, and temperature. 

The empirical nature of the data analysis techniques current- 
ly available is the direct result of the lack of understanding 
of the stress-rupture process particularly in complex engineer- 
ing alloys. Until better theoretical models of creep— rupture 
behavior are developed, the engineer or analyst is faced with 
the task of establishing- a functional relationship which will 
describe and correlate the data at hand. Regression analysis 
has been found to be a useful tool for the analysis of multi- 
factor data particularly when the physical factors which control 
the response to be predicted are understood only in general 
terms. Such is currently the case in the analysis of stress- 
rupture data. 
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The purpose of this paper is to present the results of an 
investigation to determine the applicability of multiple 
regression analysis techniques to stress-rupture data 
correlation. The particular regression techniques developed are 
first compared to several parametric methods using both simulated 
and real stress-rupture data sets. The potential of the 
developed regression techniques is further explored by subjecting 
a large number of real data sets to a preliminary analysis 
designed to select the functional form of an equation to be .used 
for detailed analysis. These results are also compared to 
several parametric methods. 

DATA FOR ANALYSIS 


Both simulated and real data sets were used to assess the 
capabilities of multiple regression analysis techniques for 
stress-rupture data correlation. 

Simulated Data 

Simulated data sets were derived from data for Timken 35-15 
stainless steel taken from reference [8], These data were fitted 
by the method of least squares to transformations of the 
following parametric expressions: 

Larson-Miller 

T(C .+ log t^) = b^ + b^ log a 
Orr-Sherby-Dorn 

log t^-AH/2.3RT = Dq + b^^log a 
Rabotnov 

a(l + A t^) = b^ + bj_/T + 
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where 

R - universal gas constant 
t^ time to rupture 
T = temperature 
a = stress 

Oj AHi A, bi b ^3 ^2 - constants determined by least 

squares 

5he Larson-Miller [1] and the Orr-Sherby-Dorn [2] expressions are 
familiar time-temperature parameters which assume that the 
parameter (left side of equation) is constant for a given stress, 
^he parameter can be considered a temperature compensated time, 
The Rabotnov [12^- 14] expression is a time-stress parameter which 
assumes that the value of the parameter (left side of equation) 
is a constant for a given temperature. The parameter represents 
a time compensated stress. Although the Rabotnov expression was 
originally developed for correlation of creep data, its use for 
creep-rupture correlation has been suggested [12] as an 
alternative to TTP methods. 

The values of the constants determined by the regression 
analysis for each parametric expression were used with the 
experimental stress and temperature levels to calculate "exact" 
times for each simulated data set. The simulated data sets are 
referred to as L-M Exact, 0-S-D Exact, and RAB Exact, Additional 
details of the fitting procedures and tabulation of the real and 
simulated data are presented in Appendix A, 

Real Data 


All real data were taken from a recent evaluation of para- 



metric methods for extrapolation ,[7]« Careful attention was 
paid to the adequacy of the data in terms of range of stress and 
temperature exposure and long times to rupture. The data in- 
cluded a wide range of materials. The material types and number 
of observations in each data set are shown in Table 1, The data 
set numbering in reference [7] has been retained in this investi- 
gation. Tables 2 and 3 present the data for the two data sets 
(4 and 16) which are analysed in detail. Reference [7] lists 
the data for the other data sets analysed. 

ANALYSIS PROCEDURES 

The three types of analysis techniques used during this in- 
vestigation (1) parametric 3 (2) minimum commitment (MCM), and 
(3) Generalised Iteracting Variables (GIVAR) are discussed in 
this section 
Parametric Methods 

A number of different parametric techniques have been 
suggested for correlating stress-rupture data. The equation 
forms used for multiple regression analysis of the parametric 
techniques selected for this investigation were as follows: 

Larson-Miller (L-M) 

Y = loe tr = + b^/Tj, + bjS/Tj, + b3 + 

bjsVTp + bgsVlj, 

'Qrr-Sherby-Dorn (0-S-D) 

y = log + b^/Tj^ + bgS + b^ + bj^S^ + b^S^ + bgS^ 

Mans on-Suc cop (M-S) 

Y = log = b^ + b^Tp + b 2 S + b^S^ + + b^S^ + bgS^ 
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Manson-Haferd CM-H) 


I = log + bjT^S + bjT^S^ + b,T„s3 + bjT„s3 + 

t’eV'' 

Rabotnov (RAB) 

Y = + b^/0Tp + b2/aTp + b^/aT^ + b^j/aTp + b^/aT^ 

where 

= time to rupture, hours 

S = log a 

o = applied stress, ksi 

Tp = temperature, 

= temperature, Kelvin 

Tp = temperature, Rankin 

= offset temperature “ 

b., T., a = constants estimated by method of least squares. 
Both the M-H and RAB techniques required the use of iterative, 
non-linear multiple regression techniques to estimate all of the 
constants . 

In all cases, some function of time to rupture was considered 
the dependent variable whose variance was minimized. High order 
polynomials which are functions of stress have often been used to 
correlate stress-rupture data [7, 8], Although a sufficiently 
high order polynomial can approximate any function, it can also 
result in unrealistic waviness, in plots -Of tbe dependent 
variable versus any one of the independent variables. For these 
reasons, the parametric model equation forms were also analyzed 
in functional forms which included only second or third order 
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polynomials in the stress function. 

In addition to estimating the required constants and predicted 
values of log time to rupture, the parametric analysis procedures 
produced the following summary values to- aid data correlation and 

parameter comparison: 

RMS = { gCOTR - , PTR)^ ,)l/2 

= / ^(OTR ~ PTR)^ .l/£ 

N - K - 1 '' 

DPAVG ^ t (PlMAX - PIMIN) 

N 

DPMAX - maximum value of PIMAX - PIMIN 
where 

OTR ^ observed log time to rupture 
PTR - predicted log time to rupture 
N - number of observations in data set 

K - number of constants in regression model 

PIMAX, PIMIN = upper and lower bounds of $^% prediction 
interval for each observation in a data set 
The root mean square (RMS) provides an overall comparison of 
data correlation including both random error and functional bias. 
It does not, however, reflect the increases in the regression 
standard deviation which can occur when high order polynomial 
terms are included in the model equation. The added high order 
terms may be highly correlated with the other independent 
variables already in the equation and consequently may not reduce 
the residual sum of squares enough to account for the loss in 
degrees of freedom [15]. For all regressions which used log time 
to rupture as the dependent variable the calculated value of STD 

8 


3 ^ 



is equivalent to the standard deviation of the regression. 

The average width (DPAVG) and the maximum width CDPMAX) of the 
95^ prediction interval are considered useful indicators of the 
expected scatter for a future observation taken from the same 
material under the same testing conditions. The prediction 
interval [16, 17] is used to make a statement about the expected 
value of the dependent variable (log time to rupture) for a 
single future observation at specific values of the independent 
variables (functions of stress and temperature). The prediction 
interval is wider than the more familiar confidence interval on 
the mean, since it includes both sampling errors and the 
uncertainties in estimating the mean value of the dependent 
variable . 

Minimum Commitment Method 

The minimum commitment method (MCM) of parametric analysis [7, 
10] was developed to minimize the dependence of the data analyst 
on the particular model equation forms of the generally used 
parameter methods. The MCM concept -is to utilize a parameter 

■ r 

model equation general enough to encompass most of the popular 
parameter methods. The parametric equation chosen has the form: 
log t (1 + AP) + P = G 
v/here 

t = time to rupture 
,A = constant 

P = function of temperature 
G = function of stress 
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The functions P and G are "station functions" which are defined 
by their values at selected levels of temperature and stress. 
Since it is not necessary for P and G to be explicitly expressed, 
there is no commitment on the part of the analyst to a particular 
parametric form, MEGA CManson-Ensign Generalized Analysis) is 
the computer program developed to Implement the MCM [13]. The 
particular version of MEGA used during this investigation 
utilized three stations of temperature to define P and three 
stations of stress to define 6. In addition, the firsu and 
second derivatives of the G function at the mid station were 
included in the analysis. The analysis, therefore. Involved the 
calculation of eight constants [7]. 

The parametric equation form which has been selected for the 
MCM does not readily lend itself to a least squares method of 
solution with log of time to rupture as the dependent variable. 
Consequently, the MEGA computer program in its current form does 
not yield least squares statistics such as the standard deviation 
of the solution (regression). The lack of appropriate statistics 
necessitated the use of RMS as the evaluator when comparing the 
MCM method to other methods of stress-rupture data correlation. 
Generalized Interacting Variables Method 

Development - The basic concept for the ^generalized 
Interacting vari ables (GIVAR) method of data correlation was- 
developed for the analysis and correlation of creep data [18], 
Simply stated, it is assumed that the functional relationship 
between the dependent variable and independent variables can be 
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described by a low order polynomial in each Independent variable. 
For stress-rupture data correlation, this concept leads to a 
model response equation of the general form: 

f(y) = g[(a^ + + c^X^)Ca2 + b2X2 + C2X2 + d2x|)] 

where y, X^, and X2 are respectively functions of time to 

rupture, temperature, and stress. Because complex Interactions 

between time, temperature, and stress are known to occur during 

the creep-rupture process, the model equation is completely 

general and allows all' interaction terms which result from the 

combination of the low order polynomials specified for each 

independent variable. Additional independent variables can be 

. 

readily introduced into the- general model form by the Inclusion 
of additional low order polynomials, 

f(y) = g[.(A) (B) (C) (D)] _ , . ' 

where A, B, C, D are low order polynomials of the independent 
correlating variables. '■ . ‘ - , 

The computer program to implement the GIVAR method includes 
provision for transformation of y and X^, For this 
investigation, the majority of data correlations were performed 
with the following transformations 
Variable Allowed Transformations 

y log t 

X^ T, 1/T, log T 

1/Q 

X^ cr, a log a 

where t, T, a are respectively time to rupture, temperature, and 
stress. After transformation of the primary variables, the model 
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equation form is expanded and new independent variables, defined 
as follows, are introduced to yield a response equation for a 
multiple regression analysis:. 


• ■ ■ , ' k 

y - a^a 2 a^ + b^Xj_ + ^ 2 X 2 + }3ib2X^X2 

The resulting model equation form for the multiple regression 
analysis is linear in the coefficients (<j>.) and is simply an 

^ u 

extension of equation forms which have been used to -determine 
optimum conditions in. multifactor environments [19]> for example, 
to determine .the conditions necessary to maximize the output of a 
chemical process. 

Application - To perform a GIVAR correlation of stress 
rupture data, the orders of the independent variable polynomials 
were selected and the general equation form expanded. A second 
order polynomial in temperature and a fifth order polynomial in 
stress were used for the majority of data correlations. When a 
transformation was selected, a sixth order polynomial in 
stress was used. Temperature and -stress interaction terms above 
third order (X^X 2 ) were deleted from the polynomial expansions. 


Next, the transformations of each prime variable which would be 
allowed were selected. The computer program, using these 
control inputs plus the original data set, then analyzed all 
combinations of the variable transforms and printed out summary 
results for each analysis. The variable transforms- which 
produced the lowes.t standard deviation of the regression were 
then resubmitted and the number of terms in the regression model 
was reduced using a technique known as>a t^ ^-directed search [ 153 . 
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When there are M potential variables in a regression model, 

M 

there are 2 possible regression equations. The t^ . directed 

search technique has been proposed as an alternative to stepwise 

regression techniques [l6] to reduce the number of variables in a 

regression model. The t^^ ^ directed search uses the ratio of 

each b^ to its standard error as follows: 

t, . = b . 

k,i 1 

where b^ and S(b^) are the values of the coefficient and the 
standard error for ith variable. Following a regression on the 
full model equat;lon, the variables in the full regression model 
are arranged in decreasing order of their t, . values. 

K , 1 

Successive regressions reduce the number of variables until a 
"basic set" is found. The program then analyzes all model 
equations which can be constructed including all of the basic set 
of variables plus all possible combinations of the previously 
dropped variables. The "best" equation is selected en the basis 
of the lowest standard deviation of the regression. 

Finally, the "best" reduced variable regression equation was 
analyzed in detail' to verify its adequacy. If the model was to 
be used for significance tests or if a statistical Interval such 
as the prediction interval were to be used, verification Included 
careful examination of residual plots [15j 16, 20] to assess 
departures from the assumptions of the linear regression model, 
REaULTS AND DISCUSSION 
Simulated Data 
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The purpose of the simulated data sets was to assess the 
functional capability of the GIVAR method and its associated 
computer program without the confusing Influence of the large 
scatter normally associated with stress-rupture data. 

The results of the simulated data set analyses are summarized 
in Table ^ which shows the calculated values of STD for each of 
the six methods of data correlation for the three simulated data 
sets. For each data set, the generalized Interacting variables 
method (GIVAR) produced the lowest value of STD. Of equal 
importance to the significantly better correlation was the fact 
that the GIVAR computer program selected the most correct of the 
prime variable transformations for the L-M and 0-S-D Exact data 
sets. The search quickly reduced the original nine term 

model equations to the correct three term equations. The value 
of STD calculated for these two cases is due primarily to 
rounding off the calculated exact times for these data sets. For 
the RAB Exact data, log t, log T, and log a were selected as the 
best prime variable transformations. In this case, the original 

r 

eleven term model equation was reduced to nine terms during the 
t^ 1 search. 

Table also illustrates the general futility of adding higher 
order polynomial terms to improve correlation for the restricted 
models. For the -four commonly used parameters, no significant 
improvement can be seen when expanding the model equation from 
four terms to seven terms (from a second order to a fifth order 
equation in stress), A similar lack of correlation Improvement 



has been reported on real data I5J, 

The correlations produced by the M-H and GIVAR methods for the 
RAB Exact data are shown in Pig. 1, The GIVAR method correlation 
is noticeably better than the M-H correlation. It is important 
to remember that in both analyses^ log time to rupture was the 
dependent variable and consequently, minimization of differences 
between observed and calculated times to rupture was the 
regression criteria. For these data, neither of the two methods 
shown had model equation forms which would exactly duplicate the 
governing equation- for the RAB Exact data generation. This is a 
comparable situation to most real data where correlation models 
seldom represent a material's behavior exactly. Since for most 
real data either correlation would probably be considered ‘ 
satisfactory, the calculation of a statistical interval such as 
the prediction interval to assess uncertainty about a future 
observation would be a natural extension of these correlations. 

The residuals of the M-H and GIVAR correlations for the RAB 
Exact data are presented in Pig. 2. The M— H residuals clearly 
exhibit curvature as a function of the predicted log time to 
rupture. The residuals are not randomly distributed with respect 
to the dependent variable (predicted log time to rupture). This 
type of behavior indicates that the regression model is 
inadequate and needs additional terms. What has happened is that 
the M-H model equation, even with a fifth order polynomial in 
stress, was functionally incapable of correctly approximating the 
Rabotnov expression which was used to generate these data. The 



random distribution of the GIVAR correlation which includes 
interaction terms does not suggest any functional inadequacy. An 
examination of the cumulative normal distribution of the 
residuals for the GIVAR correlation failed to indicate that the 
residuals were not normally distributed. Since the GIVAR 
correlation equation of these data does not appear to violate any 
of the basic regression assumptions, the calculation and use of a 
statistical Interval would be in order [l6]. 

Real Data 

The results of the GIVAR correlation on alloy 4 (a plain 
carbon steel) are presented in Fig. 3, As for all GIVAR 
correlations, log time to rupture was the dependent variable. 

The prime variable transformations selected by the computer 
program are shown. The original eleven term model equation was 
reduced to seven terms during the t^ ^ search. The GIVAR mean 
fit seems to satisfactorily correlate this complex behavior. The 
STD value of the GIVAR correlation for these data was 40 percent 
lower (0,103 versus 0.146) than a third order M-H model which was 
the best of the parameter models. 

To minimize the computer time, the 95^ prediction interval 
about each observation is normally calculated during the computer 
run which performs the regression on the model equation. The 
upper and lower bounds of the 95^ prediction are listed along 
with the calculated time to rupture. For these data, the 
calculated prediction interval called attention to a possible 
outlier, x,e,, an atypical observation. This data point is shown 
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with the filled symbol. Examination of the residual plot with 
respect to predicted log time to failure CPig» suggested that 
the residuals were randomly distributed, had a mean of zero, and 
exhibited constant variance with the single exception of the 
residual for the possible outlier. The cumulative normal 
distribution plot of these residuals (Fig, 5) also appeared 
normal with the exception of the single suspect data point. 
Although there are many schemes for outlier rejection [21, 22], 
the present purpose is to demonstrate that the prediction 
interval provided a useful tool for focusing attention on a 
possible outlier which may have otherwise been overlooked. For 
other data sets, the calculated prediction Interval has called 
attention to data transcription errors which had gone undetected 
because of large data scatter. It should be pointed out that the 
use of the prediction interval to provide a focus for possible 
outliers is not strictly .correct In the statistical sense. Its 
proper use is to make estimates of the bounds which can be 
expected from a single future observation from the same 
population. Dismissing the outlier for the moment, we can say 
that 95% of the time a future single observation will fall within 
the bounds shown in Pig, 3, The implications of this kind of 
statement for acceptance testing, quality control, or determining 
the significance of a new test variable, are obvious. 

Temperature and stress are usually considered the prime 
variables for stress rupture correlation. Some authors [5], 
however, have been able to improve correlation by the use of an 



additional variable such as elastic modulus to normalize stress. 
Table 5 summarizes the results of correlation analyses on alloy 
16 (a nickel base alloy) to evaluate the effect of additional 
variables. The listing Includes the analysis method, the prime 
variable transformations, and the calculated values for STD, 

DPAVG and DPMAX. The units of DPAVG and DPMAX are log (time to 
rupture, hours). For these data, the M~S and M-H methods were 
the best (lowest STD) of the parameter methods. However, the use 
of elastic modulus (E) to normalize stress did not significantly 
Improve the fit in either case. Using just temperature and 
stress, the GIVAR method resulted in a significantly lower value 
of STD than the best parameter method. When second order 
polynomial expressions for elastic modulus and ultimate tensile 
strength at the test temperature were incorporated into a 
generalized interacting model equation, a significant further 
correlation improvement was achieved. The significance of the 
better correlation provided by the GIVAR method is more easily 
appreciated when it is realized, that within the average 
prediction interval bounds, the predicted time to rupture varies 
by a factor of 3 for the best parameter method and by a factor of 
1.6 for the GIVAR method. For the maximum width of the 
prediction Intervals, these values are 4.5 and 1.8, respectively. 
It should be pointed out that the GIVAR model equation did not 
allow interactions to occur between .elastic modulus or ultimate 
tensile strength and temperature, since they are both highly 
correlated with temperature. In this 


case, the original 21 term 



search 


model equation was .reduced to 13 .terms during the t 

The best M-H and GIVAR correlations of the ’alloy l6 data are 
presented graphically in Pig, 6. The GIVAR fit is noticeably 
superior. Even with a fifth order polynomial in log stress, the 
M-H model equation appears to be functionally inadequate to 
correlate the complex behavior of alloy 16, This functional 
inadequacy is further demonstrated in Pig, 7 which presents the 
residuals as a function of the predicted log time to failure. 

The M-H residuals are not randomly distributed and definitely 
display a curvilinear tendency suggesting the need for 
interaction terms. The GIVAR residuals appear to be randomly 
distributed and do not suggest any inadequacies in the model 
equation form. The cumulative -normal distribution of the ■ 
residuals for the GIVAR solu-tion (not shown) did not reveal any 
gross departures from normalcy. Since none of the basic 
assumptions of the linear regression appear to have been 
violated, the making of significance statements or the 
calculation of statistical intervals for this solution would be 
in order. 

In order to further assess the generality of the GIVAR method, 
all of the data sets of reference L7] were correlated with the 
five parameter methods, the MCM method and the GIVAR method. The 
independent variables for these analyses were limited to 
functions of temperature and stress. For the parameter 
methods, second, third, and fifth order model equation forms were 
■ examined. The lowest RMS values for the five parameter methods, 

‘fs 


19 



MOM and GIVAR methods are tabulated In Table 6 and presented 
graphically in Pig. 8. RMS was selected as the basis of 
comparison in order to include the MCM analyses.. Additional 
details and other summary values for these analyses are presented 

r 

in Appendix B. . - . 

In Pig, 8, a range band is shown for the five parameter 

methods. The MCM and GIVAR method are shown with symbols.. ’ Por 

each of the twenty data sets analyzed, the GIVAR method produced 

the lowest value of RMS, The GIVAR method on the average 

• ♦ ^ 
porduced a 19^ lower RMS value than the MCM 'which was on the 

average the best of the other methods examined. Examination of 

Table 6 reveals that the GIVAR solution in several cases required 

less terms in the model equation than the best parameter model 

equation. The MEGA computer program used to implement the MCM 

required the determination of eight constants. Table 6 also 

shows that the Rabotnov method was in all cases the worst of the 

parametric methods. It should be pointed out, however, that a 

polynomial in 1/T was the only function of temperature 

Investigated and that other functions of temperature might 

provide better correlations. With the exception of the GIVAR 

method, none of the other methods consistently produced the 

lowest RMS value for all twenty alloys. The failure of any 

single method to be consistently superior was also observed in 

reference [7] where the primary emphasis was on the extrapolative 

characteristics of the various parametric methods with these sets 

of data. 



CONCLUSIONS 

An investigation has been made to assess the applicability 
of a generalized interacting variable (GIVAR) multiple regression 
analysis method for the correlation of stress-rupture data. The 
GIVAR method was compared to six other methods of stress-rupture 
data correlation on twenty sets of data. The following conclu-' 
slons are made from the analyses presented herein. 

1. For all data sets examined, the GIVAR method produced the 
best correlation (lov;est RMS value). 

2. It was shown that the GIVAR method has the functional 
generality to satisfy criteria necessary for the calculation of 
statistical intervals., 

3. The GIVAR method readily accepts the Inclusion of correla- 
ting variables in addition to stress and temperature. 

The prediction interval v;as shov/n to be useful for the detec- 
tion of possible data outliers. 



APPENDIX A 


Parametric Analysis to Establish Simulated Data Sets 

The purpose of simulated data sets was to evaluate the func- 
tional capabilities of the various correlation methods without 
the confounding influences of the large scatter normally assoc- 
iated with real data. Creep rupture data are seldom the result 
of a statistically designed experiment. The data are seldom 

t 

I, 

balanced in variable space. In addition, temperature and stress 
are often highly correlated. Because of testing economics, low 
stresses are usually associated with high temperatures and high 
stresses are usually associated v/lth low test temperatures. In 
order to include this type of Inbalance in the simulated data 
sets, the data for Timken 35-15 stainless steel [83 were fitted 
to a first order Larson-Miller and Orr-Sherby-Dorn expressions 
and to a second order Rabotnov expression by the method of least 
squares. The equation forms and the fitted coefficients were as 
follows : 

Larson-Miller 

(TxlO“‘^) (C + log t^) = + b^ log 0 

where T = test temperature, 

C = iteratively determined constant = 13 
t^ = time to rupture, hours 
b^ = 6.39038 
b^ = -0,9058it 
a - stress, psi 
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Orr-Sherby-Dorn 

log 1°S o 

^ 2.3RT 

where 

= time to rupture, hours 

AHp = apparent activation energy, iteratively calculated = 

58000 

R = universal gas constant = I .986 

T = temperature, K 

b^ = 46410 

b^ = -4,60029 
O' = stress, psl 
Rabotnov 

t®- ='b^ + b^/aT + b2/aT^ 
v/here 

t = time to rupture, hours 

a = constant iteratively determined = 0.3637 
b^ = -1.62434 
b^ = -2.44083 X 10 ^ 

o = stress, ksi 
T = temperature, °F 
b 2 = 4.88958 X 10 ^ ‘ 

The rupture times which were calculated for each of the three 
solution methods were substituted for the experimental times to 
rupture to form the "exact" simulated data sets. These calcula- 
ted times and the original data for the Timken 35-15 stainless 

I 

steel are presented in Table, 7. 
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APPENDIX B 


Supplementary Analysis of Correlation Methods 

The purpose of this appendix is to supplement the correla- 
tion method comparison presented in the main body of the paper 
on the twenty real sets of data. 

The results of the parametric correlations are summarized 
for the L-M, O-S-D^ M-S, M-H and RAB in tables 8 through 12, 
respectively. The tables present values of RMS, STD, DPAVG and 
DPMAX which were calculated for each level of polynomial model 
equation which was evaluated. For the L-M, 0-S-D, M-S, and M-H 
methods, second, third,, and fifth order expressions in stress 
required 4, 5, and 7 terms, respectively. The RAB method re- 
quired 3, 4, or 6 terms to develop' second, third, and fifth 
order expressions. Table 13 presents a summary of the GIVAR 
method for these twenty data sets. ' > 

Parametric methods 

In all cases for the L-M, 0-S-D,, M-S, arid M-H methods, a 
fifth order expression produced the lowest ■, value of RMS for a 
given alloy. In some cases, however, the high correlation of 
the power terms in stress resulted in ill-conditioned solutions 
which were not reliable (see Table 8, alloy l4, for example). 

Such was not the' case for the RA-B solutions (Table 12) where 
third order expansions (4 terms) of temperature fit better than 
fifth order in a number of cases (alloys 4, 6, 8, llA, IIB, 17A). 

The calculated values of STD, which for the L-M, O-S-D, M-S, 
and M-H methods were equivalent to the standard deviation of the 



regression, did not follow the trend of better correlation with 
increasing degree of stress polynomial. The increased STD values 
reflect the fact that added variables did not reduce the residual 
sum of squares enough to account for 'the loss in degrees of 
freedom. These cases included the following: 

Alloy Method ( 5 )' . 

1 L-M 


4 

L-M, 

0-S-D, 

M-S, 

M-H 

6 

L-M, 

0-S-D, 

M-S, 

M-H 

8 

L-M, 

M-S, M- 

-H 


llA 

L-M, 

0-S-D, 

M-S, 

M-H 

IIB 

L-M, 

0-S-D, 

M-S 


12 

L-M, 

0-S-D, 

M-S, 

M-H 


This behavior, larger values of STD with a higher order 
polynomial, was also exhibited for several of the alloys during 
the RAB method correlations (’Table 12), The poorer correlation 
provided by the higher order polynomials can be better 
appreciated when we recall that the units of DPAVG and DPMAX are 
log time. Taking the best parametric method correlation in terms 
of RMS for alloy 4 (Table 11), we see that the average predicted 
time within the 95^ prediction interval varies by a factor of 4.9 
for a seven term equation and by 4.5 for a five term equation. 

The comparable values for the maximum width of the prediction 
interval are 6,4 and 5.^. In this case the use of a 
fifth order expression has significantly degraded the 
correlation. In addition to providing more sensitivity to 


25 



changes in the "goodness” of correlation, the values of DPAVG and 
DPMAX as preliminary evaluators of ■ correlation have the feature 
of allowing all methods to be compared on an equal basis. Values 
of DPAVG and DPMAX can be backtransformed and averaged if 
necessary to accommodate different transforms of the dependent 
variable. They can thus provide the- analyst with a "feeling" for 
the scatter and uncertainty in the data and its correlation. 

It is beyond the scope of this paper to summarize the results 
of all of the analyses which were performed by the GIVAR method 
on the real data. Table 13 summarizes the "best" model equation 
results for each alloy. In most cases, the "best" equation was 
selected- after the examination of summary computer results for 
nine different model equation forms. Log T and 
transformations of temperature and stress were selected for 
several of the alloys (Table 13). These transformations are not 
suggested by any of the standard parametric methods. As 
expected, not only did the GIVAR method produce the lowest value 
of RMS for each of the alloys, but it also produced the lowest 
value of the other preliminary correlation evaluators STD, DPAVG, 
and DPMAX (Tables 8 through 13). 

It is rare that stress-rupture data have the replicated 
observations that' are necessary to provide an Internal estimate 
of data scatter. The data for alloy 13 [7] was such an 
exception. There were seventeen experimental conditions which 
were replicated. These replicated observations had an average 
standard deviation of 0,232 with a spread of from 0,024 to 0.476, 



in terms of log time. The best GIVAR correlation of these data 
(Table 13) had a standard deviation of 0,280 indicating that the 
fit was comparable to the data scatter. This value is somewhat 
lower than the best (M~H) parameter method STD of 0,293 (Table 
11 ) . 

Summary 

1. Higher order polynomial model equations do not always pro- 
vide the best correlations of stress-rupture data, 

2. The standard deviation of the regression (STD) is a better 
correlation evaluator than RMS. 

3. The average and maximum width of the 95^ prediction Interval 
(DPAVG and DPMAX) are sensitive preliminary evaluators for 
stress-rupture data correlations. 



ALLOY 


1 

2 

4- 

5 

6 

7 

8 

9 

llA, IIB 
12 

13 

14 

15 

16 

17A, 17B. 
18a, ‘I-SB 
19 


TABLE 1 — Real data ’sets examined. 


MTERIAL 

NUMBER OP 


OBSERyATIONS 

f / 

llOO-O' ALUMINUM’ ' 

64 

5454-0 ALUMINUM 

. 75 

PLAIN CARBON STEEL 

26 

lCr-lM 6 STEEL 

33 

ICr-lMo- 0.25 V STEEL ' 

26 

304 STAINLESS STEEL 

52 

30‘4 STAINLESS STEEL 

39 

316 STAINLESS STEEL 

38 

347 STAINLESS STEEL 

42,44 

A286 IRON-NICKEL 

24. 

INCO 625 IRON-NICKEL 

99 

INCO 718 NICKEL-BASE 

26 

RENE 41 NICKEL-BASE 

37 ' 

ASTROLOY ® NICKEL-BASE 

. 33 

UDIMET 500 NICKEL-BASE 

103,105 

L-605 COBALT-BASE 

100,104 

606 I-T 65 I ALUMINUM 

99 


TOTAL =20 DATA SETS 



TABLE 2 — Stress-rupture data for alloy ^ 


Temperature 

Stress 

ksi 

Time to Rupture 
Hours 

752 

40.3 

752 

752 

38.1 

1696 

75-2 

35.8 

3973 

752 

33.6 

6134 

752 

31.4 

10422 

752 

29.1 

20227 

842 

33.6 

65 

842 

31.4 

44 l 

842 

26.9 

1341 

C\J 

CO 

24.6 

3023 

842 

22.4 

3934 

842 

17.9 

12985 

842 

15.7 

18648 

842 

13.4 

34753- 

932 

22.4 

63 

932 

20 . 2 - 

247 

932 

17.9 

430 

932 

15.7 

1317 

932 

13.5 

2958 

932 

11-.2 

3202 

932 

9.0 

7558 

932 

6.7 

22707 

1022 

13.5 

43 

1022 

11,2 

' 142 

1022 

9 . 0 - 

‘4,96 

1022 

6.9 

1935 
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•"’table 3 — Stress rupture data for alloy 16 . 


Temperature 

Stress 

Time to Rupture 

Tensile Strength^^^ 

Elastic Modulus 

°P 

1 

ksl 

' Hours 

ksi 

10”^ psl 

1^100 

101.0 

■ 12.8 

150' 

25.80 

lAOO 

■ 86.0 

‘■59.0 

150 

25.80 

ii»oo 

* 80 . 0 .: 

176.6 

150 

25.80 

1400 

74.0 . 

■ •' 400.7 

150 

25.80 

i4'oo 

. 70.0... 

■ ' 577.0 

150 

25.80 

14 bo 

61.0 

2279.8 ' 

150 

25.80 

1400 

55.0 

4,063.2 

150 

25.80 

1500 

75.0 

30.5 

130 

25.05 

1500 

64.-0 

142.2 

130 

25.05 

1500 

, 56.0 

351.3 

130 

25.05 

1500 

■52.0 

712 v 0 

130 

25.05 

1500 

4'5.0- 

1228.3 

130 

25.05 

1500 

'39.0 

2227.4 

130 

25.05 

1500 

31.0 

• 4393.4 

130 

25.05 

1600 

64.0 

10.5 

lid 

24.50 

1600 

56.5 

28.8 

110 

24.50 

■1600 

46.5 

145.8 

no 

24.50 

1600 

41.0 

253.0 

no 

24.50 

1600 

37.0 

535.7 

lid 

24.50 

1600 

31.0 

888.0 

iio 

24.50 

16B0 

.24..- 5 

2899.7 

no 

24.50 

1600 ' 

19.0 

6331.0 

11-0 

24,50 

1700 

41.0 

11.5 

.80 

23 . 3.0 

IJOO 

33.5 

44-. 2 

•80 

23.30 

1700 , 

. , 29 . 0 ' 

120.9 

80 

■23.30 

1700 

24.0 

342.7 

80 

23.30 

1700 

21.0 

. 7'46.7 

80 

-23^30 

1700 

17.5 

1768.7 

80 

23.30 

1700 

14.5 

2838.7 

80 

23.30 

1800 

29.5 . 

6.1 

40 

22.15 

1800 

20.5 

49.3 

40 

22.15 

1800 

17.0 

174.0 

40 

. 22.15 

1800 

14.5 

340.7 

4o 

22.15 


•Si; 


(a) Estimated from reference [233 

(b) Prom reference [53 





IABJ.E 4— Comparison of STD values for simulated data. 



Data 

Set 

No, of 

equation 

terms 

(1 

L-M 

5 

7 

4 

0-S-D 

5 

7 

Parameter Methods 
M-S 

^ 5 7 

31 

0-S-D 

Exact 


.oil? 

.045 

.045 

-- 

- 

— 

-.092 

.095 

.036 . 


L-» 

' Exact 


— 

~ 

-- 

.048 

.047 

.047 

.052 

.054 

,055 . 


RAB 

Exact 


.117 

.118 

.119 

.145 

.149 

.150 

.083 

.084 

.085 . 



OenoraXlzed Interacting Variable^ 
OIVAR No. of varia- Prime variable 


H-H 

5 

7 3 

RAB 

4 5 


bles in “best" 
equation 

transforms t i 0 hs 

.023 

.024 .317 

.046 .003 

.00003 

3 

log t, 1/T, log 0 

,024 

.025 .202 

.052 .051 

.00002 

3 

log t, 1/T, log a 

.042 

.044 — 



.010 

9 

log t, log T* log 0 



TABLE 5 — Effect of additional variables on correlation 



alloy 16 

- Astroloj^ 



ANALYSIS 

METHOD 

PRIME 

VARIABLES 

STD 

PREDICTION 

AVERAGE 

INTERVAL 

MAXIMUM 

L-M 

1/Tj^, log a 

.1^12 

.631 

.756 

0-S-D 

1/Tg, log a 

.1^8 

.661 

.824 

M-S 

Tp, log 0 

,118 

.527 

.657 

M-S 

Tp, log a/E 

.114 

.506 

.648 

M-H 

Tp, T^^log 0 

.116 

.517 

.660 

M-H 

Tp, T^, log o/E 

.110 

a\ 

CO 

• 

• .652 

RAB 

1/Tp, a 

.373 

1.159 

4.140 

GIVAR 

' 

.061 

.279 

.353 

GIVAR 

log Tp,-^ cr^/3^ 1/g 
^ °TU» 

.044 

.213 

.256 



TABLE 6— Summary of RMS comparisons. 


(jO 

U) 


-v:s 



Number 
of Data 

L-M 


0-S- 

-D 

M-S 


Alloy 

Points 

Terms 

RilS 

Terms 

RMS 

Terms 

RMS 

1 

eit 

7 

.159 

7 

.139 

7 

.220 

2 

75 

7 

.082 

7 

.086 

7 

.160 


26 

7 

• l6l 

7 

.149 

7 

.161 

5 

33 

7 

.063 

7 

.050 

7 

.089 

6 

26 

7 

.097 

7 

.057 

7 

.124 

7A 

52 

7 

,140 

7 

.131 

7 

.179 

8 

39 

7 

.178 

7 

.131 

7 

.233 

9 

38 

7 

.111 

7 

.094 

7 

.141 

llA 

iJ2 

7 

.134 

7 

.142 

7 

.122 

IIB 

41) 

7 

.132 

7 

.139 

7 

.122 

12 

24 

7 

.183 

7 

.191 

7 

.178 

13 

95 

7 

.288 

7 

.291 

7 

.291 

lij 

26 

5^ 

.064 

5a 

.074 

5^ 

.073 

15 

37 

7 

.088 

7 

.092 

7 

.100 

16 

33 

7 

.126 

7 

.132 

7 

.105 

17A 

103 

7 

’.202 

7 

.232 

7 

.201, 

17B 

105 

7 

.200 

7 

.228 

7 

.201 

18a 

100 

7 

.216 

7 

.251 

7 

.182 

18B 

104 

7 

.214 

7 

.252 

7 

.180 

19 

100 

7 

.265 

7 

.253 

7 

.308 

Average 



.155 


.156 


.169 


^Evidence of ill-conditioned solution for seven term model. 


M- 

■H 

RAB 

MCM 

GIVAR 

Terms 

RMS 

Terms 

RMS 

A 

RMS 

Terms 

RMS 

7 

.153 

6 

.209 

0 

.127 

9 

.106 

7 

.074 

6 

.245 

0 

.077 

11 

.055 

7 

.128 

4 

.247 

-.15 

.109 

• 7 

.088 

7 

.054 

3 

.304 

0 

.054 

9 

.043 

7 

.045 

4 

.298 

-.05 

.043 

6 

. 042 

7 

.121 

6 

.256 

0 

.131 

10 

.091 

7 

.137 

4 

.291 

-.05 

.115 

10 

.074 

7 

.111 

6 

.148 

.15 

.078 

8 

.068 

7 

.113 

4 

.179 

-.10 

.109 

7 

.100 

7 

.111 

. 4 

.218 

-.05 

.111 

9 , 

.099 

7 

.178 

6 

.385 

-.10 

.175 

6 ' 

.166 

7 ■ 

.282 

6 

.363 

-.05 

.290 

8 

.268 

5^ 

.059_ 

6 

.293 

0 

.056 

10 

.037 

.7 

.088’ 

6 

.321 

0 

.096 

7 

. 068 

7 

.103 

• 6 

.344 

-.15 

.072 

9 

.052 

7 ■ 

.200 

4 

.461 

0 

.198 

9 

.191 

7 

.198 

4 

.477 

0 

,201 

' 7 

.196 

7 

.182 

6 

.381 

0 

.186 

8 

.173 

7 

.180 

6 .. 

.410 

,0 

.187 

8 

.171 

7 

.276 

6 

.452 

0 

.350 

10 

.225 


.140 


;314 


.138 


.116 


ORIGINAL PAGE IS 
OF POOR QUALITY 



TABLE 7— T-Experimental and calculated stress-rupture data 
for Timken 35-15 stainless steel. 


Experimental 


Temp . , 

Stress, 

Time to 

ksi 

Rupture 

Hours 

1200 

21.0 

120 

1200 

19.0 

i7o 

1200 

18.0 

300 

1200 

13.0 

975 

1300 

16-.0 

60 

1300 

13.0 

’160 

1300 

11,0 

■ 300 

1300 

7.5 

■ 1300 

1400 

8.5 

120 

i4o6 

7.0 

400 

l400 

6,0 

900 

15P0 

6.0 

120- 

1500 

■4.9 ^ 

300- 

1500 

3.5 

950'“ 

1600 

6,0 

-2-0 - 

1600 

4.0 " 

170 

1600 . 

3.0 

■ 500 

1600 

2.5 ■ 

1300 

1800 

3.0: - 

22 

i860 

2.6 ' 

• 100 

1800 

1.5 

5‘66 

1800 

1.3 

1000 


Calculated time to rupture, hours 


L-M 

0-S.tD 

.RAB 

Exact 

Exact 

Exact 

81,36 

I4g.l6 

77.21 

i4o,47 

236.38 

110.82 

188.67 

303.14 

134.11 

1114.10 

1354.55 

401,32 

46.89 

. 52.16 

71.44 

136.51 

135.58 

150,14 

■322.53 .. 

292.38 

265.00 

2315.56 

1702.62 

910,22 

■166,00 

122.73 

228.84 

^ 427.33 - 

299,80 ' 

434.31 

-,905.32 

- 609.27 

7li.23 

138.89 

.96.32 

252.30 

'354.14. 

• 244.54 

490.24 

■ 1676.99 ■ 

;ii49.68 

1409.44 

■■ 25.56 

• 18.21 

73,96 

152.03 

■ 117.62 

302.83 

538. .66 • 

441,81 

765.48 

‘1200.88- 

. 1022,11 

1349.54 

. . 21,84 

24.58 

30.82 

110.91 

158.75 

139.88 

■ 351.35 

596.31 

’ 369.08 

623.51 

1151.78 

585,-77 



TABLE 8 — Summary of Larson-MLHer method correlattona 


Number of terms 

4 

RMS 

5 

7 

Alloy 

1 

r .ia <*2 

0. 16.00 

0 .1537 

2 

n.m7i 

n.oovo 

n,u«i7 

i| 

c.i/yc- 

0 ,lfe2f 

n.ibOo 

c 

0.03?.'^ 

0,073? 

0 ,Pb?5 

6 

7 

8 
Q 

0.100') 

0.0V77 

0.0970 

0,17t.? 

0.1527 

0.1404 

0.1 7 VI 

0.17V1 

(1.1777 

0.13A7 

0.12?>* 

0.1114 

llA 

IIB 

12 

13 

l^t 

15 

16 
17A 
17 B 
18A 
18B 

19 

(1.13 VS 

0.1350 

n.l337 

0,1397 

0.1340 

0.1317 

0.23&O 

o,iaf.i 

0.1833 

0.331 3 

0.3057 

0.26(11 

O.Cb77 

0,0643 

3. 

0.1>)97 

0.0991. 

0 . 067.. 

n.lb03 

0.1342 

0.1257 

r.?37i' 

0.2064 

'>.2020 

0.25&A 

n.204( 

0. 1994 

0.230V 

0.2245 

0,2161 

0.22A'’ 

0.22M 

0.213° 

0.V3H. 

('.35 3t 

0,2t53 


STD 

U 5 7 't 


0. 1"03 

0,1671 

0.1660 

0.7937 

0.1514 

0.0921 

0.0656 

0.6267 

0.195? 

0.1811 

0. 1652 

0 .6497 

0.0676 

n. 1.795 

0.0704 

0.3750 

0.1164 

n.l 067 

.0.1135 

n.5153 

0.1639 

0.1606 

0.1509 

0.7714 

n.i6°i 

0.1916 

0.1962 

n.no 06 

0.1466 

n.l 31'3 

0.1233 

0.6214 

0,1467 

0.3447 

0.1464 

0.6193 

0.1465 

0.1424 

0.1436 

0.6174 

0.2561 

0.2091 

0.2178 

1 .1320 

0.3365 

0.3137 

0.2993 

1,3935 

(..0736 

6. 0716 

a 

' 0,3201 

0.1565 

0.1072 

0. n.973 

0.6729 

0.1774 

0.1457 

0.1417 

0,7077 

0.2426 

0.2116 

0.209? 

0,9951 

0,2575 

0.2090 

n.2069 

1.0554 

0.2356 

0.2303 

0.2241 

n .9679 

0.2327 

0.2336 

0.2215 

0.9541 

0,4409 

C. 7,6 26 

0.2751 

1 .6096 


DPAVG DPMAX 


5 

7 

4 

5 

7 

(1.7020 

C.7159 

0.6455 

0.7960 

0.9100 

0.3649 

0.3632 

0.6633 

0.4218 

0.4410 

0.6029 

0.(654 

0,9486 

0.9214 

1.0204 

0.3443 

0.3131 

0,4068 

0,4077 

0.4036 

0.4618 

0.5223 

0.5700 

0.5662 

0,6494 

0.6797 

0.6497 

0.6091 

0.7370 

0.7530 

0.8219 

0.8597 

0.6631 

0.9337 

1.0632 

0.5631 

n.5410 

n .7412 

0,7294 

1.7116 

0,6176 

0.6366 

0.6664 

0.6766 

0.7051 

0.6062 

0.6241 

0.7014 

0.7037 

0.7229 

0,9359 

1,0142 

1.4289 

1.2125 

1.2660 

1 .2983 
0.3172 

1.25^4 

1.5469 

0.3576 

1.5613 

0.3772 

1. 1^586 

0.4606 

0.4282 

0.7275 

0.5297 

0.5359 

0,6309 

0.6303 

0.6909 

0.7562 

0.606^ 

0.6725 

0,6711 

1.6364 

0.9426 

1.0413 

0.8607 

0.6602 

1.1363 

1.0042 

1.1166 

0.9510 

0.9344 

1.0056 

1.0292 

1.1236 

0,9631 

0.9212 

1.0206 

1.1120 

1.1542 

1.4981 

1.1473 

1.6603 

1.5931 

1.2661 


^Evidence of ±-ll-conditioned 


solution. 
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TABLE 9~-Summary of 

Orr-Sherby-Dorn 

method correlations, 






RMS 



STD 



DPAVO 



DPMAX 

Number of terms 

4 

5 

7 

4 

5 

7 

4 

5 

7 

4 

5 

Alloy 

1 

0 . 1 692 

0.J417 

0.1367 

0.1954 

0.1476 

0.1470 

0.8149 

0.6202 

0.6264 

0.8619 

0.6899 

2 

0.201'S 

0.1072 

0,0860 

(I.20P1 

0.1109 

0.0903 

0,6643 

0.4636 

0.3823 

0,9103 

0.4950 

>t 

0.1819 

0.15111 

0.149S 

0.1978 

n ,1680 

0.1746 

6.8608 

0.7452 

0.8036 

0.9502 

0.8133 

5 

0.681? 

0.0611 

0,0495 

0.0872 

0.0663 

0,0557 

0.37P7 

0.2676 

0.2482 

0.3975 

0.3029 

6 

D.U661 

0.0574 

0.0566 

0.0718 

0.0639 

0.066? 

0.3128 

0.2634 

0.3045 

0.3262 

0.3110 

7 

O.lbli? 

0.1456 

0,1306 

Q.1647 

0.1531 

0.1404 

0.6907 

0.6480 

0.6048 

0.7164 

0.6680 

8 

0.136? 

0.1374 

0,1308 

0 , 14 66 

(1.147? 

0 .1 4-*4 

0.6P09 

0,6308 

0.6329 

0.6616 

0.6944 

9 

0.12o4 

0.1077 

0,0935 

0.1337 

0.1156 

0.1035 

0.5668 

0.4962 

0.4541 

0.5925 

0.5893 

llA 

0.1511 

0.1437 

0.1424 

■0.1588 

0.1531 

0,1560 

‘ 0.6707 

0,6539 

0,6803 

0,7124 

0,6901 

IIB 

0.1566 

0.141 3 

0,1393 

0.1644 

0.1501 

0.1519 

0.6930 

0.6393 

0,6605 

0.7603 

0.7149 

12 

0.2669 

0.1940 

0,1914 , 

0.2945 

0.?180 

0.2275 

1.293) 

0.9767 

1.0607 

1.4279 

1.1902 

13 

0,34?“ 

0.3193 

■ 0.?910 

0,3503 

0,3281 

0.3023 

1.44P5 

1.3978 

1.2643 

1.5942 

1.6033 

14 

0.0873 

f..0738 

a 

0.0949 

0.0821 

a 

8.4131 

0.3642 

a 

0.4261 

0.3991 

15 

0.1742 

0.09'46 

0.0925 

O.lflOH 

0,1073' 

0.1025 

0.7830 

0.4614 

0,4516 

0.8225 

0,4657 

16 

0.1797 

0.1416 

0,1317 

0,1917 

n. 1537 

0.1484 

' 0.6191 

0.6661 

0.6605 

0.8540 

0.7447 

17A 

0,2864 

0,?370 

0.23?? 

0,?9?1 

0,24?9 

0.P40P 

1.1962 

1.0016 

1.0013 

l.?249 

1.0250 

17B 

0.311)1 

P.?329 

0,2284 

o; 1P43 

0.P3P.7 

0.P364 

1.3290 

0.9830 

0,9030 

1.3996 

1.0309 

18A 

0.2656 

0.2650 

0.2505 

0,2713 

0.2719 

0.2597 

1 .1145 

1.1226 

1,0832 

1,1424 

1.1523 

16B 

0.270? 

0 . ?6M 

0,2S?0 

0.2755 

0.2727 

0.2609 

1.1297 

1.1238 

1.0855 

1,1678 

1.1764 

19 

0.510? 

0.4067 

0.2533 

0i5207 

0.4173 

0.?6?7 

2.1390 

1.7232 

1.0957 

2.1995 

1,7714 


Evidence of ill-conditioned solution. 


7 


o,7nc6 

D.4450 

D.9169 

0.3043 

o.3oae 

0.6624 

0.7569 

0.5902 

0.7378 

0.7542 

1.3414 

1.5164 

0,5349 

0.6237 

1,1031 

1,1771 

1,2127 

1.2777 

1.1595 





TABLE 10— Suaunary of Manson-Succop method correlations 




RMS 



STD 



DPAVO 



DPMAX 


Number of terras 

i| 

5 

7 


5 

7 

4 

5 

7 

4 

5 

7 

Alloy 

1 

0.2H17 

0.2256 

0.2204 

0.2496 

0.2350 

0.2335 

1.0412 

0.9673 

0.9951 

1.1062 

1.0916 

1.2462 

2 

C.204? 

0.1679 

0.1602 

0.2096 

0.1738 

0. 1682 

0.8714 

0.7264 

0,7119 

0.9425 

0.8035 

0.8258 

H 

0.1776 

0.1620 

0.1610 

0.1930 

0.1812 

0.1 883 

0.8404 

0,803.6 

0.8666 

0.8913 

0.8661 

0.9908 

5 

O.lObS 

0.0941 

o.oeflt 

0.1126 

0.1022 

0,0996 

0.4809 

0,4429 

0.4446 

0.4986 

0.4680 

0.5449 

6 

0.12t.6 

0,1239 

0.125b 

0.1376 

0.1578 

0.1445 

0.5990 

0.6111 

0.6647 

0.6761 

0.6824' 

o.eoso 

7 

0.2023 

0.1030 

0.176ft 

0.2105 

0.1925 

0.1922 

0.8830 

0.8145 

0.8279 

0.9291 

0.8500 

0.9079 

8 

0.231'» 

0.2339 

n.2330 

0.2469 

0.2505 

0,2572 

1,0455 

1.0735 

1.1276 

1.1286 

1.1731 

1.3518 

9 

0.1639 

0.1516 

0.1409 

0.1733 

0,1627 

0.1559 

0.7346 

0,6981 

0.6644 

0.7706 

0.7975 

0,8695 

llA 

0.1300 

0.1241 

0.1219 

0.1369 

0.1322 

0.1336 

0.5779 

0.5646 

0.5826 

0,6197 

0,5990 

0.6316 

IIB 

0.1317 

0.1249 

0.1219 

0.1301 

0.1327 

0. ) 329 

0.5822 

0.5652 

0.5778 

0.6514 

0.6310 

0.6609 

12 

0.23S'- 

0.1U21 

0.1779 

0.2579 

0.2046 

0.21J14 

1.1325 

0.91.67 

0,9856 

1.2058 

1.1163 

1.2468 

13 

0,3492 

0.3050 

n.2914 

0.5560 • 

0.3146 

0.3027 

1 .4666 

1.2996 

1.2661 

1.6739 

1.5537 

1.5409 

li| 

u.oaoi- 

0.0720 

a 

0.nft75 

o.oftin 

a 

0.3810 

0.3591 

a 

0.3926 

0.3936 

a 

15 

0.1673 

0.1153 

0.0994 

n.1771 

0.1240 

0,1110 

0.7520 

0.5331 

0.4889 

0.7838 

0.5629 

0.5796 

16 

U.15!<9 

0.1157 

n.lobl 

0.1663 

0.1256 

0. )1R4 

6.7106 

0.5443 

0.5268 

0.7567 

0.6055 

O'. 6570 

17A 

['.26,5^ 

0.2104 

0.2012 

0.2736 

0.2239 

0.2085 

1.1222 

0.9231 

0.6679 

1.1497 

0.9447 

0.9559 

17B 

n.279b 

C.2205 

0.2007 

0.2849 

0.2260 

0,2077 

1.1677 

0.9307- 

0.8637 

1.2206 

0.9669 

1,0323 

18a 

0.194-J 

0,1942 

0.1022 

0.1969 

0.1992 

0.188° 

0.8171 

0.8227 

' 0.7876 

0.8376 

0.8482 

0.8827 

18 b 

0.1907 

0.1954 

0.1003 

0.2027 

0.2003 

0.1667 

0.8509 

0.82S3 

0.7768 

0,8621 

0,6620 

0.9142 

19 

0.4161 

0.364(i 

0.3076 

0,4247 

0.3735 

0.3192 

1 .7446 

1.5423 

1.3313 

1,8299 

1.6135 

1.4131 


^Evidence of ill-conditioned solution. 



TABLE ll“Summary of Manson-Haferd method correlations 


RMS 

Number of terms ^57 


Alloy 


1 

0 .212? 

0.1653 

0,1533 

2 

, 0.1893 

0.1145 

0,0741 

A 

0.1722 

0,1312 

0.1264 

5 

0,0930 

0.0597 

0.0541 

6 

0.0481 

0.0455 

0.044B 

7 

0.1475 

0.1294 

0 .121? 

8 

0.1524 

0,1390 

0.1365 

9 

0.1476 

0. 1299 

n.iiii 

llA 

0.1248 

0.1143 

0.1120 

IIB 

0,1287 

0.1144 

0.1111 

12 

0.217fa 

0,1821 

0.1779 

13 

0.3446 

0,3041 

0.2819 


0.0773 

0,0507 

a 

15 

0.1675 

0.1015 

0.0863 

16 

0.1557 

0.1148 

0.103? 

17A 

0.2677 

0.2170 

0.1996 

17B- 

0.2758 

0.2179 

0.1984 

18a 

0.1946 

0.1936 

0.1622 

16B 

0.1967 

0.1948 

0.1604 

19 

0.4129 

0.3619 

0.2760 



STD 



4 

5 

7 

4 

0,2192 

0.172? 

0.1624 

0.9144 

0.1946 

0.1189 

0.0779 

0.6060 

0.1872 

0.1460 

0.1502 

6.6148 

0,0992 

0.0649 

O.C610 

0.4?39 

0.05P3 

0,0507 

0.0525 

0.2277 

0.1535 

0.1361 

0.1303 

0.6437 

0.1609 

0.1409 

0.1507 

0,6813 

0.1560 

P.1594 

0.1231 

0.661 1 

0.1312 

0,1217 

0.1 ?35 

0.5540 

0.1380 

0.1215 

0.1212 

0.5689 

0.2384 

0.2046 

0.2114 

1 .0495 

0.3523 

0.3124 

0 .? 9?9 

1.4505 

0,0840 

0.0654 

a 

0.3695 

0.1771 

n.iogi 

0.0980 

0.7520 

0.1660 

0.1247 

0.1 16? 

0.7093 

0.2730 

n.?224 

0.?069 

1 ,1199 

0.2612 

n.?232 

0.2054 

1 ,15?5 

0,1986 

0,1986 

0.1889 

0.6158 

0.2006 

0,1997 

0.1866 

6 ,82?4 

0.4214 

0.3713 

0, ?86? 

1 .7310 


DPAVG DPMAX 


5 

7 

4 

5 

7 

0.7235 

0.6921 

0.9716 

0.7998 

0.8466 

0.4954 

0.3295 

0.8668 

0,5535 

0.385? 

0.6472 

0.6906 

0.9018 

«0.7294 

0.607? 

0.2809 

0.2712 

0.4677 

0.3406 

0,3510 

0.2247 

0.2413 

0.2468 

0.2594 

0.2988 

0.5759 

0.560B 

0.6907 

0.6429 

0.6534 

0.6360 

0.6605 

n.7561 

0.7468 

(1,8061 

0.5975 

0,5397 

0,8048 

0.7819 

0.7105 

0.5196 

0.5388 

0.5986 

0.5721 

0.5898 

0.5172 

0.5266 

0.6492 

0.6060 

0.6095 

0.9157 

0.9847 

1 .3147 

1.1877 

1.2502 

1.2929 

1.2244 

1 .6561 

1.5899 

1.6379 

0.2896 

a 

0.4131 

0.3461 

a 

0.4690 

0.4315 

0.0203 

0.5505 

0.5441 

0.5398 

0.5172 

0.7821 

0.6419 

0.6680 

0.9170 

0.6614 

1.1706 

0.9968 

1 .0388 

0,9193 

0.8537 

1.2466 

1.0843 

1.1137 

0,6200 

0.7678 

0.0540 

0.8953 

0.9548 

0.8226 

0.7769 

0.8926 

0.9539 

0.9775 

1,5331 

1 .1934 

1 .8026 

1.6363 

1.3242 


^'Evidence of ill-conditioned solution. 



TABLE 12 — Summory of Rabotnov method correlations, 


a 




HMS 



STD 



DPAVG 


DPMAX' 


5f terms 

3 

it 

6 

.3 

4 

■ 6 

3 

4 

6 

3 

4 

6 

Alloy* 

1 

0,4387 

0.2185 

0.2091 

o '. 4494 

0.2257 

0.2178 

0.0000 

0,0000 

0 . 9l'58 

o.nooD 

0.0000 

1.0270 

2 

0.3800 

0.2458 

0.2452 

0.3879 

0.2526 

0.2538 

1.8936 

1.0257 

1.0269 

3.5595 

1.6010 

1.6469 

0.3033 

0.2473 

0.2644 

0.3225 

0 .' 268 a 

0.2942 

1 .4667 

' 1.4262 

1.4309 

3.4359 

3.8611 

5.3300 

5 

0.3040 

0.3053 

0.3132 ■ 

0 . 31'86 

0.3257 

0.3400 

1.7321 

1.6602 

1.7326 

3.1460 

2.9195 

3.2385 

6 

0.3228 

0 ,' 29‘75 

0.3425 • 

0.3432 

0.3235 

0 . 381'1 

1 .6292 

1,6693 

1.6632 

4.3534 

4.2000 

4.4432 

7 

0.4090 

0,2674 

0 . 2561 . 

0;4214 

0,2783 

0 . 2694 , 

1 .9803 

1.8364 

1.3095 

6.4589 

4.1598 

3.6349 

8 

0.5885 

0.2912 

0.3257 ' 

0.6105 

0.3074 

0 . 348 '' 

2.9696 

1.6791 

1.9132 

12,1315 

3.4606 

S.6765 

9 

0.2364 

0.1523 

0.1476 

0;2464 

0.1610 

0.1611 

1,0861 

0.7106 

0.7195 

1.3500 

1,0488 

1.1497 

llA 

0.5037 

0 . 1-791 

0.1915 

0.5227 

0,1882 

0 . 2'040 

■ 1 . 476 ? 

0.6655 

0.9193 

3.7275 

2.4516 

2.8222 

IIB ■ 

0.,4127 

0,2180 

0 ^ 2326 " 

0,4276 - 

0.2287 

0-.2471 

1 . 874’8 

1.0452 

1.0775 

4.8993 

3.7546 

4.5071 

12 

0.4197 

0-,3891 

0.’3654 ^ 

0.4467 

o '. 4262 

■ 0,4332 

* 2.4168 

1.8572 

1.9490 

6.8192 

3.1121 

3.1838 

13 

0 . 814O 

0,3982 

0 . 36'26 

O.a'272 

0.4068 

0.3746 

3.7109 

2.0841 

1.9390 

7.4721 

5.0910 

3.7243 

l^i 

0,3240 

0.2945 

’ 0;2927 

0.3445 

0,3202 

■ 0.3257 

- 0.0000 

1.4057 

1.4524 

0,0000 

1.5688 

1.5985 

15 

0;3249 

0,3306 

0..3210 

0.3369 

0.3501 

0.3451 

1 . 626 ? 

1.6262 

J .4666 

3.4998 

5.7345 

4 . 536U 

16 

0.5893 

0,3604 

0.3436 . 

0.6181 

0,3845 

0.3732 

1.2650 

1.1264 

1.'1590 

6>260 

3.8236 

4.1400 

17A - ■ 

- 0.4785 

0.4608 

0,4638 

0.,4857 

0,.4700 

0.4755 

2,6635 

2.2600 

2.3512 

13.8266 

7.2185 

8.4368 


0.4892 

0.4770 

0.4770 

. 0.4963 

0.4863 

0.4887 

?,.6286 

2.4270 

2.4692 

13.1340 

10.3149 

10.8354 

18a 

0.3896 

0.3921 

0.3807 

0.3955 

0.4002 

0,3906 

1 .7884 

1.8021 

1.8257 

5.8525 

5,8298 

11.6130 

18 b 

0.4166 

0,4116 

0.4101 

0.4229 

0.4198 

0,4203 

1 . 94 ?? 

1.7542 

1 .8609 

16.5520 

5 , 757.3 

14.9883 

19 

0.6105 

0.4668 

0.4516 

0.6196 

0.4764 

0,4634 

3.4160 

2,0414 

2 , 0226 . 

10.1529 

3.4808 

3.3542 


a 


Based upon backtransformed log time values. 



TABLE 13 — Summary of GIVAR method correlations. 


Alloy . 
1 

No, of 
Observ, 

Prime Variable 
Transformation 

No, of 
Start 

Variables 

"Best” 

RMS 

STD 

• DPAVG 

DPMAX 

, 64 

Temp, 

1/T 

Stress 

12 

9 

.1060 

.1130 

.4921 

.6387 

2 

• 75 

log T 

„ l /3 

‘12 

11 

.0551 

,0592 

.2548 

.3099 

4 

26 

1/T 


12 

7 

.0883 

.1033 

.4754 

.5273, 

5 

33 

log T 

0 

11 

9 

.0426 

.0499 

,2284 

, ,2799 

6 

26, 

1/T 

log a 

9 

6 

.0418 

,‘0476 

,2150 

.2485 

7 

52 

log T 


12 

10 

.0910 

.1013 

.4471 

,5155 

8 

39 

1/T ■ 

0 

11 

10 

.0744 

.0863 

.3910- 

.4971- 

9 

38 

log T 

0 

11 

"8 

.'0677 

,0762 

.3384 

.4169. 

llA 

42 

1/T 


12 

7 

.0997 

.1092 

.4764 

.5360 

IIB 

44 

1/T 

a 

11 

9 

.0985 

,1104 

.4894' 

.5686 

12 

24 

1/T 

cr 

11 

6 

.1664 

,1-921 

.8784 

1.0377 

13 

95 

1/T 

log a 

11 

8 

.2677 

♦ 2797 

1.’1756 

' 1.4435 

14 ■ 

26 

log T 

1^1/3 

12 

10 

.0368- 

.0456 

.2213. 

' ,2547 

i5 

37 

iog T 

a 

li 

7 

.0683 

»0758 

. .3337 

.4185- 

16 

33 

1/T 

a 

11 

■ 9 

.0520 

■ .0610 

.2788 

.3529 

17A 

103 

1/T 

a 

11 

9 

.1913 

.2002 

^'8417 

,8450 

.9649' 

17B 

105 

1/T 


12 

7 

.1963 

.2032 

9-056 

18A 

100 

1/T 

0 

, 1/3 

11 

8 

.1726 

.1799 

.7541 

.’8583 

18B 

104 

' 1/T 

12 

8 

.1712 

.1782 

.7451 

.8174 

19 

100 

1/T 

jjl/3 

12 

10 

.2248 

.2369 

1.0d28 

1.0978 
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RESIDUAL = PREDICTED LOG TIME "ACTUAL LOG TIME 



PREDICTED LOG TIME 




Fig. 2 Comparison of regression residuals for Rabotnov simulated data set. 








Fig, 4— Regression residuals for alloy 4, GIVAR analysis. 




Fig. 5— Cumulative normal distribution of regression residuals for alloy 4 
GIVAR analysis. 
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COMPUTER PROGRAM FOR PARAMETRIC 
ANALYSIS OF CREEP-RUPTURE DATA^ 

By- 

Donald R. Rummler 
ABSTRACT 

A computer program which uses several parametric model 
equations to analyze creep-rupture data is presented in detail. 
The model equations include the Larson-Miller, Orr-Sherby-Dorn, 
Manson-Succop , Manson-Haferd, and Rabotnov parameter methods. 
Standard multiple regression techniques are used to analyze data 
with respect to each model equation. In addition to the usual 
regression statistics, the program calculates statistical 
intervals including confidence and prediction intervals. 

Graphical output Includes a residual plot with respect to the 
dependent variable and a cumulative distribution of the 
residuals. The computer input and output, in printed and plott.ed 
form, for sample problems are presented to aid the user in 
setting up and running the program. 



SUMMARY 


A computer program which uses several parametric model 
equations to analyze creep-rupture data is presented in detail. 
The model equations include the Larson-Miller, Orr-Sher by-Dorn, 
Manson-Succop, Manson-Haf erd, and Rabotnov parameter methods. 
Standard multiple regression techniques are used to analyze data 
with respect to each model equation. In addition to the usual 
regression statistics, the program calculates statistical 
intervals including confidence and prediction intervals. 
Graphical output Includes a residual plot with respect to the 
dependent variable and a cumulative distribution of the 
residuals. The program, its subroutines and their variables are 
listed and defined. The computer input and output, in printed 
and plotted form, for sample problems are presented to aid the 
user in setting-up and running the program. The development of 
the parameter model equations and the use of statistical 
Intervals is discussed. ' • ' 



INTRODUCTION 


The Importance of creep-rupture data analysis has led to a 
large number of papers which either propose new parametric 
analysis approaches (refs. 1, 2, 3, and k, for example) or offer 
detailed comparisons of different parametric methods (refs, 4, 5> 

t 

and 6). Most parametric methods for creep-rupture data analysis 
are empirical. Consequently, it is common practice for the data 
analyst to fit the creep-rupture data at hand to a variety of 
parametric model equations to select the most appropriate 
analysis method. 

Although several analysis methods have been presented in 
general terms (ref, 6, for example), there is no widely used, 
efficient computer program tailored specifically to the 
parametric analysis of creep-rupture data. In addition, most 
methods do not include generation of statistical intervals to aid 
in the selection of the "best” parametric model equation for a 
particular set of data. 

This paper describes the development and use of a computer 
program for the parametric analysis of creep-rupture data. The 
program includes provisions for the analysis of five different 
parameter methods. The parametric equations used and the 
statistical quantities calculated are discussed. The computer 
program input and output, in printed and plotted form, for three 
sample problems are presented to aid the user in setting up and 
running a problem with the program. 



PROGRAM description 

The computer program (PARAM) was developed to analyze and 
correlate creep-rupture data utilizing a variety of parametric 
method model equations. For each model equation, a function of 
the time to a particular creep event (such as time to 0.005 
strain) is the dependent variable. Functions of stress and 
temperature are the only correlating independent variables. The 
major features of the program are as follows: 

(1) The method of least squares is used to establish the 
coefficients for the parametric model equation selected for 
analysis. 

(2) Provisions are made for analysis with four widely used 
time-temperature methods (Larson-Miller, Orr-Sherby-Dorn, 
Manson-Suecop, and' Manson-Haferd) and one time-stress (Rabotnov) 
method, 

(3) Polynomial forms of the parametric model equations up 
to the fifth order are included, 

(4) Multiple analyses can be accomplished during a single 
computer run. 

(5) In addition to the usual regression statistics, the 
program calculates the maximum and minimum value of each 
independent variable, as well as' its range and average value, 

(6) The program also calculates the relative influence, 
contribution to the sums of squares, and warns of coefficient 
solution errors for each Independent variable, 

(7) Listings are made of the observed and fitted values of 



the dependent variable in both regression and real variable 
coordinates, 

(8) Two statistical intervals, the 95 percent confidence 
and the 95 percent prediction, are approximated and calculated 
for each observation. 

(9) Residual plots are made to indicate how the regression 
residuals are distributed over all of the fitted values of the 
dependent variable and whether they are normally distributed, 

PARAM was written in FORTRAN IV language for the Control 
Data 6000 series digital computer under the SCOPE 3.0 operating 
system. The program is dimensioned for a maximum of 5 input 
variables, a maximum of 10 derived Independent variables and a 
maximum of 200 observations for each data set. It requires 
approximately 60,000 octal locations of core storage, A source 
listing of the main program and its subroutines is presented in 
appendix A, A detailed description of the matrix equation 
solution subroutine MATINV and the plotting subroutines PSEUDO, 
DDIPLT and CALPLT are presented in appendix B. 

ANALYSIS 

The analysis utilizes standard least squares multiple 
regression analysis techniques (refs, 7 and 8) to solve 
parametric equations of the following form: 

Y = bo + b^X^ + b^X2 + + b^^X^ (1) 

where Y = fitted value of dependent variable 

Xi,Xg . . . ,, Xj^ = independent variables 



= estima-ted Y intercept when all = 0 

^1^ ^2^ — — — = estimated coefficients of Independent 

variables 

Specifically, the equation forms chosen for each of the 
parametric methods selected are as follows: 

Larson-Miller CL-M) 

Y = log t = b^ + b^/Tp + bgS/Tj, + + b^^ sVTj^ + 

bgsVTj, + .bgsS/T^ , C2) 

Orr^Sherby-Dorn (0-S--D) 

r = log t = b^ + b^/Tj^ + b 2 S + b, + ’b.S^ + b.S^ + 

^6^^ (3) 

Manson-Succop (M-S) 

Y = log t = b^ + b^Tj. + b 2 S + b^S^ + + b^S^ + 

CO 

Manson-Haferd (M-'H) 

Y = log t = b^ + b^T^ + b 2 T^S + b^T^S^ + + 

b^ToS^ + bjT^s'' + bgT„s5 (5). 

Rabotnov (RAB,)' 

Y = = b^ + bj^/oTp + bj/oT^ + bj/oll + b^/tjT^' + 

bj/oll + bg/oT® (6) 

where 

t = time to a particular creep-rupture event, rupture, 
for example ' • , 

S = log a 

a = applied stress 

f3 - 
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Tp = temperature, °F 

= temperature, Kelvin 
= temperature, Ranklne 
Tq = offset temperature = 

^1, ^ “ constants -estimated by method of least 

squares , 

Both the M-H and RAB techniques require the use of iterative-, 
non-linear multiple regression techniques to estimate all of the 
constants. 

Each parametric equation can be analyzed in' truncated form 
since the number of equation terms (LLO) is selected with input 
case control cards. 

The development .of each of the parametric method model 
equations is presented in appendix C. 

PROGRAM USAGE 

To submit a problem, information is normally entered on 
punched cards. Pour types of information cards (option, case 
control, data set identification, and. data) are the only input 
required. Output includes listings and plots. 

Input 

The option card controls both the printed and graphic output 
of the program, ' It also establishes the Initial values to be 
used for the iteratively modified .constants for the Manson-Haferd 
and Rabotnov parametric analyses. The case control cards 
determine the paraimetrie equation forms to be evaluated and their 



degree of truncation. A data identification card and the data 
cards complete the deck set up. The input card order, format, 
permitted values and comments follow: 


Option card 

Column FORTRAN Variable 
5 INPUT 


10 OUTPUT 


11 to 20 TA 


21 to 30 RA 


(215, 2F10.0) 

Value Comments 

0 No listing of input cards 

1 List data set I.D., option, 
and case control cards 

2 List 1 + data observa- 
tions 

3 List 2 + regression varia- 
bles for first five 
observations 

0 No listing of residuals 

1 List regression residuals. 

2 1 + list back transformed 
residuals 

3 2 + regression residual 
plots 

Initial value for constant 
in non-linear M-H equation; 
A value of -5000.0 is 
recommended 

Initial value for constant 
in non-linear RAB equation; 


A value of 0.2 is recommended, 



Column 


FORTRAN Variable Value 
NPAM(I) 


LLO(I) 


Comments 

Parametric expression to 
be evaluated 
Larson-Miller- 
Orr-Sherby-Dorn 

% 

Manson-Succop 

Manson-Haferd 

Rabotnov 


2to6 Number of coefficients to 

be determined for parametric 
expression selected, see 
ANALYSIS section of 


paper, 


The program is dimensioned. for a maximum of 20 case control cards. 
During a single computer run, a data set can be evaluated with 20 
different parametric model equation forms. A blank card must 
follow the last case control card. . 


Blank Card 

Data identification card 
Column FORTRAN Variable 

1 to 80 TYPE 


(8a10) 


Comment 


Data I.D. Any characters in 
‘columns 1 to 80'. This title 
is included in all listed output 
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Data . cards 


C3F12.0) 


Column 

FORTRAN Variables 

Comments 

1-12 

RS(I,1) 


Time to a particular creep event 

13-24 

RS(I,3) 


Temperature, °P 

25-36 

RSCI,2) 


Applied stress 


The program is dimensioned for a maximum of 200 observations in a 
data set. Round-off errors can be minimized by limiting the range 
of the variables. This range reduction is helpful since most 
creep-rupture data is ill-conditioned Csee refs, 7 and 8), 

Last data card must be followed by a blank card. 

Blank card 

More than one set of data may be analyzed with a single set up of 
the option and case control cards. To analyze additional data set 
during a single computer run, assemble the deck as follows: 

Option card 
Case control cards 

Blank card First data set 

Data identification card 
Data cards 
Blank card 

Data identification card 

Data cards Second data set 


Blank card 



Data identification card 


Data cards Third data set 

E-lank card 

As many data sets as desired may be analyzed during a single 
computer run with this type of deck setup. 

Output 

Examples of printed and plotted output are presented in 
the discussion of sample problems. Most of the output 
headings are self-explanatory or standard statistical terms 
(refs. 7 and 8). Some headings are abbreviations of 
standard terms and/or require additional description. These 
headings and brief descriptions ^ in the order of their 
appearance for the printed output are as follows; 

Heading Description 

STANDARD ERROR Standard error of estimate is square 

root of residual mean square, sometimes 
called residual root mean square 

MULT, CORREL. 

COEPP. SQUARED The multiple correlation coefficient 

squared, sometimes called coefficient 
of determination 
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MIN 

MAX 

Y 

X1-X(L2) 
VARIABLE 
COEP. P.I. 
S.E. COEP. 
T 

RAN X(I) 
RINP(I) 

PSUM 


The minimum value of indicated variable; 
independent variables are in tabular form 
The maximum value of Indicated variable 
Tabulated values of Independent variable 
Tabulated values of independent variables; 
L2 is number of variables in case 
Transformation required for parametric 
method being evaluated 
Calculated' coefficients for the fitted 
equation, indexed by I starting with b^ 
Estimated standard error of the coeffi- 
cient . ' - • 

t 

COEF.P(I)/S.E. COEP. ' ' ' 

Range of independent variable 

Relative influence of independent * 

variable, (COEP.P(I) (RANX(I) ) 
y RANGE 

The fraction of the total sums 
of squares explained by an inde- 
pendent variable; corrected for 
those Independent variables which 
proceed it in the listing 
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CERR 


95 PERCENT 
PREDICTION 
INTERVAL 
STATISTICS 


REAL TIME FACTOR 
RESIDUALS - 
REGRESSION 
SPACE 
RESIDUAL 

PCTERR 


The percentage difference 
between MATINV and Gaussian 
elimination solutions for 
coefficient; values in excess 
of 0.01 suggest round-off 
errors due to ill-conditioned 
normal equations 

The 95 percent prediction interval for 
a single future observation is estimated 
for each observation in regression 
variable space; these values are back 
transformed into log time space to calcu- 
late average and maximum values; values 
for the t distribution are approximated 
with a third order polynomial in log 

(degrees of freedom) 

(WIDTH) 

Values listed under this heading are 
in terms of the regression dependent 
variable coordinates 

Observed value of dependent variable-cal- 
culated value of dependent variable 
(100) (RESIDUAL) 

y 
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ORDER 


The rank order of the residual in regression 
coordinates; the rank order of the 
PCTERR in real space coordinates; 
ordered with respect to the largest 
absolute value, 

CIMIN Estimated lower limit of 95^ confidence 

Interval for the mean 

CIMAX Estimated upper limit of the 95/5 

confidence interval for the mean 
PIMIN Estimated lower limit of 95^ prediction 

interval for a single future observation 
PIMAX Estimated upper limit of the 95^ 

prediction interval for a single future 
observation 

The values of the t distribution required for the 

calculation of the statistical intervals are approximated 

♦ 

with the following expression: 

TVALUE = (10.0)^^ - . 

where 

T1 = 0.86186 - 0.98427 DP + 0.58495(DF)^ 

- 0.11594(DP)^ 

DP = residual degrees of freedom for regression. 

The graphical output of the program includes a 
plot of the residuals with respect to the calculated value 
of the dependent variable (FITTED Y) and a cumulative 
normal distribution of the residuals (ZP NORMAL), Por the 


1’4 
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ZP NORMAL plot^ the plotting points for the abscissa, 

P, are In terms of the inverse of the standardized normal 
distribution and are calculated in the following manner:, 
for PZ = 0 0.5 

ZP^ = I. 0 M 5 I + 4 . 3598 XP + 3.^6o6,(XP)^’ + 1.9088(XP-)^ 

+ 0.5446(XP)^ + 0,06o8(XP)5 

where XP = log FZ 

PZ = 'U - 3/8>/CN + 1/4) 

j =1, 2, ---N when the residuals are arranged in 
order of Increasing magnitude, 
for PZ = 0.5 1.0 

XP = log (1-FZ) 

ZP2 = -ZP^ 

The ZP expression approximates the inverse of the standard 
normal distribution. 
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SAMPLE CASES 


Three sample cases are presented to illustrate operation 
of the computer program and a method for rapidly selecting the 

most applicable parametric equation for a single set of 

♦ 

creep-rupture data. The data are for a type 316 stainless steel 
(ref. 5). The three sample cases described in this section 
required a total of 10,9 seconds of CDC 6600 CPU time to compile 
and run. 

Case 1 

For this case, all five .parametric method.s in second 
degree form were used to correlate the data. The purpose of 
this case was to quickly scan the parametric models to select a 
single parameter for further study. Output was minimized by 
using INPUT = 1 and lOUT = 0. The program input and output for 
case 1 are presented in Figures 1 and 2, respectively. 

When compared to the other four parameter methods, the 
0-S-D method had the highest MULT. CORREL. COEF. SQUARED, the 
lowest AVERAGE and MAXIMUM WIDTH of the 95% prediction Interval. 
It also had the lowest STANDARD ERROR of the four 
time-temperature parameters. 

Case 2 

Based upon the results of case 1, the Orr-Sherby-Dorn 
parameter (NPAM = 2) was selected for further evaluation. The 
purpose of this case was to quickly determine the degree of the 
O-S-D expression which would provide the best correlation of the 
data. Once again, output was minimized (INPUT = 0, lOUT = 0). 
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The program input and output for case 2 are presented in Figures 
3 and 4, respectively. 

With respect to MULT. CORREL. COEP. SQUARED, there is no 
appreciable improvement in the correlation produced by 
increasing the degree of the polynomial expression. However, 
the STANDARD ERROR shows a steady decrease as additional 
variables are added up to the fifth order expression where it 
increases slightly. The T values for this fifth order ex- 
pression clearly illustrate the inflation of the standard error 
of the coefficients which this high level of co-linearity 
produces. The CERR value for I = 2 (XCl) = LOG STRESS) suggests 
that the solution matrix was ill-conditioned because the two 
methods of solution do not agree. 

The RESIDUAL SUMS OF SQUARES for the fourth order 
expression is approximately 30 percent lower than the third 
order expression. Although significant differences between the 
other correlation indications are not apparent, the fourth order 
expression is selected for further evaluation. 

Case 3 

Pinal verification of the fourth order expression selected 
in case 2 requires the full output capabilities of the program 
(INPUT = 3» lOUT = 3). The input and output for this case are 
presented in figures 5 and 6, The output includes a listing of 
the first 5 values of the regression variables, residuals and 
statistical Intervals in regression and back transformed 
coordinates and plots of residuals with respect to the 



calculated dependent variable (Y FITTED) and with respect to the 
normal cumulative distribution. The most important part of the 
verification of the fourth order expression is the examination 
of the residual plots. These plots suggest that the residuals 
have a zero mean and are randomly distributed with respect to 
the PITTED Y and that their cumulative distribution is normal. 
These two characteristics of the residuals are necessary for the 
calculation of valid statistical Intervals, 

The method selected for determining the "best” parametric 
equation for a set of data was used primarily to demonstrate the 
capabilities of the computer program PARAM. For other methods 
see references 4, 5, and 6. For a further discussion of the use 
of statistical intervals, the reader is referred to references 7 
and 11. 

CONCLUDING REMARKS 

A computer program specifically developed for the 
parametric analysis of creep-rupture data has been discussed. 

The equations used for the analysis of five parametric methods 
and the computer program used to implement the analysis are 
given. 

The computer program is versatile, allows rapid assessment 
of parametric methods for creep-rupture data, and has a 
relatively small core storage requirement. In addition to the 
statistics which are usually calculated and output by multiple 
regression programs, the program outputs the 95^ confidence 
interval on the mean and the 95^ prediction interval for a 



future observation. Residual plots are provided to .assess the 
validity- of the calculated statistical intervals. 
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APPENDIX A 

SOURCE LISTING OF PROGRAM' PARAM 


C 

C 

C 

C 

C 

c 
c 
c 
c 
c 
c 
c 
c 
■ c 

ro 
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PROGRAM PARAM ( I NPUT » OUTPUT » PUNCH « TAPE5= I NPUT . T APE6=OUTPUT ♦ 0000000 1 

1 TAPE7=PUNCH) ' OOOOOOOH 

PARAM 00000003 

PROGRAM FOR PARAMETRIC ANALYSIS OF CREEP-RUPTURE DATA 00000004 

COEFFICIENTS FOR PARAMETRIC MODEL EQUATIONS ARE DETERMINED BY 00000005 

METHOD OF LEAST SQUARES 00000006 

Y= B0+B1X1+B2X2 00000007 

parametric METHODS INCLUDE 00000008 

LARSON-MILLER (L-M) 00000009 

ORR-SHERBY-DORN(O-S-D) 00000010 


(LI ) 


MANSON-SUCCOP { M-S ) 

MANSON-HAFERD (M-H ) 

RABOTNOV (RAB ) 

DONALD R. RUMMLER , ' 

NASA-LANGLEY RESEARCH CENTER. HAMPTON, VA. , 1976 

ARRAYS WHICH DEPEND ON NUMBER OF ' OBSERV AT I ONS IN DATA SET 
DIMENSION AA{ 200), -GY( 2'00 ) • ' C 1 m’AX ( 200), CIMIN( 200) 
DIMENSION ERRPER( 200), f \ 200, 10), I PERM ( 200), PYMAX( '200) 
DIMENSION PYM!N.{ ‘ 200 ) , RIS( 200), 'RS ( 200,5), TEMPI 200), Y< 
DIMENSION 2P(200)-. 

ARRAYS WHICH DEPEND ON NUMBER’ OF VARIABLES IN REGRESSION MODEL 
NUMBER OF INDEPENDENT VARIABLES CL2 > 

CERR'd'O),- PARCiO) 

PARI < Ip), SB(lO)i SSR(IO), SUMAIIO), 

SUMP2 ( 1 p ,4 0 )., SUMX(IO)', S'UMXY ( I O ♦ 1 0 ) 

SUMX2{l-b;), T('IO), XMAXdO), XMIN(IO) 

XMEAN{ 1.0 ) ■ 

COEFFICI-ENTS DETERMINED (L3) 

D( 1 1 , 1 1 ) , ‘dD( 1 1 , 1 1 ) , Enifl), G(ll 


DIMENSION 
D I MENS I ON 
DIMENSION , 
DIMENSION-. 
DIMENSION’ 
number OF 
D I MENS I ON 
DIMENSION 
NUMBER OF 
DIMENSION 


SUMB ( 10 ) 

. SUMXinO, 
, XRAN(IO) 


10 ) 


1 ■) , INdEX( 11,2) 


IPIVOTJ I 1 ) , X( 1 I ) 
COEFFICIENTS +1 <N3) 

A ( 12, 12 ) , B( 12, 12 > 


ARRAYS which DEPEND UPON OTHER FACTORS 


000000 1 1 
00000012 
00000013 


00000016 

00000017 

ooooooia 

200)00000019 
00000020 
00000021 
00000022 
00000023 
00000024 
00000026 
00000026 
• 00000027 

00000028 
, 00000029 
00000030 
00000031 
00000032 
00000033 
+ 
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IV) 
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C 

C 
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c 
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c 

c 
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c 

c 
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NUMBER OF CASES 

DIMENSION LLO < 20 ) * NPAM(20) 

MI SC 

DIMENSION TYPE<8)» IN(2)» VAR(30)« PAM (5) 

DATA (PAM ( I ) « I =1 »5 )/ 3HL-M* 5H0-S-Di 3HM-S« 3HM-H» 3HRAB/ 

DATA (VAR ( I ) » 1 = 1 » 30 ) /3H 1 /T * 3HS/T « 6HS**2/T ♦ 6HS**3/T < 6HS**4/T » 
l6rtS**5/T ♦ 

23H 1/t » 1 HS » 4HS»*2 » 4HS**3 4 4HS**4 » 4HS**5 1 
31HT, 1HS»‘4HS**2»4HS**3« 4HS**4 ,4HS**5» 

42HDT^4HDT*S«7,HDT*S-»*2»7HDT*S**3i 7HDT*S**4 , 7HDT*S**5 , 

55Hl/L*T,aHl/L*T*-!f2.3Hl/L*TX-*3.SHl /L*T**4 ,8H1 /L^^T**5,SHl/L*T**6/ 

LI = NUMBER OF OBSERVATIONS IN DATA SET ^ - 

LI IS determined by PROGRAM - ^ 

L2 = NUMBER OF' VARIABLES INPARAMETRIC EQUATION SELECTED 
L3 = NUMBER OF COEFFICIENTS TO BE DETERMINED, INCLUDES BO 

call plot vector file ONLY WHEN OUTPUT INCLUDES PLOTTING - _ . 

CALL PSEUDO ' ~ 

COMPLETE DATA DECK SETUP INCLUDING OPT'iON AND CASE' CONTROL CARDS 
FOR EACH DATA SET ARE REQUIRED IF 1 -.CONTINUE CARD, IS HERE _ 

1 continue p • 

READ INPUT AND OUTPUT OPTIONS AND 
initial VALUES OF M-H AND RAB .CONST ANTS 
• ipuy = INPUT LISTING OPTIONS 

0 - NO INPUT listing, 

1 - CASE CONTROL VARIABLES * 

2 - + data set observations ■■ 

3 + transformed REGRESSION VARIABLES FOR FIRST 
FIVE observations 

OUTPUT = OUTPUT Of=T I ONS 

0 - NO RESIDUALS 

1 - residuals regressed SPACE 


00000034 

00000035 

00000036 

00000038 

00000039 

00000040 

00000041 

00000042 

00000043 

00000044 

00000045 

00000046 

00000047 

00000048 

00000049 

00000065 


00000081 


00000053 


000000S8 

00000060 

+ 





o n 


ro 

ro 


C 

C 

C 

C 

C 

C 

C 

c 

c 


c 

c 

c 

c 


c 

c 

c 

c 

c 

c 


2 - 1 + real space residuals 

3 - 2 + residual plot IN REGRESSED SPACE 
READ(5*4) iNPUTtlOUT, TA»RA 

4 FORMAT{2I5»2F10«0) 

IF(EOF»5)900f9 
9 CONTINUE 

READ CASE CONTROL CARDS 

PUT BLANK CARD AFTER LAST CASE CARD 

LLO = total number OF VARIABLES FOR CASE 

NPAM = PARAMETRIC EXPRESSION TO BE EVALUATED 

1 - LARSON-MILLER (LM) 

2 - ORR-SHERBY-DORN (OSD) 

3 - MANSON-SUCCOP (MS) 

' 4 - mAnSON-HAFERD (MH) 

. 5 - RABOTNOV (RAB) 

13=1 

3 READ (‘S»2) NPAM ( 13 ) .LLO ( I 3 ) 

IF(LL0(13)) 7.8.7 

7 13=13+1 
GO TO 3 

2 FORMAT (215) 

8 13=13-1 

ONLY ONE SETUP OF OPTION AND CASE CONTROL CARDS ARE REQUIRED 
FOR MANY DATA SETS IF 1 - CONTINUE CARD IS HERE 
1 continue 

READ DATA SET I DEN I F I CAT I ON (TYPE) 

READ ( 5 * 777 ) ( TYPE ( I ) . I = I . S ) 

777 FORMAT (8A10)‘ 

IF (EOF. 5) 900.6 
6 1 = 1 

READ IN OBSERVATIONS 

IF number of CORRELATING VARIABLES CHANGES. 

change STATEMENTS 5 AND 10 

RS(I.1)= RUPTURE TIME 

RS(I.2)= APPLIED STRESS. PSI 

RS(I.3)= TEST TEMPERATURE. DEGREES F 


(10000061 
00000062 
0000005 I 
00000052 
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00000064 

00000069 

00000070 

00000071 

00000072 

00000073 

00000074 

00000075 

00000067 

00000068 

00000076 

00000077 

00000078 

00000079 

00000080 


00000082 

00000083 

00000084 

00000085 

00000086 

00000088 


00000091 

00000092 

00000093 





/(TO 


ro 

(jO 


C PUT BLANK CARD BEHIND LAST DATA CARD 

10 READ(5»5) RS(1*1)» RS(I»3). RS(I«a) 

IF (EOF* 5 > 900,901 

901 CONTINUE 

IF (RS { I , 1 )-0. ) 11,12,11 

11 1 = 1+1 

GO TO 10 
5 FORMAT(3F12.0 ) 

12 Ll=r-1 

IF( INPUT- 1 )3dl ,300,300 
C INPUT = 1 LISTING 

300 WR1TE(6,414) 

WRITE(6,220) 

WRITE(6»221 ) 

WRITE (6, 502) (TYPE( I ) , 1 = 1 ,8) 

302 FORMAT.! 10x,*DATA SET»/10x. 8A10/) 

WRITE <6, 299 ) 

299 FORMAT!* OPTION CARD*) 

WRITE! 6, 303) INPUT, I OUT, TA,RA 

303 FORMAT!* INPUT= *,I1/* IOUT= * ,11/* TA ="* , F 1 0 . 0/* RA= *,F10.4/) 

WRITE! 6,304 ) 

304 FORMAT!* CASE CONTROL CARDS*/5X,* PARAMETER CODE*,SX, 

1*N0. COEFFICIENTS*/) 

WRITE! 6, 305) !NPAM! I ) ,LLO ! I ) , 1 = 1 , 13 > 

305 FORMAT ! lOX, 15, 15X, 15 ) 

301 CONTINUE 

C INPUT = 2 LISTING 

IF! INPUT-2) 309,308,308 
308 WRITE!6,414> 

WR ITE !6,306 ) 


00000089 
00000090 
00000094 
00000095 
00000096 
00000097 
00000098 
00000099 
00000100 
00000101 
00000102 
00000103 
00000104 
00000105 
00000106 
00000107 
00000108 
00000109 
000001 1 0 
000001 1 1 
OOOOOl 12 
000001 13 
OOOOOl 1 4 
OOOOOl 15 
OOOOOl 16 
OOOOOl 17 
OOOOOl 18 
OOOOOl 1 9 
00000120 
00000121 
+ 



( 


306 FORMAT (5X«* INPUT DATA OBSERVATIONS*/ 3X » *NO. ♦ ♦ 1 4X » * TIME** 5X* 
1 *STRESS*, 5X**TEMPERATURE* ) 

WRITE<6*307) ( I *RS( I * 1 ) *RS( I *2) *RS{ I *3 ) ♦ I = l*un 

307 FORMAT! 15* 10X«F10»2* F8.0 * 4XtF10«C) 

309 CONTINUE 

C START CASE LOOP (13) 

C 13 = NUMBER OF CASES (PARAMETRIC EQUATIONS) TO BE EXAMINED 

C FOR EACH data SET 

DO 350 KK= 1*13 
NEGS3=0 
L3=LL0 (KK ) 

L2=L3-1 

LAST=0 

BMSe= 1000000, 

XN=L1 
LIM = 0 
IFG=0 

L=NPAM (KK ) 

IF (L-A )22*Hi i20 

20 CONTINUE 

c rabotnov constants 

C=RA 

DEL =0,1 
0ELMIN=0.001 
GO TO 23 

c manson-haferd constants 

21 CONTINUE 
C=TA 

DEL = 1000,0 
DELMIN=10, 

GO TO 23 

22 LAST=2 

23 CONTINUE 


00000122 

00000123 

00000124 

00000125 

00000126 

00000127 


00000130 
00000131 
00000132 
00000133 
00000134 
00000135 
OOOOOl 36 
00000137 
00000138 
00000139 
00000140 
OOOOOl A1 
00000142 
00000143 
00000144 
00000145 
00000146 
00000147 
00000148 
00000149 

oooooiso 

00000151 

00000152 

00000153 

00000154 

+ 





57 CONTINUE 

c select parametric form for regression 

L=nPAM (KK ) 

GO TO (61 »62»63»S4»65> .L 

61 CALL LM (Y»RS<F«L1 ) 

GO TO 66 

62 call 0SD(Y.RS»F»L1 ) 

GO TO 66 

63 call MS(Y«RSiF,L1 ) 

GO TO 66 

64 CALL MH( Y»RStF»Ll»C > 

GO TO 66 

65 call RAB (Y.RS«F»L1 ♦C) 

66 CONTINUE 
SSER=0,0 

C ZERO’ A»8»SUMX1 ARRAYS 

DO 473 M=l»12 
^ DO 473 J=1.12 

A(M*J)=0.0 
B(M«J)=0.0 
473 SUMXl (M* J )=0*0 
DO 105 M=1»L2 
SUMX(M )=0,0 

D0i05l = ni-1 

105 SUMX(M ) =SUMX (M )+F ( I « M ) 

DO 106 M=1*L2 

DO 106 J=1«L2 
SUMP2(M*J)=0.0 
D0106I=1 *L1 

106 SUMP2(M*J)=SUMP2(M, J)+F< I , M ) »F ( I « J ) 
SUMY=0,0 

SUMY2=0.0 
D0107I=1 »L1 


OnOOO! 55 
00000156 
00000157 
000001 SB 
00000159 
00000160 
00000161 
00000162 
00000163 
00000164 
00000165 
00000166 
00000167 
00000168 
00000169 
00000170 
00000171 
00000172 
00000173 
00000174 
00000175 
00000176 
00000177 
000001 78 
00000179 
00000180 
00000181 
00000182 
00000183 
00000184 
00000185 
00000186 
00000187 
+ 



SUMY = SUMY + Y ( 1 ) 

107 SUMY2=SUMY2+Y { I )»*2 
DO 108 M=1«L2 
SUMXY { K ) =0. 0 
DOiORI=I »L1 

108 SUMXY ( M ) =SUMXY (M )+F ( I » M )*Y ( I > 

DO 109 M=1«L2 

DO 109 J=1 «L2 

109 SUMX 1 (M » J ) =SUMP2 < J )- ( SUMX ( M ) *SUMX ( J ) ) /XN 
DO 110 M=1 »L2 

1 I C SUMXl ( M «L3 )= SUMXY ( M > - ( SUMX ( M )*SUMY ) /XN 
SUMXl (L3«L3 > =SUMY2- (SUMY**2 ) /XN 
D04?M= 1 «L3 

42 SUMXl (L3<M )=SUMX1 (M»L3 ) 

N3=L3+1 
D016I =1 »L3 
SUMXl ( I <N3 )=0.0 
D016M=1 *L3 

16 SUMXl ( I »N3 )= SUMXl ( I «N3 ) +SUMX1 ( I . M ) 

DOl 7J-1 »N3 

17 A ( 1 , J ) =SUMX1 < 1 ♦ J) 

SUMS { 1 ) =0.0 
D018J=1 .N3 

R<1.J) = A(1.J)/A(1.1) 

IB SUMR < 1 ) =SUMR ( 1 )+P ( 1 . J 1 
SUMS ( 1 > =SUMB (1 ) -B { 1 . N3 ) 

D0115I=2.L3 
DOl 1 5J= I »N3 
NIX=I-1 
temp=o.o 
DOl 161 1=1 iNIX 

116 TEMP = -A ( I 1 . I )-X-B ( I 1 . J )+TEMP 
A( I » J)=TEMP+SUMX1 (I ♦ J) 


000001 88 
00000189 
OOOOO I 90 
000001 91 
000001 92 
00000193 
00000194 
00000195 
00000196 
00000197 
00000198 
00000199 
00000200 
00000201 
00000202 
00000203 
00000204 
00000205 
00000206 
00000207 
00000208 
00000209 
00000210 
0000021 1 
0000021 2 
00000213 
00000214 
00000215 
00000216 
00000217 
00000218 
00000219 
00000220 
+ 



^ o / 


IV) 


11'=;. B t T « J ) =A { I * J ) /^A ( T » I 5 
D029I=1 «LE 

29 SSR a )=AU «L3 )*B ( I tL3 ) 

REGSS=SUMX1 (L3,L3 )-A (L3,L3 ) 
SSER=A (L3«L3 ) 
C0RC=REGSS/SUMXI (L3»L3) 
XNl=Ll-l-2-l 


C 

C 


XMSER=SSER/XN1 

Z[P=XMSER 

STD=SQRT(XMSER> 

XMRSS=REGSS/L2 


51 


FTR.9? = XMRSS/XMSER 
T0TSUM=SUMX1 <L3»L3) 

ITERATE ON LOWEST 1 -CORC FOR RABOTNOV SOLUTION 
IF (NPAM {KK ) •EQ.S ) ZIP=1.0-CORC 

LOOP AROUND ITERATION FOR L-M»0-S-D» AND M-=S SOLUTIONS 
L=NPAMCKK) 

IF(L«ET«4) go TO 1439 

IF(LAST-I) 51t52,1439 ■ . 

CONTINUE 

call ITERIC ,CBEST*ZIP , BMSE . I FG . I CT , DEL , DELM I N , L I M . LAST ) 

GO TO 57 

CONTINUE 

C=CBEST 


LAST=2 
GO TO 57 
1439 CONTINUE 

D0117I=2«L2 
SUMA { I )=0«6 
SUMB < I ) =0*0 


D01l7J=l»L3 
SUMA ( I ) =SUMA < n+A ( I « J ) 
1 1 7 SUMB ( 1 ) =SUMB ( 1 ) +B ( I » J ) 


nonoo??i 

00000222 

00000223 

00000224 

00000225 

00000226 

00000227 

00000228 

00000229 

00000230 

00000231 

00000232 

00000233 

00000234 

00000235 

00000236 

00000237 

00000238 

00000239 

00000240 

00000241 

00000242 

00000243 

00000244 

00000245 

00000246 

00000247 

00000246 

00000249 

00000250 

00000251 

00000252 

00000253 

*4- 





00801=1 *L2 
D08OJ=l *L2 

80 D( I « J)=SUMX1 ( 1 » Jl 
DD( 1 ♦ 1 ) =L1 

0072 M=1 ,L2 
I =M+1 

72 DD( I ♦ 1 ) =SUf^X (M ) 

0073 K=l ,L2 
J = K + 1 

73 DO ( 1 ♦ J )=SUMX(K ) 

0074 M= 1 tL2 

I =M+ 1 

0074 K= 1 4L2 
J = K+1 

74 00 ( I < 3 1 = O (M »K ) 

C CHANGE L9 when YOU RED I MENS I ON PROGRAM 
L9=ll 

CALL MAT I NV ( L9 »L3 » DO . 0 ♦ G* 1 » OETERM . I SCALE » I P 1 VOT « I NOEX ) 
CO 

0081 1=1 »L2 

81 E ( I « 1 1 =SUMX1 ( I tL3 } 

CALL MATINV (L9,L2«D , 1 , E . 1 » DETERM , I SCALE . I P I VOT ♦ I NDEX ) 
PAR ( 1 )=B(l2.L3) 

M3 = L2 
K3=L2 

D0113I=2»L2 
MIX=I-l 
M3=M3-1 
TEMP1=0.0 
001141 1=1 »M1X 

TEMPI =”PAR( I 1 )*B(M3«K3)+TEMP1 
1 14 K3=K3“1 

PAR < I )=T£MP1+B <K3»L3 ) 

1 13 K3=l2 


00000254 

00000255 

^0000256 

00000257 

OOOOOP58 

00000259 

00000260 

00000261 

00000262 

00000263 

00000264 

00000265 

00000266 

00000267 

00000268 

00000269 

00000270 

00000271 

00000272 

00000273 

00000274 

00000275 

00000,276 

00000277 

00000278 

00000279 

00000280 

00000281 

00000282 

00000283 

00000284 

00000285 

00000286 

+ 



D047Ni = l »L2 

47 SUMX2(M)=SLIMX(M )/XN 
SUMY3=SUMY/XN 
K4=l2 

D0216I =1 »L2 
PARI (K4 ) =PAR ( I ) 

216 K4=K4-1 
PARO=0. 0 
D0217I=1 «L2 

217 PARO=PARO-PAR1 ( I ) *SUMX2 < I ) 

PAR0=PAR0+SUMY3 

N5=l1-L2-1 

XNl =N5 

XMSER=SSER/XN1 

IF(XMSER .GT. 9.0E+100) GO TO 350 
00491 = 1 .L2 

1F(D(I»I)»I-’T*0,0)WR1TE<6.1100)(I*D(I»1)) 

IF(D( I )«LT.O.O) NEGSB=NEGSB+1 
C AVOID '^ODE 2 DUMP ABORT CASE 3/3/76 

IF (D ( I » I ) *LT. 0*0) GO TO 350 

1100 FORMAT (//«5X*****#NEGAT I VE SB ( I ) * I =* ♦ I 3 . *0 I = * » E20 • 8 ) 
49 SB { I )=SQRT(ABS(D( I . 1 ) -»XMSER ) ) 

DOl 18I=1 *L2 
T < I )=PARI ( I )/SB ( 1 ) 

118 T< ! )= ABS(T( I )’) 

C CY(I>= SOLUTION IN REGRESSION SPACE 

D0122l = l *L1. 

SUMCY =0-.0 
TEM=0.0 
•DOl 23M= 1 »L2 

123 TEM = TEM4-PAR1 (M)*F(I *M) 

SUMCY =SUMCY+TEM 
122 CY(I)=SUMCY -fPARO 

' 4 . 


00000287 
00000288 
00000289 
00000290 
00000291 
00000292 
00000292 
00000294 
00000295 
00000296 
00000297 
00000298 
00000299 
00000300 
00000301 
00000302 
00000303 
00000304 
00000305 
00000306 
00000307 
00000308 
00000309 
0000031 0 
0000031 1 
00000312 
00000313 

00000314 

00000315 

00000316 

00000317 

00000318 

00000319 

+ 
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C 


calculate min, max, range, mean 

call MINMAXIYMIN, YMAX, YRAN,YMEAN,Y,L1 ) 

DO 95 I =1 ,L2 
DO 96 K=1 »L1 
96 X(K)=F(K, I ) 

CALL MINMAX(XMIN( 1 ) , XMAX(I)« XRAn ( 1 ) < XME AN ( I ) » X i H ) 

95 CONTINUE 

PRINT REGRESSION STATISTICS 
WRITE(6,4I4 ) 

WRITE (6, 220 ) 

WRITE<6,22I ) 

220 FORMAT! 5X,* LEAST-SQUARES REGRESSION PROGRAM FOR PARAMETRIC*) 

221 FORMAT! 5X,* ANALYSIS OF CREEP-RUPTURE DATA*) 

WRITE (6,526 ) 

526 FORMAT! /* REGRESSION VALUES ■--* ) 

WRITE (6,502 ) (TYPE ( I ) , I = 1 , 8 ) 

502 FORMAT!* DATA SET *,1X,8A10) 

MM = NPAM (KK ) 


00000320 

00000321 

00000322 

00000323 

00000324 

00000325 

00000326 

00000327 

00000328 

00000329 

00000330 

00000331 

00000332 

00000333 

00000334 

00000335 

00000336 

00000337 


write (6,506 )PAM (MM ) 

506 FORMAT!* PARAMETER SELECTED 

WRITE(6,507> ( Ll ) • ' 

507 FORMAT!* N0« OF OBSERVATIONS 
WRITE!6,508) ! L2) 

508 FORMAT!* N0» OF INDEPENDENT VARIABLES 
WRITE(6,509) (XNl ) 

509 FORMAT!* RESIDUAL DEGREES OF FREEDOM 
WRITE!6i514) ! FTRSS ) 

514 FORMAT!* F - VALUE 
WR I TE ( 6 , 545 > ! XMSER ) 

545 FORMAT!* RESIDUAL MEAN SQUARE 
WRITE!6,546) ! STD) 

546 FORMAT!* STANDARD ERROR 


* , 

7X, AS) 

00000338 

00000339 

*, 

8x, 

14 ) 

00000340 

00000341 

* , 

8X, 

14) , 

00000342 

00000343 

*, 

8X, F4) 

00000344 

00000345 

*, 

FI 2. 1 ) 

00000346 

00000347 

*, 

E12.4 ) 

00000348 

00000349 

*, 

E 1 2 . 4 ) 

00000350 

00000351 


WRITE(6,547) ! SSER) 


00000352 


+ 





(JO 

H 


c 


c 


547 format {* RFSIDUAL SUM OF SOUARFS 
write (6 » 548)(T0TSUM) 

548 FORMAT (^<- TOTAL SUMS OF SQUARES 

WRITE (5* 549)(C0RC) 

549 format (* MULT* CORREL* COEF . SQUARED 
WRITE(6»320 ) 

320 FORMAT </) 

IF(MM.EQ.4) WRITE<6»432) CBEST 
IF(MM*EQ*5) WRITE(6.433) CBEST 

432 format (* MANSON - HAFERD CONSTANT <TA) 

433 FORMAT (* RABOTNOV CONSTANT ( RA ) 
WRITE(6t492) ( YMI N < YM AX ♦ YRAN « VMEAN ) 

492 F0RMAT(5 X** MIN Y =■)<■» E 1 1 . 2 , 3X ♦ * MAX Y 
1 E11*2«3X* MEAN Y =*»EIl*2/> 

INPUT= 3 LISTING 
IF{ INPUT-3 131 1 »310»31 1 
31C WRITE(6*312) 

312 FORMAT!/ 5X«*FIRST 5 OBSERVATIONS - 
15X»*Y* ♦ 18X*» XI - X(l2) *) 

DO 313 1=1 *5 

DO 315 J=1»L2 
315 TEMP! J)=F ( I ♦ J) 

WR1TE(6*314) (Y<I )< (TEMP(J)*J=1»L2>) 

313 CONTINUE 

314 F0RMAT(8E15.5) 

WRITE(6t320 ) 

311 CONTINUE 

WRITE(6»422) 

422 FORMAT (3Xi* 

1 4X» *T*» 

2*»3Xt »RAN X( I 
WRITE(6.535) (PARO ) 

535 format (6X»*0*, 1 IX «E14.4 ) 


00000353 
00000354 
00000355 
00000356 
00000357 
00000358 
00000359 
00000360 
00000361 
00000362 
00000363 
00000364 

2,3X** Y RANGE =*. 00000365 

00000366 
00000367 
00000368 
00000369 

transformed VARIABLES*/ 0000037Q 

00000371 
00000372 
00000373 
00000374 
00000375 
00000376 
00000377 
00000378 
00000379 
00000380 
*S.E«COEF.*« 00000381 
3X**MAX X( I >00000382 
00000383 
00000384 
00000385 
+ 


*, F12.4) 

*0 E12.4) 
*, F12.4/) 

=*.F10.1/> 
=*,FI0.5/ ) 

= *,E1 1 


I *»2X»*VAR1 ABLE* « 4X 4 *COEF.P ( I ) * . 3X » 
5X4 *MEAN X(I>*4 3X4 *MIN X(I)*4 

4X4*RINF*4 3X4 #PSUM*4 3X 4 *CERR* > 



LO 

ro 


'S3 


M=NPAM (KK )*6-6 
DO 420 1=1 »L2 

CERR ( I 1 = 1 00.0* ( (PARI ( I )-E ( I ) )/PARl < I ) ) 

RINF= (XRAN C I )*PAR1 ( I 1 + 1 .OE-30 ) /YRAN 
SSRR = SSR( I l/REGSS 

WRITE<6*421 ) ( I .VAR( I +M ) »PAR1 (lit SB ( 1 ) » T C I ) .XMEANC I ) ♦ 

1 XMIN( I ) »XMAX ( I ) .XRAN ( I ) » R I NF « SSRR ♦ CERR ( T ) ) 

420 CONTINUE 

421 FORMAT! 17. 4X.A8.1X. E12.4. Ell»2. F7 • 2 » 

1 E13.3. 3E11.2. F3.2. F7.3. F7.2) 

WRITE(6.424) 

424 FORMAT! /* VARIABLE CODE*/ 1 OX t*S=LOG STRESS*/ 1 OX. *T=TEMPERATURE* 
1/1 0x.*DT=T-TA*/1 OX.*L=STRESS*/ ) 

CUMERR=0.0 
RRI S=0.0 
SRIS2=0.0 
EMAx=0.0 
EMAxP=0 .0 
NZER0=0 
SDP=0,0 
SSDP=0*0 
DPMAX = -10.0 
XDF=ALOG10(XN1 > 

T6=0. 86 1 8559 -0.98427 1 5*XDF+0 • 5849466*x0F**2- • 1 1 59365*XDF**3 
T6=l 0.**T6 

CCCC **** START 333 LOOP***** 

D0333M=1 «L1 
X! 1 1=1.0 
DO 92 K=1 .L2 
I=K+1 

92 X! I )=F !M.K1 

DO lOO J =1.L3 
TEMP!J) =0.0 


00000386 
00000387 
00000388 
00000389 
00000390 
00000391 
00000392 
00000393 
00000394 
00000395 
00000396 
00000397 
00000398 
00000399 
00000400 
00000401 
00000402 
00000403 
00000404 
00000405 
00000406 - 
00000^07 
00000408 
00000409 
0000041 0 
0000041 1 . 
00000412 
00000413 
00000414 
00000415 
00000416 
00000417 
00000418 
+ 



Q // 


(jO 

uo 


no 1 00 I =1 fLO 

TEMP(J)= TEMP(J>+ X(I>* DD(I*J) 

100 CONTINUE 
ANS =0,0 
DO 200 J=1 fL3 

ANS= ANS +TP’MP( J)*X( J) 

200 CONTINUE 

XmfR = XMSrR 
x:(vie:r= abs(xmer) 

ANS = ABS(ANS) 

C CALCULATE 95 PERCENT STATISTICAL 

DELTA=T6*SQRT ( XMER^ANS ) 

Cl MAX ( M )=CY ( M )+DELTA 
CIMIN (M ) =CY(M J-DELTA 
DELTA = T6*SQRT <XMER* ( 1 +ANS ) ) 

PYMAX(M }=CY (M)+DELTA 
PYMIN(M)=CY(M)-DELTA 
RI5(M)= CYIM)- Y(M) 

C AVOID DUMP WHEN Y=0 2/25/76 

IF ( Y (M ) •EQ,6»0 > Y (M ) =0 • 00000 1 
ERRPER ( M ) =R I S ( M ) /Y ( M ) * 1 00 • 

IF( ABS(RIS(M) ) *GT. ABS ( EMAX ) ) EM AX 
RIS2 =RIS(M}**2 
SRIS2=SRIS2+RIS2 
RRIS = RRIS+ABS<RIS(M) ) 
IF(ABS(ERRPER<M) ) .GT. ABS(EMAXP) ) 
CUMERR=CUMERR+ABS (ERRPER (M ) > 

333 CONTINUE 

CCCC ****** END 333 LOOOP SSSS** 

C FIND OBSERVATIONS OUTSIDE 

401 FORMAT</* OBSERVATIONS OUTSIDE OF 
15X«*0BS.*, 5X«*CALC Y*» 5X«*PYMIN*» 
I3AD=0 


noonna l o 
00000420 
0000042 1 
00000422 
00000423 
00000424 
00000425 
nn0004?6 
00000427 
00000428 

intervals 00000429 

00000430 
00000431 
00000432 
00000433 
00000434 
00000435 
00000437 
00000438 
00000439 
00000440 

= RIS(M> 00000441 

00000442 


00000443 
00000444 

FMAXP = FRRPERIM) 00000445 

00000446 
00000447 
00000448 

OF 95 PERCENT PREDICTION INTERVAL 00000449 

95 percent PREDICTION 1 NTERVAL*/000 00450 

5X«*PYMAX* /) 00000451 

00000452 
+ 



c 


c 

c 


399 

402 

400 

403 


404 


405 


406 


407 

408 


II 

O 

o 

o 

o 


00000453 

IF <CY< I )-PYMlN( I ) ) 402«399.399 


00000454 

CONTINUE 


00000455 

I F (CY ( I I-PYMAX ( I > ) 400»400«402 


00000456 

IBAD= IBAD+1 


00000457 

IF ( IBAO.EO. 1 ) WRIT£(6*40U 


00000458 

WR 1 TE (6 *403 ) ( I» CY ( I ) ♦ PYM I N ( I ) * PYMAX ( I ) > 


00000459 

CONTINUE 


00000460 

format { I5.3E16.6) 


00000461 

qEXERMINE dp statistics in terms of log time to RUPTURE 


00000462 

DPSUM=0 «0 


00000463 

dpmax= -IOO.O 


00000464 

if (NPAM <KK ) -5 ) 404,406^404 


00000465 

DO 405 I =1 ,L1 


00000466 

DP=PYMAX( I )-PYMlN( I ) 


00000467 

IF(DP»GT»DPMAx ) DPMAX=DP 


00000468 

DPSUM= DPSUM+DP 


00000469 

CONTINUE 


00000470 

GO TO 408 


00000471 

rabotnov dp 


00000472 

DO 407 I =1 »L1 


00000473 

AVOID NEGATIVE PY DUMP 


00000474 

IF(PYMIN( I ).LT*0.0) PYMI N( I ) =1 »0 


00000475 

1 F ( PYMAX ( I ) «LT» 0 • 0 ) PYMAx ( I ) = 1 • 0 


00000476 

RPl =PYMAX ( ( 1 ,0/CPFST ) 


00000477 

rP2=PYMIN( I )**{ 1 .0/C8EST) 


00000478 

DP= ALOG10(RP1 )-AL0G10(RP2) 


00000479' 

IF (DP*GT.DPMAx >dpmax=dp 


00000480 

DPSUM=DPSUM+DP 


00000481 

CONTINUE 


00000482 

DPAVE = DPSUM/LI 


00000483 

RPl =1 0.**DPAVE 


00000484 

RP2=10,**DPMAX 


00000485 


+ 
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WRITE (6. 41 0 ) 

410 FORMAT! 5X»* 95 PERCENT PREDICTION INTERVAL S T AT I ST I CS*/25X » 
l*LOG TIME** 10 x«*REAL TIME FACTOR (ANTILOO WIDTH)*/) 

WR ITE ( 6 * 409 ) < DPAVEfRPl «DPMAX «RP2 ) 

409 FORMAT!* AVERAGE WIDTH *, 5X* F10«3*19X* F10ol/>^ MAXIMUM WIDTH** 

1 6X* FI 0.3. 19X*FI 0* 1 ) 

C order residuals - LARGEST TO SMALLEST 

LINEG =-Ll 
DO ?1 00 1 =1 *L1 

TEMP! I ) =0.0 

2100 TEMP!I)= A8S!RIS!I)) 

CALL A0RDER!TEMP * LINEG, I PERM) 

no I ?02 T =1 ,L1 

J= IPERM ! I ) 

TEMP! J) = I 
1202 CONTINUE 

OUTPUT = 1 OR GREATER 
RESIDUALS IN REGRESSED SPACE 
IFIIOUT -1) 413,412,412 

412 CONTINUE 
WRITE!6*414> 

WRITE! 6* 415) 

414 FORMAT! IHl ) 

415 FORMAT!* RESIDUALS - REGRESSION SPACE*/) 

WRITEI6.416) 

WRITEI6.417 ) ! I , Y ! I ) *CY ! I ) *RI S! I > .ERRPER! 1 ) , TEMP! I ) , 

1 Cl MAX! I ), Cl MIN! I ).PYMIN! I ) .PYmAX! I ) * 1 = 1 *L1 ) 

417 F0RMAT!’I5* 1X.3E12.3. ’ ,Fl0.1i 8X,F5, 4X* 4E12.3) 

416 FORMAT! 2X**OBS*. 7X.*Y 03S* , 6X * *Y CALC*.5X. *RES I DUAL* *6X * 

1 *PCTERR *♦ 1 X ,*ORD£R* *7X **CIMIN*, 7X ♦ *C I MAX* * 7X ♦ 

2 *PIMIN*,7X ,*PIKAX*/> 

413 CONTINUE 
XMRSS=REGSS/L2 
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FTRSS=XMRSS/XMSER 
DPAvE= SDP/L I 

qPSIG=(L1*SSDP-SDP**2>/(LI*(L1-1 .0) ) 

DPSIG=SQRT (DPS IG) 

STD=SQRT(XMSER) 

PLOTTING ROUTINE 

PLOT RESIDUALS WITH VARIAN ON LINE PLOTTER 
IF(IOUT -3)445»440.4A0 
440 CONTINUE 

IN( 1 )= 5HPARAM 
1N(2)= 4HPL0T 
N=L1 

ISYMD=12 
IEC=l 

CALL MINMAX(YL»YH»YRAN»YMEAN»RIS*L1 ) 

YL=1 .8*YL 
^ YH=i.8»YH 

XL=0.0 
XH=0.0 
NXM= I 
NYM=1- 

YN0TE= lOH RESIDUAL 
XN0TE5 = 10H2P NORMAL 
XN0TE6= I OH FITTED Y 

call VDIPLT(IEC« in* N, CY(l)f RIS(1>« XL » XH * YL » YH . NXM , 
1 XN0TE6* NYM» YNOTE« ISYMD) 
call A0RDER(RIS«L1 » IPERM) 

DO 430 I =1 .LI 
j=1PERM ( I ) 

TEMPI I )= RIS(J) 

XI = I 

FZ= (Xl-»375 )/ (LI +*25 > 

IF(FZ-0.5)570.570.571 
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570 XX=ALOGl 0(FZ ) 

ZP(I)= 1.04505 + 4.35979*XX + 3 . 46057*XX*-*-2+ 1 « 90878*XX**3 
1 + 0.54456*XX**4+ 0.0608*XX**5 
GO TO 572 

571 XX=ALOG10< 1 .0-F2 ) 

2P(I)= 1.04505 + 4.35979*XX + 3 . 46057*XX**2+ 1 . 90878-K-XX**3 
1 + 0.54456*XX**4+ 0.0608*XX**5 
ZP ( I ) =-ZP ( I ) 

572 CONTINUE 
430 CONTINUE 

YL=0.0 
YH=0.0 

call VDIPLT(IEC. IN, N, ZP ( 1 ) . TEMP (1 ) , XL . XH , YL « YH . NXM , 

1 XN0TE5, NYM, YNOTE, I SYMD ) 

445 CONTINUE 

OUTPUT. = 2 OR GREATER 
real SPACE RESIDUAL OUTPUT 

BACKTRANSFORM SOLUTION AND PREDICTION INTERVALS 
MX=NPAM (KK ) 

DO 441 M=1.L1 

GO t 0<201 ,201 ,201 ,201 i203) , MX 
201 CY(M)= 10.0**Cy(M) 

PYMAX'I M )=10.0-lHfPYMAx(M ) 

PYMINIM )=10.0**PYMIN(M) 

C I MAX’! M ) = 1 0 . 0**C I MAX < M ) 

CIMIN(M)=10.0**CIMIN(M) 

GO TO 205 
203 CONTINUE 
C AVOID NEGATIVE TO A POWER DUMP 

IFIPYMINIM) .LE. 0.0) PYMIN(M)=1.0 
IF(PYMAX(M).LE. 0.0) PYMAX(M)=1.0 
IF (Cl MAXIM ) .LE. 0.0 ) CIMAX(M) = 1 .0 
IF(CIMIN(M).LE.0,0) CIMIN(M)s1.0 
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IF {CY ( M ) .LE . O .0 ) CY{M) = 1.0 

CY(M)= CY(M)**(1 •0/C8EST) 

PYMAX(M)= PYMAX <M )** ( 1 . 0/CBEST ) 

PYMIN(M)= PYMlN(M>-tt*( 1 *0/CBEST) 

Cl MAX<M )=C IMAX(M ) *•» ( 1 •0/CBEST ) 

CIMIN(M)=CIM1N(M)#*(1 .0/CBEST > 

205 continue 

R1S(M)= RStM.l )-CY(M) 

ERRPER <M)=(RIS(M)/RS(M,1 )>*100. 

441 CONTINUE 

C ORDER real SPACE RESIDUALS 

DO 425 I =1 .L 1 
TEMP (I > =0.0 

425 TEMP ( I ) =ABS (ERRPER ( I ) ) 

CALL AORDER (TEMP.Ll NEC, I PERM ) 

DO 1203 1=1 .LI 

J= IPERM ( 1 ) 

TEMP(J)=I 
1203 CONTINUE 

IF( IOUT.LT.2 ) GO TO 350 
V/RIT= (6.414 ) 

WR1TE(6.431 ) 

431 FORMAT!* 8ACKTRANSF0RMED RESIDUALS - REAL SPACE*/) 
WRITE (6. 416) 

WR ItE ( 6. 41 7 ) ( I ,RS ( I . 1 ) . C Y ( I ) i R 1 S ( I ) .ERRPER! I ) «TEMP( I ) . 
1 C IM IN ( I ) .Cl MAX ( I ) . PYMIN C I )’-tPYMAX ( I ) . 1 =1 .Ll ) 

350 CONTINUE 
GO TO 1 
900 CONTINUE 

C CALL CALPLT ROUTINE ONLY WHEN PLOTTING 

IFdOUT .GE.3) CALL CALPLT ( 0 . 0 . 999 ) 

STOP 

END 

SUBROUTINE AORDER (AA.N. IPERM) 
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C THIS SUBROUTINE ORDERS VALUES IN AA AND STORES ORDER IN I PERM 

C N IS NUMBER OF VALUES IN AA 

C I PERM IS ORDERED WITH RESPECT TO LOCATION OF VALUFS IN AA 

C IF N IS POSITIVE IPERMCl) HAS LOCATION IN AA OF SMALLEST VALUE 

C I PERM (N) HAS LOCATION OF LARGEST VALUE IN AA 

C IF N IS NEGATIVE IPERM IS ORDERED BY LOCATION OF LARGEST TO 

C SMALLEST VALUES IN AA 

C ARRAY A A IS NOT CHANGED 

DIMENSION AA(1>» IPERM(l) 

LOGICAL SWITCH 
nabs = IABSIN) 

DO I 00 1 = 1 «NABS 

1 00 IPERM ( I ) = I 

IF< NABS .LT.2) RETURN 
200 SWITCH = .false. 

00500 1= 2.NABS 
I 1= IPERM ( i-n 
JJ= I PERM! I) 

IF( N.LT. 0) GO TO 400 
1F(AA( I n.LE.AAI JJ> ) GO TO 500 
300 ITEMP= IPERM ( I- 1) 

IPERM<I-1> = IPERM (I) 

I PERM ( I ) = ITEMP 
SWITCH = .TRUE. 

GO TO 500 

400 IF{ AA { I I ) .LT.AA { JJ ) ) GO TO 300 
500 continue 

IF( SWITCH) GO TO 200 
900 RETURN 
END 

SUBROUTINE LM (Y.RS.FiLl) 

CONVERTS TIME. STRESS. AND TEMPERATURE TO FORMAT REQUIRED 
FOR LINEAR SOLUTION OF LARSON-M I LLER EXPRESSION 
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C SOLUTION ALLOWS FIFTH ORDER EXPANSION OF LOG STRESS 

C Y= BO + B1 (XI )+Ba CX2 ) e6CX6) 

C WHERE Y = LOG (RUPTURE TIME) 

C S = APPLIED STRESS IN PSI 

C T = TEST TEMPERATURE IN DEGREES F 

C Xl= l/(T+460> 

C X2= LOG(S >/ (T+460 ) = S/TK 

C X3= S**2/TK 

C X4= S**3/TK 

C X5= S**4XTK 

C X6= S**5/TK 

C C«B0-B6 = CONSTANTS DETERMINED BY LINEAR LEAST SQUARES METHOD 

C B0= OPTIMUM L-M CONSTANT (C) 

DIMENSION Y(200)« RS(200,F), FC20n,10) 

DO 10 1=1. LI 

Y( I )= AL0G10(RS( I . 1 ) ) 

S= ALOG10(RS( 1.2)) 

T= (RS( I .3)+460.0 ) 

F ( I ♦ 1 ) = 1 .0/T 

F ( I ,2 ) = S/T 
F( 1.3) = S**2/T 
F(1.4) = S»*3/T 

F( I ,5) = S**4/T 
F(1.6) = S**5/T 
10 CONTINUE 
RETURN 
END 

SUBROUTINE OSD ( Y . RS . F . L 1 > 

C CONVERTS TIME. STRESS. and TEMPERATURE TO FORMAT REQUIRED 

C FOR LINEAR SOLUTION OF ORR-SHERB Y-OORN EXPRESSION 

C SOLUTION ALLOWS FIFTH ORDER EXPANSION OF. LN STRESS 

C where Y = LOG (TIME TO CREEP EVENT) 

C S = APPLIED STRESS I N PS I 
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C 

r 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 
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T = TEST TEMP I N DEGREES F 

XI = l/TK 
X2= LN(S> = SL 
X3~ SL**2 
X4= SL**3 
X5= SL**4 
X6= SL**5 

B0-B6= CONSTANTS, DETERMINED BY LINEAR LEAST SQUARES METHOD 
Bl= DELH/R 

DELH= APPARENT ACTIVATION ENERGY 
R= UNIVERSAL GAS CONSTANT^ 

DIMENSION Y(200>» RS(200»51» F(200*10) 

DO lO I=1«L1 

Y( I )= ALOGlOiRSn »1 ) ) 

S= AL0G10<RS< 1*2)) 

T= {5./9. )*{RS( I *3)-32. ) +273. 

F f I , 1 )= 1 .0/ T 

F < T * 2 ) = S 
F(I.3)= 

F ( I » 4 ) = S**3 

F ( I , 5 ) = S**4 

FCI,6)= S**5 

10 CONTINUE 
RETURN 
END 

SUBROUTINE MS ( Y ♦ RS . F * L 1 ) 

converts TIME.STRESS.AND temperature TO FORMAT REQUIRED 
FOR LINEAR SOLUTION OF MANSON-SUCCOP EXPRESSION 
SOLUTION ALLOWS FIFTH ORDER EXPANSION OF LOG SfRESS 
B1 = OPTIMUM M-S constant (C) 

DIMENSION Y(200). RS (200.5 >« F(200.l0) 

DO 10 I =1 »L1 

Y< I )= ALOGIO(RS( I .1 ) > 
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S= ALOGIO (RS( 1 ) ) 
T= R5( I «3) 

F( I ,1 )=T 


F ( U2 ) = S 
F(I,3)= S**2 

F(I,4)= S**3 

F(I,5)= S**4 

F(I,5)= S**5 


10 CONTINUE 
RETURN 

end 

SUBROUTINE MH( Y*RS^F »L 1 »CMH) 

C FOR NONLINEAR SOLUTION OF MANSON-HAFERD EXPRESSION 

C CMH = TEMPERATURE OFFSET (TA) 

C Xl= T-CMH =DT 

C X2= DT*S 

C • X3= 0T*S**2 

C X4= DT*S**3 

C X5= DT*S**4 

C X6= DT*S**5 

C SOLUTION IS ITERATED TO FIND CMH WHICH PRODUCES BEST FIT 

DIMENSION Y(200)» RS (200,5 )» F(200,10) 

•DO 10 1=1, LI 

Y ( I )= ALOGl 0 (RS( 1,11) 

S= ALOGl 0 <RS ( I ,2 1 ) 

DT= RS( I ,3)-CMH 
F ( I , n= DT 
F( I ,2)= DT*S 
F(I,3)= DT*S*-»^2 
F(l,4)= DT*S**3 
P(I,5)= DT*S**4 
F(I,6)= DT*S**5 
10 CONTINUE 
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RETURN ■ ^ * 
END 000007*02 

SUBROUTINE ITER(C« BC* X ♦ PBX< I FG * I CT < OEL» DELMIN» LIM*lAST) 00000753 
ITERATES CONSTANT <C> TO MINIMIZE VALUE (X) 0000075A 
BC = VALUE OF CONSTANT ASSOCIATED WITH LOWF.ST ( BEST ) X VALUE 00000755 
PBX= BEST PREVIOUS VALUE OF X 00000756 
IFG =FLAG to control INCREASING OR DECREASING C FOR NEXT I TERAT I ON00000757 
DEC ^CONTROLS SIZE OF C INCREMENT 00000758 


VALUE 


ICT =ALLOWS C TO INCREMENT BEYOND BC BEFORE CHANGING 

INCREMENT SIZE. ITERATION STOPS WHEN DEL .LE. DELM I N 
LIM = COUNTER FOR ITERATIONS 
LAST=END ITERATION FLAG 
1F( IFG) 5.5.30 
INCREASING C . 

5 LIM = LIM+1 

NEXT CARD prevents NEGATIVE X FROM BEING BEST X VALUE 
IF(X*LE. 0.0) GO TO 10 
IF(PBX .GT. X) 8C=:C 
IF(P8X .GT. X) PBX=X 
IFIPBX .EQ. X) ICT=0 
1F(X .OT. PBX) ICT=ICT+1 
IFIDEL .LE.DELMIN) GO TO 40 
1F( ICT.LT. 2) GO TO 10 
IF( ICT.LT. 10) GO TO lO 
DEL = 0.3*DEL 
ICT = 0 
C= C-DEL 
C=BC+5. 0*DEL 
IFG=1 
GO TO 50 
10 C= C+OEL 
GO TO 50 

DECREASING CONSTANT 
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30 LIM =LIM+l 

C NEXT CARD PREVENTS NEGATIVE X FROM BEING BFST X VALUF 

IFIX.LE. 0.0) GO TO 35 
IFIPBX .GT. X) BC=C 
IFCPBX *GT. X) PBX=X 
IF(P3X .EQo X) 1CT=0 
IF(X .GT. PBX) ICT=ICT+1 
IF(DEL.LE.DELMIN) go TO 40 
IF( ICT.lt. 10) GO TO 35 
DEL=0.3*DEL 
ICT = 0 

C = BC-5.XH^DEL 
IFG=0 

35 C=C-DEL 
GO TO 50 
40 LAST=1 
50 CONTINUE 
RETURN 
„ END 

SUBROUTINE RAB ( Y . RS » F » L 1 . A ) 

C FOR NONLINEAR SOLUTION OF RABOTNOV EXPRESSION 

C SOLUTION ALLOWS FIFTH ORDER EXPANSION OF TEMPERATURE FUNCTION 

c where y= 'rupture time **A 

C' T= TEST temperature IN DEGREES F 

C Xl= 1/ST 

C A= ITERATED CONSTANT 

C S= STRESS IN PSI 

DIMENSION Y<200). RS ( 200.5 >» F (200. 10) 

DO 10 i =l .LI 

Y( I ) = (RSI I . 1 ) )*»A 

S = RS( I .2 ) 

T = RS( I .3) 

F(l,l)= 1.0/(S*T) 
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r(I,2)= i.0/(S*T**2) 

F(J,3)= l./(S*T**3) 

F(I,4)= l./(S*T**4) 

F(I.5)= 1 ./ <S*T**5 > 

F(I,6)= 1*/(S*T**6> 

10 CONTINUE 
RETURN 
END 

SUBROUTINE M INMAX ( CM IN» CMAX » CRAN * CMEAN* C .N ) 
CALCULATES MINIMUM* MAX I MUM * RANGE ♦ AND MEAN OF C(I) 
WHERE N= NUMBER OF OBSERVATIONS 
dimension C ( 1 ) 

CMAx=- 1 -OE+l 00 
CMIN=1 *0E+10O 
CSUM=0*0 
DO 5, I = 1 »N 
CSUM=CSUM+C ( I ) 

IF< C<I)-CMIN) 2*3*3 

2 CMIN = C( I ) • 

3 IF( C( I )-CMAX)5\5*4 

4 CMAX=C( I ) ’ 

5 CONTINUE 
CMEAN=CSUM/N ' 

, CRAN=CMAx-CMIN 
7 CONTINUE 

RETURN ‘ 

END 
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APPENDIX B 


LAIGUACE ! 
PURPOSE! 

USB! 


LANGLEY RESEARCH CENTER SYSTEM SUBROUTINES 


SUBROUTINE MATINV 
* * 


FORTRAN 

> * ^ 

To invert a real s quare ’ matrix A, solve 'the matrix 

equation AX = B, where B is a matrix of constant vectors., 

and by an opti'on evaluate the determinant. 

’ w ' V ' V ' 

• % 

’ f * r 

CALL MATINV (MAX , N, A , M ,B , lOP , t)ETERM' , ISCALE , IPI VOT , INK ) 

i 

V 

t 

f 

MAX An input integer specifying the maximum 

order of A as stated in the dimension 
statement of the calling program. 

N An input integer specifying the order 

of A; 1 < N <_ MAX. 

A An input /output two-dimensional array of the 

coefficients. On return to the calling 
program, A ^ is stored in A. A must be 
dimensioned in the calling program with 
first dimension MAX and second dimension 
at least N. The original A matrix is 
destroyed. 

M An input integer specifying the number of 

column vectors in B. M = 0 signals that 
the subroutine is used solely for inversion; 
however, in the call statement an entry 
corresponding to B must be present . 
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lOP 


DETERM 


ISCALE 


An input /output tvo-dimensional array of the 
constant vectors. On return to the calling 
program, the solution X is stored in B.. 

B should have its first dimension MAX and 
its second dimension at least M. The 
original B matrix is destroyed. 

Compute the determinant option. 

lOP = 0, Compute the determinant. 

lOP = 1, Do not compute the determinant. 

For lOP = 0, in conjunction vith ISCALE, 
represents the value of the determinant of 
A as follows: 


DET(A) = (DETERM)10^°°^^®^^^^ 


For I0P=1, the determinant is set to 1. The 
determinant is set to zero for a singular 
matrix, for both lOP = 0 or 1 option. Upon 
return from MATIBV, DETERM should he tested 
or written out in the calling program. 

(See Other Coding Information) 

For lOP = 0, the scale factor is computed hy 
the- subroutine to avoid overflow or under- 
flow in the computation of the quantity, 
DETERM. For lOP == 1, ISCALE may be a 
dianmy -argument . 



METHOD: 


ACCURACY : 


REFERENCE : 
STORAGE; 


IPIVOT , ■ A one-dimensional array used by the subprogram 
’ » to store pivotal information. It should be 

dimensioned at least N. In general the user 
does not need to make use of this array. 

IWK An integer array of temporary storage, 

dimensioned at least 2 x N. 

Jordan's method is used to reduce a matrix A to the identity 

matrix I through a succession of elementary transformations: 

1,1 , . . . , JL, . A = I. If these transformations are 

n n-1 1 

simultaneously applied to I and to a matrix B of constant 
vectors, the results are A ^ and X where AX = B. Each 
transformation is selected so that the largest element is 
used in the pivotal position. 


Total pivotal strategy is used to minimize the rounding 
-errors; however, the- accuracy of the final results depends 
upon how well-conditioned the original matrix is. A return 
with DETERM ^ 0 does not guarantee accuracy in the 
solutions or inverse. 

Fox, L., An Introduction to Numerical Linear Algebra. 

Oxford University Press, New York, 1965* 

5l6n locations 
o 
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SUBROUTINE DDIPLT 


LANGUAGE : 
PURPOSE: 


USE: 


FORTRAN 

To provide a one-call method of preparing plotting. This 
routine was originally designed for recording plots on 
the DD80 plotter only; however, it has heen. redesigned 
to use on any plotter. This one -call routine should not 
he used on any new ,1ohs ; new jobs requiring one-call dis- 
plays should use INFOPLT. 

These displays will not meet specifications for final 
figures . 

' CALL DDIPLT (IEC,IN„N,XDATA,YDATA,XMIN,XMAX,YMIN,YMAX, 

' NXM.,XM,NYM,YM,ISYMD} 

where ’ , 

. ^ r 

lEC is the code for terminating the frame 

• 0 frame incomplete 

, . 1 frame complete with this data. The frame 

change is built in and the plotter will be 



spaced for the next frame. 


IN 

is a two-word array. Each word contains 10 
Hollerith characters for plot identification. 

N. 

is the number of points to be plotted. 


XDATA 

is the name of the array containing the 
point values of X to be plotted. 

floating 

YDATA 

is the name of the array containing the 
point values of Y to be plotted. 

floating 

YMIN 

is the minimm value for X. 


.XMAX 

is the maximum value for X. 


YMIN 

is the minimum value for Y. 
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YMAX 


is the maximum for Y. 


The routine checks for the first call only to determine 
if either (XMAX-XMIN) or (YMAX-YMIN) is equal to zero. 

When either is zero', the routine will scan the X and/or 
Y array to determine the limits. For multiple curves 
per display, the limits must he specified on the first 
call to include all curves since the limits from the 
first call will be used for all curves . 

If any data falls outsidd the limits, it will be elimi- 
nated; but a count will be kept of all points dropped 
and written at top of the plot. 

Minimum/maximum values are next checked to see that the 
range is not zero. When it is, the specified values are 
adjusted by 10 percent of the minimum or set equal to 
il.O in cases where minimum and maximum are equal to zero. 

NXM is the number of central memory words in the 

message for the horizontal annotation. Maxi- 
mum number of words is 13j each word contains 
10 characters. If NXM and WYM are both neg- 
ative, tic marks will be generated instead of 
grid. 

XM is the name of array containing the label for 

the horizontal annotation. 

NYM is the number of words in the message for the 

vertical .annotation. Maximimi number of words 
is 13. 

YM ‘is the name of array containing the label for 

the vertical annotation. 

ISYMD is the integer code specifying the symbol or 
mode to be used for plotting the data values. 

1 Circle Q ^ Fan '9 

2 Square □ 9 Long diamond 0 

. 3 Diamond O 1° House Q 

U Triangls ^ 11 Circled dot Q 

5 Right Triangle 12 X 

6 Quadrant 13 Dot 

7 Dog House D l^t Vectors 
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RESTRICTIOlfS: 


METHOD: 


ACCURACY: 

REFERENCES: 

STORAGE: 

SUBPROGRAMS USED: 

OTHER CODING 
INFORMATION : 


The following arrays must he specified in a DIMENSION 
statement of the calling program: IN(2), XDATA(-N), YDATA(n), 

:XM(NXM), YM{NYM). 

.Each curve on a display requires a separate entry 'to 
the routine. X and Y coordinates for plotting must he in 
separate arrays of single precision, floating point data. 
Frame control is specified hy the lEC code in the calling 
sequence for the routine. 

Data are scaled and plotted; axes are drawn and annotated, 
and grid lines or tic marks are generated. 

Minimum /maximum values are adjusted to provide a range 
when all values of an array are equal. Adjustment is 
also made where needed to irrorove the appearance of the 
plot. 

Approximately three significant figures may he read in 
either direction. 


3021g locations 

CALPLT, NOTATE, NUMBER, PNTPLT, NFRAME 

A call to PSEUDO (l.k.l) must precede the first call 
to DDIPLT. An entry called VDIPLT with the same para- 
meters as DDIPLT is available which packs 8 6" xr 6" 
plots per frame for the Varian postprocessor. 



SUBROUTINE PSEUDO 


LMGUAQli 

purpose i 


USl! 


RESTRICTIONS: 


COMPASS 

To create and -write an appropriately named Plot Vector 
File. Through linkages set up by an initial call to 
PSEUDO > all subsequent graphics data generated by the 
user -will be routed through one of the PSEUDO entry 
points and written on the Plot Vector File. The PSEUDO 
processor is designed for use -with the frame dependent post- 
processors described in Section 1 . 5 , Volume IV, of the 
Computer Programing Manual. 

CALL PSEUDO 

or 

CALL PSEUDO (FN.) 

PN file name left-justified with zero fill. 

Default file name is SAVPLT. 


Example : 

CALL PSEUDO 

. 'This will establish a Plot Vector File named 
- SAVPLT. 

CALL PSEUDO (6LMYFILE-) 

This will establish a Plot Vector File named 
MYFI-LE. 


NOTE; The Plot Vector File (or Files) will usually be 
written to disk (as opposed to tape) and may be 
' • postprocessed following user program termination 

via appropriate specification of one or more 
PLOT control cards (see Section 1.3, Volume 
IV, Computer Programing Manual). 

(l) An initializing call to PSEUDO (with or without a 

file name argument) must be made prior to any calls 
to CALPLT or any other graphics -output routine. 



METHOD; 


ACCURACY: 

EEFERENCES: 

STORAGE: 

SUBPROGRAMS 


( 2 ) ’ Every Plot Vector File should be terminated 

with a 999 pen code, CALL CALPLT(0.0,0.0,999) • 

The transmission of the 999 code will cause an 
EOF write on the Plot Vector File, and the file 
will temporarily be closed. Thus, any given 
Plot Vector File will contain only one 999 pen 
code and/or one EOF. 

( 3 ) To continue plotting execution following trans- 
mission of a 999 code to a current Plot Vector 
File, the user program must call the PSEUDO pro- 
cessor to create new Plot Vector File (i.e., 

CALL PSEUDO ( 6LMYFIL2 ) ) . 

In addition to entry PSEUDO,, this processor contains two 
other entry points, namely PLT9999 an<3- PLT9998. An initial- 
izing call to PSEUDO will set PLT9999 into the processor 
switching mechanism (PLOTSW). Subseq,uent plot data gen- 
eration will then be routed via CALPLT, PLOTSW, and 
PLT9999 and written on the Plot Vector File. The entry 
PLT 9998 is used to record special purpose data from 
routines NFRAME and PLTSTOP, 


See Section 1 .^, Volume IV, Computer Programing Manual. 
215 5g locations total for direct subprograms 
USED: NUMARG, PLOTSW 
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APPENDIX C 


DEVELOPMENT OP PARAMETRIC MODEL EQUATIONS 
This appendix presents the development of the parametric 
model equations used in the computer program PARAM. 

The Larson-Miller j Orr~Sherby-Dorn, Manson-Succop, and 
Manson-Haf erd expressions are familiar time-temperature 
parameters. These parameters assume that the value of the 
parameter (a function of stress) is a constant for each 
value of the temperature compensated time parameter. The 
Rabotnov parameter (refs, 9 and 10) i& a time-stress 
parameter which assumes that the value of the parameter (a 
function of temperature) is a constant for each value of 
the time compensated stress parameter. Time to a given creep 
event and a polynomial in the parameter function (stress 
or temperature) were respectively the dependent and inde- 
pendent variables all regression model equation forms 
used in PARAM. The following presents the development of 

these five equation forms: 

■ ' 

Larson-Miller Parameter 
P = Tj^ (log f + C) = f (a)' : ‘ . 

Tj^(log t + C) = b^ + b2 log'a’+ b^dog cr)^ + b-i^Clog o)^ 

+ b^'dog o)^ + bgd'og cr)'^ 
assuming = -C • , . 

log t = b^ + b^/Tj^ + bg log cf/Tj^ + b^ (log ct)VTj^ 

+ b|j (log a)^/T^ + b^ (log o)^/Tj^ + bg (log ct)^/Tj^ 



where P = the Larson-Miller parameter 

t 

= temperature, °R ^ 

t = time to a particular creep event 

C = Larson-Miller constant 

* > / ”, 
t 

0 = applied stress 

t) _ _ t),' = coefficients estimated by method of 

o ' b ■ ■ 

least squares; 

■» 

Orr-Sherby-Dorn Parameter 
P = t exp (-AH/RTj^) = g (a) 
log t - K (AH/RTj^) = f (log a) 
assuming = KAH/R 

P 3 

log t = b^ +■ + i>2 log 0 + b^ (log o) + b^ (log 0)“^ 

+ b^ (log 0)^ + bg (log 0)^ 
where P = Grr-Sherby-Dorn parameter 

t = time to a particular creep event 
AH = apparent activation energy 
R = universal gas constant 
T,, = temperature, Kelvin 
0 = applied stress 

b bzr = coefficients estimated by method of 

o b 

least squares, 

Manson-Succop Parameter 
P = log t + CTp = f Co) 
log t = -C Tp + f (0) 
assuming b^^ = -C 
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2 3 

log t = + b^Tj, + ^2 log 0 + b^ Clog 0) + b-j^ (log 0) 

+ b^ (log 0)^ + bg (log 0)^ 
where? = Manson^Succop parameter 

t == time to a particular creep event 
C = Manson-Succop constant 
Tp = temperature, 

0 = applied, stress 

b — b/: = coefficients estimated by method of 
o o 

least squares, 

Manson-Haferd Parameter 
P = (log t - log t^)/(a?p - Tj^) = f(a) 
log t = log t^ + (Tp - 0?^) f(0) 
assuming b^ = log t^ 

O 

D = Tp - 

2 

log t = b^ + b^D + t>2 D log 0 + b^ D (log 0) 

+ b^ D (log 0)^ + b^ D (log 0)^ + bg D (log 0)^ 

where P = Manspn-Haferd parameter 

t = time to a particular creep event 

t = offset time 
a 

Tp temperature, 

= offset temperature, °P 
0 = applied stress 

b — b/r = coefficients estimated by method of 
o D 

least squares which iteratively searched values 
of to determine best fit. 
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Rabotnov Parameter 


P = a(l + At“) = fCT) 

= -1/A + 1/Aa [C^ + C^/T + C^/T^ + C^/T^ + C^/T^ 

+ C^/T^ + Cg/T^] 
assuming = -1/A 
b^ = C^/A 

t“ = bo + b^/aT + b2/aT^ + b^/oT^ + bj^/oT^ + b^/aT^ 

+ bg/aT^ 

where P = Rabotnov parameter 
a = applied stress 
A, a = constants 

t = time to a particular creep event 
T = temperature, 

bo — bg = coefficients estimated by method of least 
squares which iteratively searched values of a 
to determine best fit. 
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4 
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alloy 9 316 STAINLESS STEEL 


3142.90 
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25.00 
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- 213.00 
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28.00 

656.20 
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25.00 
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22.00 
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20.00 
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1200.00 

1 8.50 

15790.80 

1200.00 

17.00 

290.90 

1225.00 

25.00 

186.50 

1250.00 

25.00 

81.50 

1275.00 

25.00 

36.50 

1300.00 

25.00 

104.10 

1300.00 

22.00 

228.20 

1300.00 

20.00 

258 . 1 0 

1300.00 

19.00 

319.00 

1300.00 

18.00 

377 . 50 - 

1 300 .00 

17.00 

753.70 

1300.00 

16.00 

785.30 

’ i - 300.00 

16.50 

1232.50 

1300.00 

15.00 

185 . 4.60 

1 300-.00 

13.60 

■ 2421.00 

1 300.00 

13.00 

4078.30 

1300.00 

12.00 

62 ' 58 . 1 0 
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11.00 
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25.00 

’ 9.90 

1 350.00 

25.00 

2.70 

1400.00 

25.00 

, 83 '. 30 

1400.00 

15.00 

> 251 .20 

1400.00 

12.50 

921 ,-. 00 

/ 14 , 00 .. 00 

10.00 

' . 27.90 

1450.00 

15.00 

75.20 

r450.00 

12.50 

5. '00 

1 - 500.00 

16.40 

40.60 ' 

1500.00 

12.50 

87.90 

1500.00 

10.00 

1 70.40 

■ 1500 . 00 . 

9.00 

614.90 

1500*00 

7.00 

28 i 70 

■ 1 - 550.00 

10.00 

Figure 

!,<!- Input data 

for Case 1 
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L EAST-SQUARFS REGRESSION PROGRAM FOR PARAMETRIC 
ANALYSIS OF CREEP-RUPTURE DATA 


DATA SET 
OPTION CARO 
INPUT= 1 
inUT= 0 

TA= -5000 

RA= .2000 

CASF control cards 
PARAMETER CODE 


ALLOY 


316 STAINLESS iSTEEL 


NO. COEFFICIENTS 


L 

2 

3 

4 

5 


4 

4 

4 


Figure 2. -Output for Case 1. 



/ 


LEAST-SQLJAR6S REGRESSION PROGRAM FOR PARAMETRIC 
ANALYSIS OF creep-rupture DATA 


ON 

ru 


— regression values 

DATA SET ALLOY 

PARAMETER SELECTED 
,ND. OF OBSERVATIONS 
NO. OF INDEPENDENT VARIABLES 

residual oegre s of freedom 

F - VALUE 

RESIDUAL MEAN CQUARE 
STANDARD ERROR 
RESIDUAL >SUH OF SOUARFS 
TOTAL SUMS OF SOUARFS 
MULT. CORREL. COEF. SQUARED 


9 316 STAINLESS STE=L 

L-M 

38 
3 
3<f 
A76.8 
2. IA95E- 02 
1. A661E-01 
7. 3083E-01 
3.1A766F01 
.9768 


MIN 

Y = 4.31F 

-01 MAX Y 

» 4.20e<-00 

V 

RANGE = 3. 

77E«-00 MEAN Y = 2. 

45EFOO 




I 

0 

VARIABLE 

COEF.PJI ) 
-1.8792E»01 

S.E.CDEF. 

T 

MFAN X( ! ) 

HIN X( I > 

MAX XI 1 1 

RAN XI I ) 

RINF 

PSUM 

CERR 

1 

1/T 

4.7642E+04 

2.98EF03 16. 

00 

5.613E-04 

4. 98E-04 

6. 12E-04 

1 .14E-04 

1.44 

.335 

.00 

2 

S/T 

-3.7957E«-03 

4.14E+03 

92 

6.929E-04 

4.31E-04 

8190E-04 

4 .596-04 

-.46 

.663 

.00 

3 

S**2/T 

-3.29795F03 

1.69E*03 1. 

95 

8.6B9E-04 

3. 64E-04 

1.31E-03 

9.50E-04 

-.83 

.003 

-.00 


VARIABLE CODE 

, S=LOG STRESS 
T=TEHPFRATURE 
OT=T-TA 
L=STRESS 


95 PERCENT PREDICTION INTERVAL STATISTICS 




LOG TIME 

REAL TIME FACTOR (ANTILOG WIDTH) 

AVERAGE 

MIOTH 

.621 

4.2 

MAXIMUM 

WIDTH 

-T41 

5.5 


Figure 2 » -Continued, 
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LEAST-SOU?PES REPRESSION PRf)GRA>1 FOR PARAMETRIC 
ANALYSIS OF CREEP-RUPTURE DATA 


regrfssion values -- — 

DATA ‘SET ' -ALLOY. 

PARAMETER-.SELECTEO - - ■ 

NO. OF OBSERVATIONS. - 

NO. OF INOEPENDENT VARIABLES 

RESIDUAL DEGREES OF FREEDOM 

F r VALUE ' - 

RFSIDUAL M=AN SQUARE ^ 

STANDARD ERROR 
RESIDUAL SUM OF- SQUARES 
TOTAL SUMS OF SQUARES. 

HOLT. CORREL. ,COEF,. , squared. ‘ - 


9 316 STAINLESS, STEEL 

O-SrO ■' ■ 

. •‘ 38 - 

. 3 - . 

' 3A 

' " ■ -S75.9 ■ 
i.78 67E-02 ‘ 
i."3367E-0l ’ . 

6.07A8E-01 
3. 147.6Et-01 • 

.9807 


MIN Y 


A.31‘^-01 • MAX V = 4.26E«-00 Y RANGE 


3.77E*00 


MEAN Y 


2.45E+00 


1 
0 
1 

2 
3 


VARIABLE COEF.PUl 

• • , -1.5630Et0'l 

l/'T - 2illl7Et04 

S l.5077E*00 

S**2 -3. 3333E*-00 


S.E.COEF. T 

5.27F»0’2 40.07 
2.02FFO0 .75 
8.36E-01 3.99 


MEAN X( I \ 

1.0 llE-03 
l.228E*00 
1.533E+00 


MIN XI I ) 

8.96E-04 

B.ASEtOL 

7.14E-01 


MAX XI n 

UlOE-03 
1.48E+00 
2. 18E+00 


RAN X( J > 

2.06F-04 
6.32Er0l 
i .47EF00 


RINF. 

1.15 

.25 

-1.30 


PSUM 

.333 

.658 

.009 


CERR 

-.00 

.00 

.00 


VARIABLE CODE 

S=LOG STRESS 
T=TEHPSRATURE 
OT=T-TA 
' L=STRESS 


95 PERCENT PREDICTION INTERVAL STATISTICS 

LOG TIMP real TIME 


FACTOR lANTILOG WIDTH) 


AVERAGE WIDTH .567 

maximum w'iOTH .593 


3.7 

3.9 




Figure 2. -Continued 
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LEAST-SQUARES REGRESSION PRCGRAH POR PARAHETRIC 
ANALYSIS Of- CREEP-RUPTURE DATA 


REGRESSION VALUES 

DATA SET ALLOY 

PARAMETER SELECTED 

NO. OF OBSERVATIONS 

NO. OF INDEPENDENT VARIABLES 

RESIDUAL OEGREES'OF FREEDOM 

F - VALUE 

residual mean square 

STANDARD ERROR 
residual sum of SQUARES 
TOTAL 'SUMS OF SQUARES 
HULT. CORREL. COEF. SQUARED 


9 316 stainless steel 

H-S 

38 
3 
3A 
338.2 
3.0019E-02 
1. 7326E-01 
1.0206E*-00 
3.r476E + 01 
.9676 


min 

Y = 4.31E 

-01 MAX Y 

= 4.20EfOO 

Y RANGE = 3. 

77E<-00 1 

1EAN Y = 2. 

45E+00 




I 

VARIABLE 

COEF.PIU 

S.E.COEF. 

T MFAN X( I 1 

MIN X(I ) 

MAX XIII 

RAN X( I ) 

RINF 

PSUH 

CERR 

0 

1 

T 

2.2556E+01 
-1. L658E-02 

3.80E-O4 30 

.69 1.328E+03 

1.17E+03 

1, 55E + 03 

3.75E<-02 

-1.16 

.330 

.00 

2 

S 

-7.9643E-01 

2.64E+00 

.30 l.228E*O0 

8.45E-01 

1.4BE*-00 

6.32E-01 

-.13 

.665 

.00 

3 

S**2 

-2.3794E+00 

l.09E*00 2 

.19 1.533E+00 

7.14E-01 

2.18E + 00 

1.47E+00 

-.93 

.005 

-.00 


^ Q 
8 8 

S ^ 


VARIABLE CODE 

S=LOG STRESS 
T=TEMP'^RATURE 
OT=T-TA 
L=STRESS 


95 PERCENT PREDICTION INTERVAL STATISTICS 

LOG TIME REAL TIME 


FACTOR < antilog WIDTH) 


AVERAGE WIDTH .735 ' 

MAXIMUM WIDTH .771 


Figure 2 .-Continued 



LFAST-SQUftPES REGRESSION PROGRAM FOR PARAMETRIC 
ANALYSIS OF CREEP-RUPTURE DATA 


REGRFSSION VALUES 

DATA SET ALLOY 

PARAMETFR SELECTED 

NO. OF OBSERVATIONS 

NO. OF INDEPENDENT VARIABLES 

RESIDUAL OEGREFS OF FREEDOM 

F - VALUE 

RESIDUAL MEAN SQUARE 
STANDARD ERROR 
RESIDUAL SUM OF SQUARES 
TOTAL SUMS OF SQUARES 
MULT. CORR>=L. COEF. SQUARED 


9 316 STAINLESS STEEL 

M-H 

38 

3 

3A 

A19.8 

2. A33AE-02 
1. 5599E-01 
8.2737'E-OI 

3. 1476E>01 

.9737 


MANSCN - HAFERC CONSTANTITAJ = 


300.0 


os 

MIN Y = 4.31E 

-01 MAX Y 

= 4.20Et00 

Y RANGE = 3. 

77E+00 MEAN Y = 2. 

45E*00 




or 

I VARIABLE 

0 

COEF.PlI I 
1.4509E tOl 

S.E.CQEF. T 

MFAN X(I) 

MIN X(I ) 

MAX x( n 

RAN X(I ) 

RINF 

PSUH 

CERR 


1 C" 

-1.0686E-02 

I.26E-03 8.48 

l.O28E<-03 

8i75E«-02 

1.25EF03 

3.75E+02 

.-1.06 

.328 

-.00 


2 OT*S 

4.8648E-03 

2. 15E-03 2.26 

1.25lE«-03 

l.OlE+03 

1.54E4-03 

5.24E*02 

• 6 8 

.651 

— . 00 


3 OT*S«'*? 

-4.6243E-03 

9.02E-04 5.13 

1.548E*03 

8.57EF02 

2. 15E+03 

1.29E*03 

-1.59 

.021 

-.00 


VARIABLE CODE 

S=LOG STRESS 
T=TFMPERATURE 
DT=T-TA 
L=STRESS 


95 


PERCENT PREDICTION INTERVAL STATISTICS 

LOG TIME real TIME FACTOR (ANT 


ILQG WIDTH) 


AVERAGE WIDTH .661 

MAXIMUM WIDTH .805 


Figure 2. -Continued. 
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LGAST-SQUAPES REGRESSION PROGRAM FOR PARAMETRIC 
ANALYSIS OF CREEP-RUPTURE DATA 


REGRESSION VALUES 

DATA SET ALLOY 

PARAMETER SELECTED 

NO. OF OBSERVATIONS 

NO. OF INDEPENDENT VARIABLES 

RESIDUAL DEGRE'^S OF FREEDOM 

F - VALUE 

RESIDUAL MEAN SQUARE 
STANDARD ERROR ' 

RESIDUAL SUM OF SQUARES 
TOTAL SUMS OF SQUARES 
MULT. CORREL. COEF. SQUARED 


S 316 STAINLESS STEEL 
RAB 


38 

2 

35 ' 

220. A 

2.1137E-03 . 
4.5975E-02 
7.3979E-02 
1.0058EF00 

.9269 , 


RABOTNQV CONSTANT (RA) = .05630 


HIN Y = 1.06EtOO MAX Y = 1.72EF00 " Y RANGE = 6 


I 

D 

VARIABLE 

COEF.PU ) 
8.7599E-01 

S.E.COEF. 

T 

MEAN XII) 

I 

1/L»T 

-3.6886EF09 

2.06£«-03 

17.99 

9.693E-05 

2 

1/L*T«2 

6.3869E+07 

3. 19E+06 

20.01 

3.505E-08 


66E-01 MEAN Y = 1.38E*00 


HIN XII 1 

MAX X< 1) 

RAN XU) 

RINF 

PSUH 

CERR 

2.78E-05 

9.52E-05 

6.75E-05 

-3.79 

.092 

0.00 

2.09E-08 

6.35E-08 

9.31E-08 

9.13 

.908 

0.00 


g| 

tag 

Ss 

- 

3 I— I 
^ CO 


VAB IA0LE CODE 

S=LOG STRESS 
T=TEHPERATURE 
DT=T-TA 
L=STRESS 


95 PERCENT PREDICTION INTERVAL STATISTICS 

LOG TIME REAL TIME FACTOR (ANTILQG WIDTH) 


AVERAGE WIDTH 1.088 

MAXIMUM WIDTH 1.350 


12.3 

22.9 


Figure 2, -Concluded 



-5000*0 


0.2 


0 O 
2 3 
2 4 
2 5 
2 6 
2 7 


ALLOY 9 

316 STAINLESS 

STEEL 

3142*90 

1 175.00 

25.00 

74,60 

1 200,00 

30.00 

213.00 

1 200.00 

28,00 

656.20 

1200.00 

25.00 

3476.10 

1200.00 

22.00 

6825.30 

1200.00 

20.00 

1 0076.50 

1200.00 

18.50 

1 5790.80 

1200.00 

17.00 

290,90 

1225.00 

25.00 

186.50 

1 250.00 

25.00 

81,50 

1275.00 

25.00 

36,50 

1 300.00 

25.00 

104.10 

1 300.00 

22,00 

228,20 

1 300.00 

20.00 

258,10 

1300,00 

19,00 

319.00 

1300.00 

18.00 

377.50 

1 300.00 

17.00 

753.70 

1300.00 

16.00 

785.30 

1300.00 

16.50 

1 232.50 

1 300*00 

15.00 

1854,60 

1 300.00 

13.60 

2421 ,00 

1 300.00 

1 3.00 

4078,30 

1300.00 

12.00 

6258,10 

1 300.00 

11.00 

21 ,50 

1 325.00 

25.00 

Q,90 

1 350.00 

25.00 

2.70 

1400.00 

25.00 

83.30 

1400.00 

1 5,00 

251 .20 

1 400.00 

12.50 

921 ,00 

1400.00 

10.00 

27.90 

I 450.00 

15,00 

75,20 

1450.00 

12,50 

5.00 

1500.00 

16.40 

40.60 

1500.00 

12.50 

87.90 

1500.00 

10.00 

170,40 

1500.00 

9.00 

614.90 

1500.00 

7.00 

28.70 

1550.00 

1 0, 00 


Figure 3 .“Input data for Case 2, 
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LEAST-SOUABES REGRESSION PROGRAM FOR PARAHET.RIC 
ANALYSIS OF .CREEP-RUPTURE DATA 


o\ 

OO 


REGRESSION VALUES — t-t-t 

data set alloy 

OARAMETER SELECTED 
NO. OF OBSERVATIONS 38 

NO. OF INDEPENDENT VARIABLES ' 2 

RFSIOUAL DEGREES OF FREEDOM 35 

F ^ VALUE ■ . 600.1 

»=S1'DUAL MEAN SQUARE '2.5A81E-02 

STANDARD ERROR . " li5963E-01 

RESIDUAL SUM OF SQUARES , 8.9183E-01 

TOTAL SUMS OF SQUARES 3.1A76Et01 

MULT. CORREL. fOEF. SQUARED .9717 


316' STA INteSS STEFL 
0-S-D 


MIN y = 4.31E-01 


MAX Y = 9.20E+00 


I VARIABLE 

0 . 

1 l/T 

2 S 


COEF.P(I) 
-1. 1091E+01 
2.133SE +09 
-6.5279E «-00 


S.E.CQEF. 

6.26EF02 
2. 31E-01 


Y.’RANGE = 3.77E*qO 

MEAN'XI I ) 


mean y. = 


2.45E*00 


34.09 

28.22 


l.0nE.-:O3 

1.228E»00 


MIN X'(K 

MAX X( I » 

RAN X( I ) 

RINF 

PSUM 

8.96E-04 

I. ioe-03 

2-06E-04 

1.16 

.336 

8:45E-01 

l.'48E«-00 

6.32E-01 

-I.IO 

• 664 


VARIABLE CODE 

S=LDG STRESS 

t=tempfrature 

OT=T-TA 

L=STRESS 

95 PERCENT PREDICTION INTERVAL STATISTICS 

LOG TIME REAL TIME FACTOR (ANTILOG WIDTH) 


AVERAGE WIDTH 
MAXIMUM WIDTH 


.66 8 

.703 


4.7 

5.0 


CERR, 

0.00 

.00 


Figure i|. -Continued 



LEASTi-S QUAKES REGRESSION PROGRAM FOR PARAMETRIC 
•ANALYSIS OF CREEP-RUPTURE DATA 


— r . REGRESSION VALUES 

..PATA SET ' ' ^ ALLOY 

PARAMETER SELECTED 
NO. OF OBSERVATIONS 
. NO. OF - INDEPENDENT VARIABLES 
RESIDUAL DEGRE = S OF FREEDOM 
' F - VALUE 

RESIDUAL MEAN SQUARE 
STANDARD, ERROR 
RESIDUAL SUM OF SQUARES 
■"TOTAL ' SUMS OF SQUARES 
• - MOLT. CDRREL. 'COFF-. SQUARED 


9 316 STAINLFSS STE^L 

O-S-D 
38 ■ 

3 
3A 
575.9 
1..7867E-02 
1. 3367E-01 
6.07A8E-01 
3. 1A76E«-01 
.9807 


os 

*• 

HTN Y = 4.-3 IE 

-01 MAX Y 

= 4.20E«-00 

Y 

range = 3- 

77E400 MEAN Y = 2. 

45E400 




VO 


r VARIABLE 

COEF.PH I 

S.E.COEF." 

T 

MEAN X( I I 

MIN X( I 1 

MAX X( I ) 

RAN X( t ) 

RINF 

PSUH 

C-ERR 



0 «• 
1 l/T 

-1.5630E *01 
2. illTE^OA 

5.27E<-02 40 

.07 

l.OllE-03 

8.96E-04 

l.lOE-03 

2.06E-04 

1.15 

.333 

-.00 



2 S 

i.5077E*00 

2.02E*00 

.75 

1.22BF+00 

8.45E-01 

1.48Ef00 

6 .32E-01 

•?5 

.6 58 

• 00 



3 S’»»2 

-3.3333E»00 

8.36F-01 3 

.99 

1.533 = 4-00 

7.14E-01 

2- 1BE4-00 

1 .47E4-00 

-1.30 

.009 

♦ 00 


variable CODE 

S=LOG STRESS 
T=TFMPRRATURE 
DT=T-TA 
L=STRFSS 


05 


PERCENT PREDICTICN INTERVAL STATISTICS 

log time real TIME FACTOR (ANT 


ILOG WIDTH) 


AVERAGE WIDTH .567 

MAXIMUM WIDTH .593 







Figure -Continued 



UeAST-SQUAR£S REGRESSION PROGRAM FOR PARAMETRIC 
ANALYSIS OF CREEP-RUPTURE DATA 


REGRESSION VALUES 

DATA SEr ALLOY 9 316 STAINLESS STEEL 


PARAMETER SELECTED 


O-S-0 








NO. OF OBS = RVATIONS 


33 








NO. OF INDEPENDENT 

VARIABLES 

4 








RESIDUAL DEGREES OF 

FREEDOM 

33 








F - VALUE 


580.4 








RESIDUAL MEAN SQUARE 

1.33 67E-02 








STANDARD ERROR 


1. 1562F-01 








RESIDUAL SUM OF SQUARES 

4.4113F-01 








TOTAL SUMS OF SQUARES 

3. 1476EF01 








MULT. CriRREL. COEF. 

SQUARED 

.9860 








MIN Y = -V.BIE 

-01 MAX Y 

= A. 20=400 V 

RANGE = 3. 

77E*00 M=AN Y = 2. 

45E40Q 




1 VARIABLE 

COEF.P ( I > 

S.F.COEF. T 

MEAN Xt I ) 

MIN X(I 1 

MAX X(I) 

RAN XM ) 

RINF 

PSUM 

CERR 

0 

6.7051E+00 









1 1/T 

2.0982E 4-04 

4. 57E4 02 45.87 

1.0 llE-03 

8.96E-04 

l.lOE-03 

2.06E-04 

1.14 

.331 

.00 

2 S 

-5.7228E+01 

1.67E401 3.42 

1 .228E400 

8.45S-01 

1. 48E400 

6.32E-01 

-9.60 

.654 

-.00 

3 S**2 

4.7587E401 

1.45E401 3.29 

1.533E400 

7.14E-01 

2.18E400 

1.47E 400 

18,58 

.009 

-.00 

4 S**3 

-1.4563E401 

4.13F400 3.53 

1.942E+00 

6.04E-01 

3. 22E400 

Z.62E400 

-10.13 

.0 05 

-.00 


VARIABLE CODE 

S=LOG STRESS 
T=TEMPERATURE 
DT=T-TA 
L=STRESS 

95 PERCENT PREDICTION INTERVAL STATISTICS 

LOG TIME REAL TIME FACTOR (ANTILOG WIDTH) 


AVERAGE WIDTH .A96 

MAXIMUM WIDTH .569 


3.1 

3.7 


Figure 4, -Continued 


ORIGINAL PAGE IS 
OF POOR QUALiry 



i/’l 


LEAST-SQUftPFS REGRESSION PROGRAM FDR PARAMETRIC 
ANALYSIS OF CREEP-RUPTURE DATA 


~<I 

H 


REGRESSION VALUFS 


DATR set A,UL OY 9 316 STAINLESS STEEL 


PARAMETER SELECTED 

, 

0-S 

-D 








NO. OF OBSERVATIONS 



38 








NO. OF INDEPENDENT 

VARIABLES 

. 

5 








RESIDUAL DEGREES OF 

FREEDOM 

- 

32 .. 








F - VALUE 


594 

-a 








RESIDUAL MEAN SOUAR 

c 

1.0472E- 

02 








STANDARD ERROR 


1. 02.33E- 

01 








RESIDUAL SUM OF SQUARES - 

3.3509E- 

01 








TOTAL SUMS OF SOUAR 

= 5 

3.1476ET01 








MULT. CORREL. COEF, 

SQUARED 

-9894 








MIN Y = ^k.31E 

-01 max Y 

= 4.20E+00 Y 

RANGE = 3. 

77ET00 H 

EAN Y = 2. 

45ETOO 




T VARIABLE 

COEF-PI! ) 

S.E.COEF. 

T 

MEAN X( I ) 

MIN X(I ) 

max X( I) 

RAN X(I ) 

RINF 

PSUM 

CERR 

0 

1 l/T 

-l.0736E*02 
2. 1161E T04 

4.09E+02 

51.77 

l.OIlE-03 

8.96E-04 

1. IOE-03 

2.06E-04 

1.15 

.330 

-.00 

2 S 

3.4852E+02 

1.28E+02 

2.72 

1.228ET00 

8.45E-01 

1.48EF00 

6.32E-01 

58.47 

.6 52 

— .00 

3 S**2 

-4.87B0E+02 

1.69FT02 

2.89 

1.533E+00 

7.14E-01 

2.1BET00 

1.476*00 

-190.05 

.009 

-.00 


2.9592ET02 

9. 76E+01 

3.03 

1.942E+00 

6.04E-01 

3.22ETOO 

2.62E*00 

205.77 

.005 

-.00 

5 

'6.6787E+0I 

2.10E+01 

3.18 

2 .492F+00 

5. lQE-01 

4.76ET00 

4.25E*00 

-75.36 

.003 

-.00 


VARIABLE CODE 

S=LOG STRESS 
T=TEMPERATURE 
DT=T-TA 
L=STRPSS 

95 PERCENT PREDICTION INTERVAL STATISTICS 

log time real time FACTOR (ANTILOG WIDTH) 

AVERAGE WIDTH .444 2-8 

MAXIMUM WIDTH .560 3.8 


Figure 4, -Continued, 
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f=4ST-S0UARES REGRESSION PROGRAM FOR PARAMETRIC 
ANALYSIS OF CREEP-RUPTURE DATA 



REGRESSION VALUES 

DATA SET ALLOY 

PARAMETER SELECTED 

NO. OF OBSERVATIONS 

NO. OF INDEPENDENT VARIABLES 

5'SIPUAL DEGREES OF FREEDOM 

F - VALUE 

’ESIDUAL MEAN SQUARE 
STANDARD frrhR 
=FSIDUAL SUM OF SQUARES 
TOTAL SUMS OF SQUARES 
MULT. CORR-L. COEF. SQUARED 


9 316 STAINLESS STE=L 

0-S-D 
38 
6 
31 
ABA. 8 
1. 07086-02 
1.03A8E-01 
3.319AE-01 
3. IA75E»01 
.9895 


« A.31S 

-01 MAX Y 

= A.20E+00 

V 

RANGE = 3. 

77FTOO mean Y = 2. 

A5E+00 




V&'^TABLE 

COEF. PU ) 

S.E.COEF. 

T 

MEAN XI t 1 

MIN XU ) 

MAX XUI 

RAN X( 1 > 

RINF 

PSUH 

CERR 

; / T 

A.7690ET01 

2.1181E+0A 

A.15E*02 51 

.03 

l.OllE-03 

8.95E-0A 

1. lOE-03 

2.06E-0A 

1.16 

.330 

.00 

s 

-3.3766E*02 

1.27E4-03 

.27 

1.228ET-00 

8.ASE-01 

1.A8E+00 

6.32E-01 

—56.65 

.652 

• 01 

S**2 

7. 1A82E «-02 

2.22EV03 

.32 

1.533ET00 

7.1AE-01 

Z.18E4-00 

1.47E«-00 

278.50 

.009 

• 00 

S**3 

-7.A798EF02 

1.93EF03 

.39 

1.9A26 + 00 

< OAE-01 

3.22E*00 

2.62EF00 

-520.10 

• 005 

• OO 


3.8218E4-02 

e. 28 E «^02 

.A6 

2.A92ET00 

5.106-01 

A.76EF00 

A.25EFOO 

A31.2A 

.003 

.00 

S**5 

-7.6572E»0l 

1.AIE«^02 

.5A 

3.235E*00 

A.31E-01 

7.03E + 00 

6.60E*00 

-13A.18 

.000 

• 00 


VARIABLE CODE 

S=LOG STRESS 
T=TEMPFRATURE 
OT=T-TA 
L=STRESS 


95 


PERCENT PREDICTION INTERVAL STATISTICS 

LOG TIME real TIME 


FACTOR (ANTILDG WIDTH) 


average width .ASA 

MAXIMUM WIDTH .590 


2.8 

3.9 




Figure 4. -Concluded 



3 

2 


3 

6 


- 5000*0 


0.2 


ALLOY 9 

316 STAINLESS 

STEEL 

31 42.90 

1 175.00 

25.00 

74.60 

1200.00 

30.00 

213.00 

1200.00 

28.00 

656.20 

1200*00 

25.00 

3476. 10 

1200*00 

22.00 

6825.30 

1200.00 

20.00 

1 0076.50 

1200.00 

18.50 

1 5790.80 

1200.00 

17.00 

290.90 

1225.00 

25.00 

186 .50 

1250.00 

25.00 

81 .50 

1275.00 

25.00 

36.50 

1300.00 

25.00 

' 104.10 

1300.00 

22.00 

228.20 

1300.00 

20.00 

258. 1 0 

1300.00 

19.00 

319.00 

1300.00 

18.00 

377.50 

1 300*00 

17.00 

753.70 • 

1300.00 

16.00 

785.30 

1300.00 

16.50 

1232.50 

1300.00 

15.00 

1854.60 

1 300.00 

13.60 

2421.00 

1 300.00 

13.00 

4078.30 

1 300.00 

12.00 

6258 . 10 

1 300.00 

11.00 

21.50 

1 325.00 

25.00 

9.90 

1 350.00 

25.00 

2.70 

1400.00 

25.00 

83.30 

1400.00 

15.00 

251.20 

1400.00 

12.50 

921 .00 

1400.00 

1 0.00 

27.90 

1450.00 

15.00 

75.20 

1 450.00 

12.50 

5.00 

1500.00 

16*40 

40.60 

1 500.00 

12.50 

87.90 

1500*00 

10.00 

170.40 

1 500.00 

9.00 

614.90 

1500.00 

7 . 00 

28.70 

1550.00 

1 0.00 


Figure 5. -Input data for Case 3. 
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LEAST-SQUARES REGRESSION PROGRAM FOR PARAMETRIC 
ANALYSIS OF CREEP-RUPTURE DATA 
DATA SET ALLOY 9 316 STAINLESS STEEL 

OPTION CARD 
INPUT= 3 
IQUT= 3 

TA= -5000 

RA= ,2000 

CASE CONTROL CARDS 

PARAMETER CODE NO. COEFFICIENTS 

2 6 


Figure 6, -Output for Case 3 



INPUT 

DATA OBSERVATIONS 

NO. 

TIME 

1 

3142.90 

2 

74.6 0 

3 

213.00 

4 

656.20 

i:: 

3476. 10 

6 

6825.30 

7 

10076, 50 

8 

15790.80 

0 

290.90 

10 

186..50 

IL 

81.50 

! 2 

36.50 

13 

104.10 

14 

228.20 

15 

258, 10 

16 

319.00 

17 

377.50 

18. 

753.70 

19 

785.30 


1232.50 

21 

L 854. 6 0 

22 

2421.00 

23 

4078.30 

24 

6258.10 

25 

21.50 

26 

■ . 9.90 

27 . 

2.70 

28 

83.30 

29 

251.20 

30 

921.00 

31 

27.90 

3? 

75.20 

33 

5.00 

34 

40.60 

35 

87.90 

36 

170.40 

37 

614.90 

38. 

28.70 


STRESS 

TEMPERATURE 

25 

1175 

30 

1200 

28 

1200 

25 

1200 

22 

1200 

20 

1200 

18 

12 00 

17 

1200 

25 

1225 

25 

1250 

25 

1275 

25 

1300 

22 

1300 

20 

1300 

19 

1300 

18 

1300 

17 

1300 

16 

1300 

16 

1300 

15 

1300 

14 

1300 

13 

1300 

12 

1300 

11 

1300 

25 

1325 

25 

1350 

25 

1400 

15 

1400 

13 

1400 

io 

1400 

15 

1450 

13 

1450 

16 

1500 

13 

1500 

10 

1500 

9 

1500 

7 

1500 

10 

1550 


Figure 6. -Continued* 
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/S3. 



lpast-souabes regression program for parametric 

ANALYSTS OF CREEP-RUPTURE DATA 


regression values 

DATA SET ALLOY 

PARAMETER SELECTED 

NO. OF OBSERVATIONS 

NO. OF INOEPENPENT VARIABLES 

RESIDUAL DEGREES OF FREEDOM 

F - VALUE 

RF^IOUAL «=AN SQUARE 
STANDARD ERROR 
“FSIOUAL SUM OF SQUARES 
InTAL SUMS OF SQUARES 
MULT. CORREL. COEF. SQUARED 

MIN Y = 4.31E-01 MAX Y = 


? 316 STAINLESS STEEL 

O-S-0 
33 
5 • 

32 
59A.8 
1.04T2E-02 
1.0233E-01 
3. 3509E-01 
3. 1^76E+01 
.989'» 


4.20F+00 


Y RANGE = 3.77E*-00 


mean Y = 


2.45E«00 


P c 


§ ^ 
gty 


FIRST 5. nBSTRVATIOMS - TRANSFORMED VARIABLES 


3.49733=»00 
1.87274= *09 
2. 32838= *-00 
2-81704r*00 
3.54109-*-00 


XI - XtL21 


1.10132E-03 

1.08473E-03 

1.38473E-03 

1.08473E-03 

1.08473E-03 


1.39794E+00 

1.47712F*00 

1.44716E*00 

1.39794E*00 

1.34242E*00 


1.95424E*00 
2, i8189ET0a 
2.09427E+00 
1.9S424£*00 
1.80210E+00 


2.73191E»00 

3.2229lE*-00 

3.03073E*00 

2.73191ET00 

2.41918E«-00 


3.81904E4-0a 

4.76063E»00 

4.3a59SE«-00 

3.81904E*^00 

3.24756E*00 


VARIABLE 

COEF.PU I 
-1.0736E*02 

S.E.COEF. 

T 

MEAN XIII 

MIN XIII. 

MAX XIII 

RAN XI I 1 

RINF 

PSUH 

CERft 

1/T 

2. 1161E*04 

4.09E*02 

51.77 

l.OllE-03 

8.96E-04 

1. lOE-03 

2.06E-04 

1.15 

.330 

-.00 

S 

3.4852= *02 

1.28E*02 

2.72 

1.228E*00 

8.45E-0L 

1.48E*00 

6.32E-01 

58.47 

.652 

-.00 

S*’!=2 

-4.8780E*02 

l.69E*02 

2.89 

l.533F*00 

7.14E-01 

2. 18E*00 

1.47E*00 

-190.05 

.009 

-.00 

S**3 

2.9592= *02 

9. 76 E* 01 

3.03 

1.942E*00 

6.04E-01 

3.22E+00 

2.62E*00 

205.77 

.005 

-.00 

S**4 

-6.6787E*0i 

2.10F*01 

3.18 

2.492F*00 

5. lOE-01 

4.76E*00 

4.25E*0O 

-75.36 

.003 

-.00 


VARIABLE CODE 

S=LOG STRESS 
T=TEHPCRATURE 
DT=T-TA 
L=STRESS 

95 PERCENT PREDICTION INTERVAL STATISTICS 

log TIME REAL TIME FACTOR < ANTILOG WIDTH! 


AVERAGE HIOTH 
maximum width 


Figure 6. -Continued 



T^S'^ 


0=SinUALS - REGRESSION SPACE 


ns 

Y OB' 

Y CALC 

RESIDUAL 

PCT5RR 

1 

3.4975«-00 

3. 266=+00 

-2.3125-01 

-6.5 

?. 

1.873F+00 

1.879C+00 

5.4055-03 

.3 

1 

2.328= <-00 

2.331F+00 

2.687E-03 

.1 

4 

2.817E400 

2.915E+00 

9.799E-02 

3.5 

r 

3.541F*00 

3.402F+00 

-1.391E-01 

-3.9 

6 

3.S34=+00 

3.688F+00 

-1.462F-01 

-3.8 

1 

4. 003F400 

3.B92F+00 

-1. 109E -01 

-2.8 

8 

4.198F400 

4.0984+00 

-1. 004E-0 1 

-2.4 

c. 

2. 46 4f*00 

2.5744+00 

1, 106F-01 

4.S 

;o 

2.271E«-00 

2.244F+00 

-2. 704E-02 

-1.2 

11 

1.911F«-00 

1.922E+00 

1. 1305-02 

.6 

’ 2 

1. 562F*00 

1.610F+00 

4.811E-02 

3.1 

13 

2. -iWFtOO 

2.0974+00 

7.9905-02 

4.0 

14 

2.358F<-00 

2.383F+00 

2.501E-02 

1.1 

1? 

2.41 2= ♦90 

2.5204*00 

1.082E-01 

4.5 

lb 

2. 50 4= + 00 

2.656F+00 

1.519E-01 

6.1 
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